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SUMMARY 


Ve  describe  the  development  of  the  ZNOVAETC  (X5TC)  Code,  e  model  of 
interior  ballistic  phenomena  baaed  on  a  numerical  solution  of  the  governing 
equations  for  one-dimensional,  multi-phase  flow,  XETC  is  an  extension  of 
the  previously  developed  XNOVAX  Code.  The  extensions  include  revisions  to 
the  representation  of  reactive  sidewalls  such  as  combustible  case  elements, 
the  modeling  of  monolithic  charges,  the  analysis  of  charge  increments  bonded 
to  the  tnbe  or  the  projectile,  a  representation  of  a  ballistic  control 
device  intended  to  reduce  the  temperature  coefficient  of  the  charge,  and  the 
incorporation  of  logic  to  treat  end-burning  traveling  charge  increments. 
Moreover,  the  tank  gun  and  traveling  charge  features  have  been  fully  linked 
to  the  chemistry  models. 

The  XXTC  Code  is  applied  to  the  simulation  of  traveling  charges  with 
finite  reaction  2ones  to  assess  the  extent  to  which  the  ballistic 
perfoxmance  benefit  of  the  traveling  charge  is  degraded  as  the  reaction  zone 
thickness  increases.  It  is  concluded  that  reaction  zone  thicknesses  of 
several  calibers  can  be  tol  era  ted  without  a  significant  loss  of  perfomance. 
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1.0  iNTXODUCnON 


The  par  pose  of  this  report  is  to  do  case  nt  the  steps  taken  to  create  the 
XNOVAKTC  (XKTC)  Code  from  other  recent  versions  of  the  NOVA  Code.  The  XKTC 
Code  is  intended  to  provide  digital  simulations  of  the  interior  ballistics 
of  a  vide  range  of  gun  propelling  charges.  Like  all  veraions  of  the  NOVA 
Code,"  XKTC  is  based  on  a  numerical  solution  of  the  governing  equations  for 
the  macroscopic,  qua  si- one -dimensional  flow  defined  by  the  solid  propellant 
and  its  products  of  combustion.  The  XKTC  Code  has  been  developed  us  an 
amalgam  of  the  previously  developed  XNOVAT^  and  NGVAXC^  Codes  which 
respectively  address  details  of  tank  gun  and  traveling  charges.  Our 
intention  in  this  report  is  to  describe  certain  additional  features  which 
have  been  encoded  into  XKTC  and  to  provide  a  completely  updated  description 
of  the  use  of  the  code.  However,  we  do  not  provide  a  complete  description 
of  the  governing  equations  or  the  method  of  solution. 

This  introduction  contains  two  sections.  In  Section  1.1  we  provide 
some  background  information  concerning  the  various  versions  of  the  NOVA  Code 
which  have  led  to  the  development  of  XKTC,  In  Section  1.2  we  summarize  the 
new  features  and  cross-linkages  which  are  particular  to  XKTC.  Analysis 
pertinent  to  the  new  features  is  provided  in  Chapter  2.0.  In  Chapter  3.0  we 
describe  the  application  of  XKTC  to  the  traveling  charge i  we  investigate  the 
extent  to  which  the  ballistic  benefits  of  the  end-burning  traveling  charge 
would  be  compromised  by  a  reaction  zone  of  finite  thickness.  In  the 
Appendix  we  provide  a  complete  description  of  the  use  of  the  XKTC  Code. 

1.1  Background  Information 


Our  intention  in  this  section  is  simply  to  clarify  the  nomenclature  for 
the  various  versions  of  the  NOVA  Code  without  going  into  the  detailed 
differences  between  them.  The  earlier  reports  cited  here  may  be  consulted 
for  further  discussion.  The  NOVA  Code  was  originally  developed  to  provide  a 
means  of  analyzing  the  aspects  of  charge  design  which  contribute  to  the 
formation  of  longitudinal  pressure  waves  in  the  chamber  of  a  gun.  Although 
several  earlier  versions  were  developed,  we  aider  stand  the  NOVA  Cede  to  be 
defined  by  the  version  described  in  Reference  1. 


1  G  P.  S.  "The  NOVA  Code:  A  User's  Manual* 

Indian  Head  Contract  Report  IHCR  80-8  1980 

^  Gough,  P.  S.  'XNOVAT  -  A  Two-Phase  Flow  Model  of  Tank  Gun  Interior 

Ballistics"  Final  Report,  Task  Order  I,  Contract  DAAK11-85- 0-0002  1985 

^  Gough,  P.  S.  "A  Two-Phase  Model  of  the  Interior  Ballistics  of  Hybrid 
Solid-Propellant  Traveling  Charges” 

Final  Report,  Task  I,  Contract  DAAK11-82-C-0154  1983 
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Briefly,  the  NOVA  Code  was  bated  on  the  balance  eqnations  for  macro- 
scopically  one-dimensional  tr ~-ph* ic  flow.  The  state  variables  were  to  be 
thought  of  as  averages  of  the  local  valnes  or  microproper ties.  Intractable 
aiicroflow  details  such  as  drag,  heat  transfer  and  propellant  combustion  were 
assumed  to  be  related  to  the  macroscopic  variables  by  means  of  empirical 
correlations.  The  NOVA  Code  allowed  the  simulation  of  a  broad  class  of 
conventional  charges  consisting  of  grannlar  or  stick  propellant  arranged  in 
several  increments.  The  governing  equations  were  solved  by  the  method  of 
finite  differences  with  an  explicit  allowance  for  the  discontinuities  in  the 
state  variables  at  the  internal  boundaries  defined  by  the  ends  of  the 
increments. 

The  JNOVA  Code^  was  developed  to  take  advantage  of  more  efficient 
computational  procedures  which  had  been  established  during  work  on  a  two- 
dimensional  interior  ballistics  code.6  Prom  a  modeling  standpoint  INOVA 
retained  most  of  the  features  of  NWA,  only  certain  esoteric  and  seldom  used 
options  being  deleted  to  produce  a  compact  code.  However,  INOVA  also 
contained  a  modeling  extension  relative  to  NOVA  in  that  a  dual-voidage 
representation  of  perforated  stick  charges  was  admitted  according  to  which 
interstitial  propertie  s  were  distinguished  from  those  within  the 
perforations. 

The  XNOVAX  Code6  was  an  extension  of  INOVA  in  which  the  products  of 
combustion  of  the  propellant  ..nd  igniter  were  permitted  to  react  chemically. 
Whereas  earlier  code  versions  had  always  assumed  combustion  to  proceed  to 
completion  locally  and  a  haul taneonsly  with  regression  of  the  surface  of  the 
burning  propellant,  INO'/AK  adopted  the  viewpoint  that  the  products  of 
combustion  consisted  of  a  homogeneous  mixture  of  gases,  droplets  and 
particles  in  which  a  number  of  chemical  reaction*  could  occur. 

XNOYAK  was  itself  extended,  as  described  in  the  previous  task  report, 
to  become  XNOVAT.^  The  INOVAT  Code  incorporated  numerous  features  pertinent 
to  the  modeling  of  tank  gun  propelling  charges,  including  case  combustion 
and  projectile  afterbody  intrusion.  While  the  chemistry  options  of  XNOVAK 
were  retained,  they  were  not  linked  to  the  new  tank  gun  options  of  INOVAT. 


*  Gough,  P.  S.  "INOVA  -  An  Express  Version  of  the  NOVA  Code" 
Final  Report  Contract  N00174-82-M-8048 

6  Gough,  P,  S,  "Modeling  of  Rigidized  Gun  Propelling  Charges" 
Contract  Report  ARBKL-CR-0051 8 

6  Govgh,  P.  S.  "Theoretical  Modeling  of  Navy  Propelling  Charges" 
Final  Report,  Contract  N00174-83-C-0241  FGA-TR-84-1 
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Prior  to  the  development  of  XNOVAK,  the  INOVA  Code  was  used  to  oreate 
the  NOV ATC  Code  which  added  to  the  features  of  XNOVA  the  possibility  of 
modeling  all  or  part  of  the  charge  aa  an  end-burning  traveling  charge. 
Combustion  of  the  traveling  charge  waa  treated  consistently  with  that  of  the 
conventional  propellant.  kegre salon  of  the  rear  face  of  the  traveling 
charge  was  assumed  to  yield  final  products  of  combustion  at  an  infinitesimal 
distance  fro*  the  surface. 

1.2  Objectives  and  Suaaary  of  Results 


The  objectives  of  the  present  effort  have  been  two-fold.  First,  we 
have  formed  a  new  code  by  the  replacement  of  INOVA  by  XNOVAT  in  NOV  ATC. 
Second,  we  have  added  certain  new  features  and  encoded  cross-linkages  of  the 
various  options  to  produce  the  code  which  we  refer  to  as  the  ZNOVACIC  Code, 
or  XKTC.  Third,  we  have  used  XKTC  to  explore  the  extent  to  which  the 
ballistic  performance  of  a  traveling  charge  would  be  compromised  if  the 
formation  of  final  products  of  combustion  were  completed  over  a  finite 
length,  rather  than  at  an  infinitesimal  distance  from  the  base  of  the 
traveling  charge,  as  assumed  in  NOV  ATC. 

Figure  1.1  illustrates  three  types  of  propelling  charge  which  can  be 
modeled  by  XKTC.  Figure  1.1  (a)  represents  a  typical  multi- increment  charge 
of  the  type  originally  addressed  by  NOVA  or  XNOVA.  Each  increment  may 
consist  of  granular  or  stick  propellant.  Unslotted  perforated  stick 
propellant  is  given  a  dual-voidage  representation.  Figure  1.1  (b) 
represents  a  multi-increment  tank  gun  charge.  The  projectile  afterbody  is 
allowed  to  intrude  into  the  region  occupied  by  the  charge  and  reactive 
sidewall  components  are  admitted.  It  is  also  possible  to  model  the  presence 
of  increment  endwalls  as  reactive  layers  which  resist  penetration  by  the 
combustion  products.  The  increments  may  also  be  described  as  parallel 
packaged  with  appropriate  formulations  of  the  flow  resistance  and  heat 
transfer  corrections.  Figure  1.1  (c)  represents  a  multi-increment  charge 
in  which  some  of  the  increments  burn  in  a  traveling  charge  mode,  in 
successive  planar  layers  from  the  rear.  Reactive  sidewalls  are  also 
admitted,  in  the  region  occupied  by  the  conventional  increments. 

An  effort  has  been  made  to  link  all  the  code  options  in  s  physically 
complete  manner.  However,  it  is  assumed  that  there  are  no  compactible 
filler  materials  present  if  the  afterbody  intrudes  into  the  chamber  or  if 
the  traveling  charge  option  is  exercised.  It  is  also  assumed  that  the 
projectile  does  not  have  an  afterbody  if  the  end-burning  traveling  charge 
option  is  used. 

Apart  from  transfering  the  traveling  charge  model  from  NOV  ATC  to 
XNOVAT,  a  number  of  revisions  and  extensions  were  added  in  the  development 
of  XKTC.  These  are  described  in  full  in  Chapter  2.0.  They  may  be 
summarized  as  follows. 
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(•)  Artillery  or  Navy  Cate  Gun  Charge 


(b)  Tank  Gun  Charge 


(O  Charge  Co.pri.iag  End  Burning  Traveling  Charge  Increment. 


Figure  1.1  Charge  Configuration.  Repre.ented  by  XETC  Code. 
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First,  we  have  extended  the  representation  of  reactive  sidewalls  to 
admit  a  variation  of  thermochemical  and  aechanioal  properties  with  axial 
potition.  Second,  we  have  encoded  a  fora  function  for  a  nonolithic  charge 
which  is  bonded  to  the  tube  and  barns  only  on  the  snrface  of  a  single 
internal  perforation.  The  projectile  afterbody  is  peraitted  to  intrude  into 
the  perforation.  Third,  we  have  encoded  logic  to  represent  any  charge 
increment  as  bonded  either  to  the  tube  ox  to  the  projectile.  This  feature 
may  be  used  to  describe  a  traveling  charge  increaent  which  burns  in  depth 
rather  than  at  the  rear  surface,  and  also  admits  intrusion  of  the  projectile 
afterbody.  Fourth,  we  have  encoded  logic  to  describe  a  ballistic  control 
device  whose  purpose  is  to  reduce  the  temperature  coefficient  of  the  charge 
through  the  use  of  a  separately  burned  sub-charge.  Finally,  we  have  linked 
the  chemistry  options  to  the  sidewall  and  endwall  reactivity  aodels  and  to 
the  combustion  of  the  traveling  charge. 
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2.0  REVISIONS  AND  EXTENSIONS  TO  EQUATIONS 


As  discussed  in  Chapter  1.0,  most  of  the  governing  equations  for  the 
XKTC  Code  have  been  documented  in  previous  reports  and  it  is  not  an 
objective  of  the  present  report  to  provide  a  comprehensive  statement  of  all 
the  model  details.  However,  we  do  discuss  those  equations  and  linkages 
which  are  new.  This  chapter  contains  four  sections.  In  Section  2.1  we 
discuss  the  revised  representation  of  the  reactive  sidewalls  and  enthralls. 

In  Section  2.2  we  discuss  the  representation  of  the  monolithic  charge.  In 
Section  2.3  we  discuss  the  treatment  of  a  charge  increment  which  is  bonded 
either  to  the  projectile  or  to  the  tube  of  the  gun.  Finally,  in  Section  2.4 
we  discuss  the  analysis  of  a  ballistic  control  device  which  is  intended  to 
reduce  the  temperature  coefficient  of  the  charge. 

2.1  Revised  Representation  of  Reactive  Sidewalls  and  Endwalls 

In  the  previous  report, ^  which  described  the  development  of  XNWAT,  we 
discussed  the  representation  of  reactive  sidewalls  and  endwalls.  The 
sidewalls  were  intended  to  represent  combustible  case  components  and/or 
ignition  elements  and  were  understood  to  be  attributes  of  any  or  all  of  the 
following:  the  tube*  the  centerlinei  the  projectile  afterbody.  The 
sidewalls  were  characterized  by  local  value  a  of  thickness  and  surface 
regression  rate.  Both  ignition  and  compressibility  were  taken  into  account 
and  provision  was  made  for  a  layer  of  deterrent.  However,  each  sidewall  was 
considered  to  have  the  same  mechanical  and  the nso chemical  properties  over 
its  entire  length.  The  sidewall  on  the  tube  was  permitted  to  have  different 
properties  from  that  on  the  centerline  or  the  projectile  afterbody  hut  it 
was  not  considered  to  consist  of  a  number  of  segments  of  differing 
propertie  s. 

In  XKTC  the  representation  of  the  sidewalls  has  been  revised  so  that 
each  sidewall  may  be  characterized  as  consisting  of  up  to  three  segments  as 
shown  schematically  in  Figure  2.1.  All  the  mechanical  and  themochemical 
properties  may  vary  from  segment  to  segment.  The  initial  thickness  of  the 
layer  remains  an  arbitrary  function  of  position.  However,  discontinuities 
in  thickness  are  not  recognized  explicitly  by  the  numerical  method  of 
solution.  We  also  do  not  track  the  segment  boundarie s  with  precision.  The 
properties  of  the  sidewall  at  each  mesh  location  are  those  of  the  segment  in 
which  the  mesh  point  lies  and  no  attempt  is  made  to  average  sidewall 
properties  when  the  mesh  point  is  close  to  a  segment  boundary. 

An  additional  modif  ication  in  2ETC  is  concerned  with  the  treatment  of 
heat  transfer  to  the  tube.  At.  a  tube  wall  location  which  is  covered  by  the 
sidewall,  the  heat  transfer  is  assumed  to  be  zero  until  the  sidewall  is 
completely  burned  through. 


We  have  also  completed  the  linkage  of  sidewall  reactivity  to  the 
chemistry  submodels  of  XNOVAK.^  In  the  previous  report^  we  had  noted  the 
governing  equation  for  the  rate  of  change  of  T^»  the  mass  fraction  of 
species  i  in  the  form 
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(2.1.1) 


where  Yflj£,  YSei  and  Ysii  are  the  mass  fractions  of  species  i  in  the 
products  of  combustion  of  the  igniter,  the  outer  sidewall  and  the  inner 
sidewall  respectively;  Yij.o  ^  the  mass  fraction  of  species  i  produced  by 
the  near  field  (fizz)  reaction  of  propellant  jj  and  rik  is  the  rate  of 
production  of  species  i  by  reaction  k.  Ve  also  have  e,  porosity*  p, 
density;  <p,  mSe>  ms^  and  mj  as  rates  of  production  per  unit  volume  of 
igniter  products,  outer  sidewall  products,  inner  sidewall  products  and  near 
field  products  of  propellant  j  respectively.  Finally,  w ±  is  the  rate  of 
loss  of  species  i  due  to  deposition  on  the  surface  of  the  solid  propellant. 

In  the  present  code,  the  values  of  Y#ci  and  Y«ii  are  fully  supported 
when  the  user  elects  to  exercise  both  the  tank  gun  and  chemistry  options  at 
the  same  time. 

A  similar  extension  applies  to  the  reactive  enthralls.  Whereas  the 
analysis  of  the  endwalls  was  previously  unlinked  to  the  chemistry  option, 
XiCl'C  requires  that  the  composition  of  the  products  of  combustion  of  each  of 
the  substrates  be  specified  when  the  chemistry  option  is  in  effect.  The 
internal  boundary  conditions  are  then  solved  subject  to  the  additional 
balance  laws 
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where  mx  is  the  mass  flax  from  the  mixture  .region  at  the  boundary  point 
inside  the  endwall,  a,  is  the  mass  flux  at  the  boundary  point  outside  the 
endwall,  us^  is  the  rate  of  reaction  of  the  j-th  sublayer,  and  Yj^,  Yi#, 

A  are  mass  fractions  of  species  i  corresponding  to  i,  and 
respectively. 

2.2  Representation  of  Monolithic  Charge 


In  Figure  2.2  we  illustrate  a  propelling  charge  which  caaprises  a 
monolithic  increment.  Such  an  increment  is  assumed  to  be  bonded  to  the  tube 
and  inhibited  on  its  end  surfaces  so  that  oombustion  is  confined  to  the 
surface  of  the  single  central  perforation.  Moreover,  we  consider  the 
possibility  that  the  afterbody  of  the  projectile  may  penetrate  the 
perforation  of  the  monolithic  increment. 

Ve  represent  the  monolithic  charge  in  ZKTC  as  a  single  voidage  stick 
bonded  to  the  tube.  Since  end  burning  of  stick  propellant  is  neglected  in 
the  code,  the  inhibition  of  the  ends  is  automatically  captured.  However,  it 
is  necessary  to  take  care  with  the  definition  of  the  porosity  and  the  for* 
functions  in  order  to  model  properly  the  rate  of  heat  transfer  daring  the 
ignition  phase  and  the  subsequent  rate  of  pressurization  due  to  combustion. 

Let  d0  be  the  initial  diameter  of  the  perforation  and  let  d  be  the 

local  surface  regression.  Ve  assume  that  d0  is  a  function  of  position,  as 

suggested  by  Figure  2.2.  Let  At  and  A*  respectively  denote  the  cross- 

sectional  areas  of  the  tube  and  the  afterbody,  corrected  for  the  presence  of 
any  reactive  sidewalls.  The  cross- sectional  area  for  the  two-phase  flow  is 
considered  to  be 


A  =  Aj  —  Aa  •  (2.2.1) 

Let  Ap  =  n/4  (d0  +  2d)1  be  the  area  of  the  perforation.  Then  the  porosity, 
or  fraction  of  the  flow  cross-section  occupied  by  the  products  of 
combustion,  is 


(2.2.2) 


It  follows  that  if  we  define 


S_  =  n  <d~  +  2d) 
v 


(2.2.3) 
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Outside  conforms  to  tube  geometry*  is  bonded 
and  does  not  burn* 


10 


Figure  2.2  Representation  of  Single  Perforation  Monolithic  Charge  by 


(2.2.4) 


Vp  =  A  -  (Ap  -  A.) 

then  the  ratio  Sp/Vp  will  yield  the  correct  rate  of  pressurization  of  the 

perforation  when  combustion  occurs.  However,  the  heat  transfer  to  the  solid 
propellant  daring  the  ignition  phase  must  be  based  on  the  hydraulic  diameter 


4(A  -  A  ) 

D  =  - 2 - L_  (2.2.5) 

p  (S  +  S  ) 

P  * 


where  Sp  is  given  by  (2.2.3)  and  Sa  is  the  circumference  of  the  afterbody. 

We  have  already  commented  on  the  fact  that  we  treat  the  monolithic 
grain  as  bonded  to  the  tube.  We  discuss  the  analysis  of  bonded  grains  in 
the  n'xt  section  where  we  also  consider  the  possibility  of  bond  rapture. 
However,  rupture  of  the  bond  which  attaches  the  monolithic  grain  is 
presently  assumed  never  to  occur.  If  we  do  consider  rapture  we  have  also  to 
consider  the  possibility  that  the  outer  surface  of  the  grain  will  no  longer 
conform  with  the  inner  surface  of  the  tnbe.  This  will  affect  the 
c.r  eolation  of  the  form  functions  and,  more  importantly,  will  lead  to  a 
consideration  of  the  formation  of  an  outer  annular  ullage  region  and  the 
'issibility  of  ignition  of  the  outside  of  the  grain.  Since  combustion  of 
the  outer  surface  could  well  result  in  the  consideration  of  mass  addition  in 
a  region  of  arbitrarily  small  flow  cross- sect  ion,  we  have  regarded  this 
topic  as  being  beyond  the  scope  of  the  present  effort. 

If  a  reactive  layer  is  attached  to  the  tube  wall  in  the  region  occupied 
by  the  consolidated  charge,  it  is  assumed  to  be  thermally  insulated  until 
the  charge  has  completely  burned  through. 

2.3  Analysis  of  Bonded  Charge  Increment 


Figure  2.3  illustrates  a  charge  increment  which  is  bonded  to  the 
projectile.  Rather  than  writing  a  momentum  equation  for  the  projectile  and 
bonded  charge  as  a  system,  we  consider  the  equations  of  motion  for  each  and 
introduce  an  explicit  force  of  bonding  which  ensures  that  they  remain  in 
contact.  This  approach  is  computationally  convenient  as  it  allows  us  to 
make  use  of  existing  coding  structures  with  only  minor  modifications.  Also, 
the  explicit  computation  of  the  bonding  force  allows  us  to  consider  a 
rupture  criterion  according  to  which  the  charge  may  separate  from  the 
projectile.  We  first  consider  a  charge  bonded  to  the  projectile  as  in 
Figure  2.3.  Subsequently,  we  comment  on  the  case  when  the  increment  is 
bonded  to  the  tube  of  the  gun. 
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Hie  equation  of  motion  of  the  projectile  it 
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(2.3.1) 


where  Mp  is  the  mass  of  the  projectile*  vp  is  the  projectile  velocity,  pa 

and  a t  are  the  pressure  and  intergranular  stress  at  the  base  of  the 

projectile  and  Aa  is  the  flow  cross-sectional  area  at  the  sane  point  and  is 

defined  by  Equation  (2.2.1).  Ve  assume  that  the  charge  is  bonded  to  the 

afterbody  over  the  interval  [  zx ,  2j]  and  that  t  is  the  bonding  force.  The 

symbol  d  denotes  the  force  due  to  gas  pressure  and  intergranular  stress  over 
that  part  of  the  afterbody  to  the  rear  of  the  bonded  increment.  F  is  the 
bore  resistance.  We  use  the  constant  g0  to  reconcile  nnits  of  measurement. 
As  in  the  previous  section,  A#  is  the  cross-sectional  area  of  the  afterbody. 

Ihe  equation  of  motion  of  the  solid  propellant  is 

Du 

(1  -  e)p  — 2.  +  g  (l  -  e)  »  — -  f  --  (2.3.2) 

P  Dt  °  3  z  °  3  x.  SA 
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where  pp  is  the  density  of  the  propellant,  up  is  the  velocity,  fg  i*  the 
interphase  drag  and  D/Dtp  convective  derivative  along  the  propellant 

streamline. 

Now  consider  auxilliary  variables  vp  and  tp  such  that 


dv  *  *  dA 

M  — —  —  g(p  +  o’  )  A  +  g  j(p'^o)  —  dz  +  d  —  F  ,  (2.3.3) 
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Ve  may  multiply  (2.3.4)  by  A,  integrate  over  a  and  add  to  (2.3.3)  and 
perfc  a  similar  operations  on  (2.3.1)  and  (2.3.2)  to  obtain  the  physically 
expected  result 


dv  *• 

Do 

dv  ** 

Du 

M  — 2-  +  (1  -  e)o  A  — 2  dz  = 

U  — 2  +  ( 1  -  a )  p  A 

—P  dz 

P  dt 

P  Dt 

P  dt  J  P 

Dt 

P 

P 

Since  the  condition  of  bonding  requires  v„  *  Up  ft  follows  that  the  updated 
quantities  obey 

vp+l  K +  s a  -  e)pPAdz | =  vrx  +  j  (i  -  e)*>pK+ld2  •  (2-3-6) 

The  computational  algorithm  therefore  reqnires  that  we  first  integrate 
(2.3.3)  and  (2.3.4)  to  get  vj}+1  end  p+l.  Then  (2.3,6)  yields  Vp+X  end 
hence  u“+  * .  We  have  assumed  thus  far  that  the  bond  between  the  propellant 
and  the  projectile  does  not  rupture.  The  force  of  bonding  may  be  determined 

as 


F  = 


gAt 


M 

P 


(2.3,7) 


where  At  =  At  on  the  predictor  step  and  At/2  on  the  corrector  step  of  the 
finite  difference  integration.  Separation  of  the  propellant  from  the 
projectile  is  assumed  to  occur  if  F  exceeds  a  predetermined  value. 

The  intergranular  stress  a  follows  from  the  usual  constitutive  1  aw 
according  to  which  it  is  an  irreversible  function  of  porosity.  Since  axial 
strain  cannot  occur  for  the  bonded  charge,  the  stress  will  be  controlled  by 
combustion  of  the  propellant  and  variations  in  tube  area  with  travel.  The 
boundary  values  of  o  are  likewise  determined  from  the  constitutive  law  and 
not  from  the  characteristic  forms. 
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The  analysis  of  an  increment  bonded  to  the  tube  is  analogous  except 
that  an  obvious  simplification  arises  since  the  tube  is  assumed  to  remain 
stationary.  In  place  of  (2.3.7)  we  evaluate  the  bonding  force  from 


f  u“+  p  A(1  -  e )dz 
J  p  p 

F  =  - -  (2.3.8) 

*oAt 


We  also  note  that  the  attribute  of  attaclnent  to  the  projectile  or  to  the 
tube  applies  to  the  increment  as  a  whole.  If  the  increment  consists  of  a 
mixture  of  propellants  or  of  several  parallel  packaged  bundles,  all  the 
species  or  bundles  are  taken  to  be  bonded  until  the  rupture  condition  is 
achieved. 

2.4  Representation  of  Balllccic  Control  Device 

In  Figure  2.4  we  illustrate  a  ballistic  control  device  who-.e  intended 
purpose  is  the  reduction  of  the  temperature  coefficient  of  the  propelling 
charge.  The  small  control  chary..  is  burned  prior  to  or  during  the  ignition 
of  the  main  charge.  Since  thrust  is  supplied  to  the  projectile  via  the  bate 
of  the  afterbody,  the  control  charge  has  the  effect  of  varying  the  position 
of  the  projectile  during  or  prior  to  the  ignition  of  the  main  charge.  At 
higher  temperatures  the  displacement  of  the  projectile,  at  the  time  of 
ignition  of  the  main  charge,  is  expected  to  increase  thereby  offsetting  the 
increased  quickness  of  the  main  charge  and  reducing  the  temperature 
coefficient. 

In  our  schematic  illustration  we  show  features  of  the  ZJCI'C  representa¬ 
tion  which  may  or  may  not  be  present  in  actual  designs.  We  show  a  combus¬ 
tion  chamber  within  the  device  whose  diameter  differs  from  that  of  the 
propulsion  tube  into  which  the  projectile  afterbody  intrudes.  We  also  show 
sidevents  along  the  device  throngh  which  an  ignition  stimulus  to  the 
propelling  charge  may  be  induced  prior  to  the  uncorking  of  the  afterbody. 

The  XXTC  Code  also  allows  the  exterior  of  the  device  to  have  an  arbitrary 
shape . 

The  model  oi  the  ballistic  control  device  includes  the  following 
details.  Conditions  vithin  the  device  are  presmed  to  be  uniform  since  the 
device  is  not  expected  to  be  very  long  end  the  resolution  of  axial  structure 
according  to  a  continuum  model  does  not  seem  worthwhile.  The  governing 
equations  for  the  state  of  the  gas  within  the  device  are  therefore 
statements  of  the  balance  of  mass  and  energy  supported  by  a  burn  rate  law 
for  the  control  charge  and  a  covolume  equation  of  state.  The  control  charge 


15 


is  assumed  to  be  granular  and  to  consist  of  one  of  the  following  fonts: 
sphere;  cylinder;  monoperf orated  with  or  without  outside  inhibition;  seven- 
perforation.  The  control  charge  is  assumed  to  be  ignited  after  a 
predetermined  delay. 

The  equation  of  motion  of  the  projectile  is  modified  to  take  into 
account  the  thrust  due  to  the  control  charge  until  the  instant  when  the 
projectile  uncorks. 

The  external  geometry  of  the  device  is  used  to  correct  the  values  of 
the  cross-sectional  area  of  the  two-phase  flow  defined  by  the  main  charge. 
For  the  present  we  assume  that  the  device  is  short  enough  that  the  vented 
gases  can  be  coupled  to  the  reponse  of  the  main  charge  through  representa~ 
tion  as  a  basepad  attached  to  the  breechface.  The  rate  of  venting  per  unit 
vent  area  is  deduced  from  the  pressures  within  the  control  device  and  at  the 
rear  of  the  main  charge  according  to  an  isentropic  flow  law  with  an 
allowance  for  choking.  The  total  rate  of  flow  follows  from  the  total  vent 
area  which  is  exposed  as  the  projectile  moves.  This  total  area  consists  of 
the  sidewall  ntribution  and.  when  the  projectile  uncorks,  the  area  of  the 
propulsion  tube. 

The  total  flux  so  computed  is  used  to  construct  a  surface  source  term 
which  is  expressed  as  an  attribute  of  the  breechface.  Reversed  flow  into 
the  control  chamber  is  not  presently  considered.  In  subsequent  work  it  is 
intended  to  provide  the  option  of  representing  the  flax  from  the  control 
device  as  a  sidewall  source  attributed  to  the  internal  boundary  of  the  two- 
phase  flow.  This  will  permit  the  representation  of  longer  devices  than  the 
present  method. 

The  governing  equations  for  the  state  of  the  gas  within  the  device  are 
easily  derived  and  we  simply  summarize  them  here.  The  balances  of  mass  and 
energy  are 


V  m 

*  dt  p 


in 


“  pA_v 
r  c  p 


de 

dt 


m  ( e 
P  P 


e ) 


-  pA  v 
c  p 


(2.4.1) 


(2.4.2) 


where  p  is  density;  p,  pressure;  e,  internal  energy;  Vg,  volume  available  to 
gas;  Ac>  cross-sectional  area  of  propulsion  tube;  vp,  projectile  velocity! 
mp,  rate  of  combustion  of  control  charge;  pp,  density  of  control  charge; 
ep,  chemical  energy  of  control  charge;  me>  total  masa  flux  to  main  charge, 
assumed  to  be  always  positive  or  zero. 
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The  volume  available  to  the  (tea  ii  given  by 


M 

V  =  V  +  A  2  -  -£■  (2.4.3) 

g  o  c  p 

PP 

where  V0  is  the  volume  of  the  conbnstios  chamber;  zp  i*  the  displacement  of 
the  base  of  the  afterbody  relative  to  the  entrance  to  the  propulsion  tube; 
Bp  is  the  mass  of  the  unburned  propellant.  We  note  that  XETC  admits  an 
initial  gas  volume  VgQ  defined  by 

v80  -  vo  *  Vp.  (2-4-4> 

where  zpo  is  an  initial  standoff  distance.  The  initial  volume  may  be 
defined  through  either  Vc.  zpo  or  both. 

The  rate  of  combustion  of  the  propellant  is 


where  Sp  is  the  surface  area  of  a  grain;  VPo  the  initial  volume  of  a 
grain;  Mpo  is  the  initial  mass  of  the  charge  and  d  is  the  surface  regression 
rate  and  is  assumed  to  obey  the  usual  exponential  dependence  on  pressure. 

We  al  so  have 


M 

P 


(2.4.6) 


where  Vp  i$  the  current  volume  of  a  grain.  The 
related  to  the  total  surface  regression  through 
functions. 


values  of  Vp  and  Sp  are 
the  usual  geometrical  form 
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The  functional  dependence  of  i»e  on  tlie  pressures  within  the  control 
device  and  at  the  breech  of  the  nain  chaaber  ia  given  by  the  iaentropic  flow 
law  corrected  for  covolnme  aa  in  Reference  1. 

Provision  is  Bade  for  a  deterrent  layer  in  the  control  charge  but  there 
is  no  present  linkage  to  the  chemistry  options. 
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3 .0  EFFECT  OF  FINITE  FLAMS  THICKNESS  ON  TRAVELING  CHARGE  PERFORMANCE 


In  conventional  propelling  charges  the  propellant  tends  to  be 
distributed  in  a  nearly  uniform  fashion  over  the  length  of  the  tube  and  the 
velocity  distributions  of  both  the  propellant  grains  and  the  products  of 
combustion  tend  to  be  nearly  linear  functions  of  axial  position  as  suggested 
by  Lagrange.^  The  kinetic  energy  of  the  propulsion  gas  is  proportional  to 
that  of  the  projectile  and  represents  a  loss  of  ballistic  efficiency.  For 
artillery  weapons  operating  with  a  muzzle  velocity  of  approximately  1  km/sec 
the  ratio  of  propellant  mass  to  projectile  mass  (C/M)  is  about  0.2  and  the 
loss  is  not  very  important.  However,  there  is  a  current  interest  in  weapons 
operating  at  a  muzzle  velocity  of  3  km/  sec.  Estimates  of  the  required  value 
of  C/M  to  achieve  such  velocities  range  from  3  to  8  and  the  kinetic  energy 
of  the  propulsion  gas  therefore  represents  a  significant  loss. 

The  end-burning  traveling  charge  has  been  proposed  as  a  propulsion 
scheme  whereby  the  loss  due  to  the  kinetic  energy  of  the  propulsion  gas  may 
be  reduced. ^  A  model  of  the  traveling  charge  has  been  developed^  and  the 
theoretical  advantages  of  this  scheme  have  been  demonstrated.  However, 
theory  has  incorporated  the  assumption  that  the  final  products  of  combustion 
are  formed  an  infinitesimal  distance  from  the  regressing  rear  surface  of  the 
traveling  charge.  This  assumption  has  not  been  supported  by  experimental 
studies  of  those  f  ormul  ations  which  presently  show  the  most  promise. 
Combustion  has  been  observed  to  involve  a  complex  series  of  steps  which  are 
strongly  dependent  on  composition,  confinement  and  density  of  the  sample. 
Rather  than  consisting  of  a  region  of  unburned  propellant  separated  from  a 
region  of  final  combustion  products  by  a  thin  reaction  zone,  the  combustion 
process  was  seen  to  involve  a  preliminary  penetration  and  partial 
consumption  of  the  entire  sample  by  a  convective  flamefront  which  was  then 
followed  by  a  relatively  homogeneous  consumption  of  the  remainder  of  the 
propellant  accompanied,  in  some  cases,  by  deconsolidation. 
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The  objective  of  the  present  study  is  to  investigate  the  ballistic 
consequences  of  a  finite  reaction  zone  at  the  base  of  the  traveling  charge. 
We  do  cot  attempt  to  model  directly  the  complex  phenomena  reported  by  White 
et  al.*®  The  theoretical  model  is  sufficiently  broad  that  it  does  offer  the 
prospect  of  future  simulations  of  combustion  mechanisms  of  the  type 
described  by  these  authors.  In  the  present  study,  however,  the  model  is 
simply  exercised  to  describe  a  two-step  combustion  process  in  which 
regression  of  the  base  of  the  traveling  charge  yields  a  mixture  of  final 
(TC-F)  and  intermediate  (TC-I)  combustion  products.  The  intemediate 
products  react  to  completion  at  a  finite  rate  over  an  extended  region.  The 
thickness  of  the  combustion  zone  is  varied  by  varying  the  ratio  of  final  to 
intermediate  products  formed  in  the  first  step  and  the  rate  of  reaction  of 
the  intermediate  products. 

We  provide  a  summary  of  the  physical  content  of  the  model  in  this 
introduction.  The  governing  equations  are  summarized  in  Section  2.1. 

Before  discussing  the  model  we  comment  further  on  the  differences  between 
conventional  and  traveling  charges. 

We  have  already  noted  that  the  kinetic  energy  of  the  propellant  and  its 
products  of  combustion  represent  a  ballistic  loss  whose  importance  increases 
with  increa sing  muzz le  velocity.  The  original  concept  of  the  traveling 
charge  seems  to  be  due  to  Langweiler^  who  proposed  the  development  of  an 
end-burning  charge  attached  to  the  projectile  base  with  a  burn  rate  designed 
to  yield  products  of  combustion  at  rest  relative  to  the  tube.  Apart  from 
the  purely  technological  problem  of  producing  a  propellant  with  the 
necessarily  enormous  burn  rates  and  the  required  me chanical  strength,  it  was 
observed  by  Vinti^  that  the  proposed  burn  rates  would,  in  general,  require 
the  development  of  a  strong  deflagration  wave,  one  for  which  the  products  of 
combustion  would  be  supersonic  relative  to  the  flame  front.  The  strong 
deflagration  wave  is  believed  to  be  thermodynamically  unstabl e^  and 
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therefore  to  be  incapable  of  existing  in  a  steady  flow.  Although  the 
traveling  charge  is  burned  in  an  inherently  unsteady  manner,  the  strong 
deflagration  limit  should  nevertheless  apply  provided  that  the  combustion 
zone  is  sufficiently  thin  that  the  rates  of  change  of  mass,  momentum  and 
energy  within  it  remain  negligible  by  compari son  with  the  fluxes  of  these 
quantities  through  its  bounding  surfaces.  We  also  note  that  even  if  the 
strong  defl  agration  were  achievable,  it  would  not  necessarily  represent  a 
useful  state  from  an  engineering  standpoint  due  to  the  magnitude  of  the 
concomitant  pressure  drop  across  the  reaction  zone.  At  the  theoretical 
limit  of  sonic  or  choked  combustion  the  pressure  on  the  unreacted  side  of 
the  flame  is  approximately  twice  that  on  the  reacted  side.  The  ratio  of 
pressures  increases  indefinitely  with  Mach  number  and,  for  the  Langweiler 
proposal  at  least,  it  implies  indefinitely  increasing  stress  on  the  barrel 
throughout  the  combustion  of  the  charge. 

The  foregoing  objections  are  particular  to  the  Langweiler  concept. 

They  do  not  rule  out  the  possibility  of  improved  ballistic  performance 
through  a  more  general  traveling  charge  concept  in  which  the  rate  of  burning 
is  simply  required  to  be  great  enough  to  induce  a  substantial  rearward 
blowing  of  the  products  of  combustion. 

An  additional  phenemenon  to  be  considered  is  the  rarefaction  formed  at 
the  instant  the  traveling  charge  burns  out.  The  local  pressure  drop  may  be 
so  large  that  there  is  no  further  significant  propulsion  of  the  projectile 
following  burnout.  The  projectile  may  even  decelerate  due  to  the  resistive 
forces.  The  rarefaction  also  has  the  result  that  the  velocity  distribution 
of  the  combustion  products  is  relaxed  to  the  conventional  linear  Lagrange 
distribution  with  the  result  that  the  kinetic  energy  of  the  propellant  gas 
is  restored  to  the  conventional  value  and  the  benefit  of  the  traveling 
charge  is  apparently  lost.  It  may  be  expected  therefore  that  optimium 
traveling  charge  performance  will  involve  burnout  timed  to  occur  just  prior 
to  muzzle  exit. 

Although  the  initial  motivation  for  the  traveling  charge  appears  to 
have  stemmed  from  a  consideration  of  the  velocity  field  of  tho  propulsion 
gas.  it  is  our  view  that  attention  is  better  directed  towards  the  pressure 
distribution.  Associated  w  ith  the  linear  Lagrange  velocity  distribution  is 
a  parabolic  pressure  di stribution  who se  gradient  serves  to  accelerate  the 
propulsion  gas  down  the  tube.  The  pressure  at  the  projectile  base  is  less 
than  that  at  the  breech.  Accordingly,  propulsion  of  the  projectile  is  due 
to  a  lower  pressure  than  that  which  the  tube  must  withstand.  We  illustrate 
the  conventional  charge  in  Figure  3.1.  We  show  the  mixtt.re  region  separated 
fxca  the  projectile  base  by  a  small  region  of  ullage  —  anally  no  more  than 
one  or  two  calibers  —  which  is  due  to  the  inability  of  the  propellant 
grains  to  match  exactly  the  projectile  velocity.  We  also  sketch  the 
pressure  distribution.  According  to  the  appoximate  theory  of  Lagrange,  the 
ratio  of  breech  to  base  pressure  is  given  by  1  +  C/2M,  where  C  is  the  charge 
mass  and  M  is  the  projectile  mass.  If  C/M  =  0.2  as  is  typically  the  case 
for  artillery  weapons  firing  at  1  km/sec  then  the  ratio  of  pressures  is  1.1 
and  the  loss  of  efficiency  is  small.  On  the  other  hand,  if  a  value  of 
C/M  =  8  is  required  to  achieve  velociti>  of  the  order  of  3  km/sec,  then  the 
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ratio  becomes  5  and  it  is  clear  that  tho  projectile  is  receiving  little 
propulsive  benefit  from  the  pressure  exerted  on  the  tube.  It  should  he 
said,  however,  that  the  Lagrange  distribution  may  provide  a  very  poor 
characterization  of  the  pressure  field  in  a  conventional  gun  when  C/M 
exceeds  unity  and  the  pressure  drop  may  not  in  fact  be  as  large  as  1  +  C/2H, 
particularly  in  the  earlier  stages  of  the  propulsion  cycle. 

In  Figure  3.2  ve  illustrate  the  situation  for  an  ideal  end-burning 
traveling  charge.  Ve  assume  that  the  traveling  charge  is  ignited  following 
the  complete  combustion  of  a  booster  charge.  Therefore  the  unreacted 
propellant  is  separated  from  a  region  of  single-phase  flow  by  a  thin 
reaction  zone.  We  show  a  wave  front  moving  to  the  rear.  This  compression 
front  would  not  arise  in  the  Langweiler  cycle  but  would  be  expected  for 
other  types  of  burning  schedules.  The  compression  wave  may  also  reflect 
from  the  breech  and,  at  a  later  time,  be  observed  traveling  in  the  opposite 
direction.  In  the  ideal  representation  of  Figure  3.2,  all  the  chemical 
energy  of  the  traveling  charge  is  released  in  a  thin  layer.  We  accordingly 
represent  it  as  a  discontinuity  and  we  show  a  discontinuous  drop  in  pressure 
as  we  pass  from  the  unreacted  to  the  reacted  side  of  the  combustion  zone. 

We  also  show  the  pressure  field  dropping  as  we  move  through  the  unreacted 
traveling  charge  towards  the  base  of  the  projectile.  This  pressure  drop  is 
expected  to  be  linear,  if  the  traveling  charge  is  sufficiently  rigid,  and  is 
analogous  to  the  parabolic  Lagrange  pressure  drop  which  occurs  in  the 
propulsion  gas  in  a  conventional  charge.  We  therefore  note  that  while  the 
pressure  distribution  of  Figure  3.2  is  clearly  different  from  that  of  Figure 
3.1,  both  represent  the  propulsion  of  the  projectile  as  due  to  a  pressure 
which  is  less  than  the  spacewise  maximum. 

When  the  traveling  charge  is  compared  with  the  conventional  charge  in 
terms  of  the  relative  ratios  of  the  base  pressure  to  the  spacewise  maximum 
it  is  not  obvious  that  the  one  concept  is  necessarily  superior  to  the  other. 
Moreover,  elementary  interior  ballistic  theory  is  not  much  help  since  the 
Lagrange  characterization  of  the  pressure  is  not  expected  to  be  accurate 
even  for  conventional  charges  at  the  values  of  C/M  of  interest. 

It  is  clear  that  theoretical  comparison  of  conventional  and  ideal  end¬ 
burning  traveling  charges  can  only  be  conducted  by  reference  to  a  continuum 
model  in  which  the  pressure  gradient  is  developed  as  a  natural  part  of  the 
solution.  The  BRLTC  Code9  was  developed  to  permit  such  theoretical 
comparisons.  The  products  of  combustion  were  modeled  as  an  inviscid,  non¬ 
reacting  one-dimensional  gas  flow  subject  to  the  covolume  equation  of  state. 
The  unreacted  traveling  charge  was  modeled  as  either  rigid  or  as  a  one- 
dimensional  elastic  continuum.  The  reaction  zone  was  represented  as  a 
discontinuity  across  which  the  solid  propellant  was  transformed  to  final 
products  of  reaction.  A  number  of  combustion  laws  were  encoded.  The 
regression  rate  could  be  specified  as  a  function  of  pressure  or  tailored  to 
yield  a  predetermined  value  of  preasure  on  the  unreacted  side  or  of 
projectile  acceleration  or  of  the  Mach  number  of  the  combustion  products 
relative  to  the  regression  front.  Branching  between  the  various  laws  was 
also  admitted.  The  code  was  subsequently  extended  to  incorporate  multiple 
increment  traveling  charges  and  to  provide  a  representation  of  booster 


combust  ion. ^  Hie  booster  was  treated  as  a  homogeneous  mixture  and 
computations  of  the  pressure  gradient  were  not  expected  to  be  reliable  as 
the  mass  of  the  booster  was  Increased  to  become  comparable  to  that  of  the 
traveling  charge.  Accordingly,  the  NOVATC  Code  was  developed^  to  provide  a 
fully  two-phase  treatment  of  the  booster  propellant  and  its  products  of 
combustion.  The  ideal  representation  of  the  traveling  charge  combustion 
zone  as  a  discontinuity  was  nevertheless  retained  in  N0VATC. 

Figure  3.3  illustrates  the  problem  of  interest  here.  Ve  consider  a 
hybrid  charge  consisting  of  a  conventional  booster  increment  and  a  traveling 
charge  increment.  Our  approach  is  applicable  to  all  values  of  the  masses  of 
each  of  the  increments  relative  to  the  projectile  mass.  The  combustion 
model  for  the  traveling  charge  is  extended  relative  to  that  of  NWATC.  We 
still  assume  the  existence  of  a  thin  reaction  layer  at  the  base  of  the 
traveling  charge.  However,  this  reaction  zone  yields  a  mixture  of  final 
products  of  combustion  and  intermediate  products.  The  intermediate  products 
react  at  a  finite  rate  with  the  result  that  the  traveling  charge  flame 
thickness  also  becomes  finite.  By  varying  the  reaction  rate  we  may  vary  the 
thickness  of  the  reaction  zone.  When  short,  the  zone  should  approximate  the 
behavior  of  the  ideal  traveling  charger  of  Figure  3.2.  When  sufficiently 
extended,  the  zone  should  cause  the  release  of  energy  by  the  intermediate 
products  to  yield  a  pressure  distribution  similar  to  that  of  a  conventional 
charge.  Assuming  that  we  have  identified  an  ideal  traveling  charge  which  is 
ball i stical ly  superior  to  a  conventional  equivalent  subject  to  the 
assumption  of  a  thin  reaction  zone,  we  may  then  allow  the  reaction  zone  to 
become  finite  and  determine  how  the  performance  advantage  of  the  traveling 
charge  is  eroded  as  the  reaction  zone  increases  in  length. 

Numerical  simulations  of  the  flow  illustrated  in  Figure  3.3  are 
performed  using  the  XNOVAKTC  (XKTC)  Code.  The  region  between  the  breechface 
and  the  base  of  the  travel  ing  charge  is  modeled  as  a  heterogeneous, 
multiphase  flow  which  is  macroscopic  ally  one-dimensional.  The  flow  in  this 
region  is  considered  to  consist  of  the  solid  booster  propellant  and  a 
mixture  of  combustion  products.  The  combustion  products  include  those  of 
the  booster  propellant  and  both  intermediate  (TC-I)  and  final  (TC-F) 
products  of  combustion  of  the  traveling  charge.  We  distinguish  between  the 
velocities  and  temperatures  of  the  solid  propellant  and  those  of  the 
products  of  combustion.  We  also  have  as  a  field  variable  the  porosity  or 
fraction  of  a  unit  volume  occupied  by  the  mixture  of  combustion  products. 

The  mixture  of  combustion  products  is  multiphase  but  homogeneous,  all 
species  having  the  same  velocity  and,  except  as  specifically  noted 
otherwise,  the  same  temperature.  An  arbitrary  scheme  of  chemical  reactions 
is  permitted  to  occur  in  the  mixture  of  combustion  products.  The  reactions 
may  be  either  of  the  Arrhenius  type  or  pressure  dependent. 
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Figure  3,5  Structure  of  Flow  for  Hybrid  Charge  and  Finite  Reaction  Zone 
For  Traveling  Charge  Products  of  Combustion 


Hie  balance  equations  are  partial  differential  equations  and_have  a 
>artially  hyperbolic  structure.  They  are  integrated  using  a  two-step  finite 
lifference  scheme  of  the  MacCormack  type^  supplemented  by  characteristic 
;orms  at  the  boundaries.  The  chemical  reaction  rate  equations  are 
.ntegrated  using  a  simple  predictor/corrector  scheme  which  is  adequate 
provided  that  the  rates  are  comparable  to  the  hydrodynamic  time  scales. 

The  traveling  charge  may  be  modeled  either  as  rigid  or  as  a  one- 
limensronal  elastic  continuum.  Integration  in  the  latter  case  is  also  by 
:he  MacCormack  scheme.  The  boundary  conditions  at  the  base  of  the  iguited 
traveling  charge  consist  of  finite  balances  of  mass,  momentum  and  energy 
together  with  a  burn  rate  law.  Prior  to  ignition  the  boundary  conditions 
ire  simple  statements  of  contact.  The  projectile  is  taken  to  move  as  a 
rigid  body  opposed  by  friction  dee  to  the  tube  wall  and  the  pressure  in  the 
compressed  air  in  front  of  the  projectile. 

It  may  be  seen  that  with  suitable  data  to  characterize  the  regression 
rate  and  tha  ratio  TC-F:TC-I  the  model  can  be  made  to  simulate  the  first 
step  of  the  process  described  by  White  et  al.*®  A  suitable  reaction  rate 
nodel  then  permits  the  simulation  of  the  second  step,  provided  that 
deconsolidation  is  not  a  dominant  mechanism.  Such  data  are  not  presently 
available  and  it  is  not  an  objective  of  the  present  study  to  attempt  such  a 
direct  simulation.  When  more  precise  simulation  is  required  it  will 
probably  be  appropriate  to  model  the  traveling  charge  increment  according  to 
the  two-phase  analysis  presently  used  for  the  booster.  This  would  allow  a 
more  natural  treatment  of  the  porous  burning  phase  and  the  subsequent 
deconsolidation.  However,  the  increased  fidelity  of  repre  senta  ticn  woul  d 
require  considerably  greater  computer  resources  and  the  solution  would 
involve  a  great  deal  of  numerical  stiffness  which  might  well  require 
algorithm  revisions. 

The  major  limiting  assumption  in  the  present  study  is  that  of  the 
homogeneity  of  the  mixture  of  combustion  products.  If  the  TC-I  species 
consists  of  particles,  they  are  required  to  be  sufficiently  small  that  their 
mechanical  relaxation  times  are  negligible  by  compari  son  wi  th  the 
hydrodynamic  time  scales.  As  is  discussed  by  Wallis.  ^  the  characteristic 
time  for  the  equilibration  of  the  velocity  of  small  spherical  particles,  in 
a  gas  stream  is  given  by 
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“The  Effect  of  Viscosity  in  Hypervelocity  Impact  Cratering" 
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where  d  is  the  particle  diameter,  ps  is  the  particle  density  and  Hg  i •  the 

viscosity  of  the  gas.  We  may  estimate  the  gas  viscosity  at  100CC  as 

7.4  X  10"4  gm/cm-sec,*  The  value  of  ps  will  be  approximately  1.6  gm/ cm* . 

Thus  we  have  x  ~  1.2  X  10*d*  sec  and  since  a  characteristic  hydrodynamic 
time  scale  is  1  msec,  it  follows  that  d  must  be  less  30  microns  for  the 

assumption  of  mechanical  equilibration  to  be  satisfied.  We  note  that  the 
thermal  relaxation  time  is  expected  to  be  of  the  same  order  as  the 
mechanical  relaxation  time. 

When  we  discuss  the  numerical  solutions  in  Section  3.3  we  will  consider 
TC-1  particle  diameters  considerably  larger  than  30p.  Failure  to  consider 
particle  slip  is  not  considered  to  be  a  serious  omission  in  the  context  of 
the  present  study.  The  general  relationship  between  reaction  zone  thickness 
and  ballistic  performance  is  not  expected  to  be  influenced  strongly  by  the 
assumption  of  homogeneity.  The  assumption  of  homogeneity  may  be  of  greater 
concern  if  attempts  are  made  to  simulate  more  directly  the  behavior  reported 
by  White  et  al.^  Apart  from  the  previously  mentioned  possibility  of 
representing  the  traveling  charge  as  a  two-phase  region  by  means  of  XKTC,  we 
note  that  research  is  in  progress  at  BRL*®  and  in  France*^  to  model  the 
combustion  of  the  traveling  charge  on  a  more  fully  non-equil  ibrium  basis. 

3.1  Governing  Equations 

We  confine  our  discussion  to  a  statement  of  the  balance  equations,  the 
equation  of  state  of  the  mixture  of  combustion  products,  the  reaction  rate 
law  used  in  the  present  study,  and  the  boundary  conditions  at  the  base  of 
the  traveling  charge.  Reference  will  be  made  to  previous  reports  for 
further  discussion,  particularly  in  respect  to  the  constitutive  laws.  Our 
main  interest  here  is  to  show  the  difference  between  the  non-equil ibrium 
treatment  of  the  heterogeneous  mixture  consisting  of  the  solid  propellant 
and  its  produ'.ts  of  combustion  and  the  equilibrium  treatment  of  the 
combustion  products  which  are  viewed  as  a  homogeneous  multiphsse  mixture. 
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3.1.1  Balance  Equations  for  the  Mi  it  ore  of  Combustion  Product » 

These  balance  equations  were  developed  in  the  previous  report.^  We 
reproduce  them  in  full  here  even  though  they  include  certain  terms  which  are 
not  used  in  the  present  study.  We  emphasize,  however,  that  the  full 
equations  stated  here  are  completely  linhed  to  the  traveling  charge  boundary 
conditions. 

The  mixture  of  combustion  products  is  assumed  to  be  homogeneous.  All 
species  are  assumed  to  have  the  same  velocity.  In  the  previous  report^  we 
also  assumed  all  species  to  have  the  same  temperature.  As  we  discuss  in 
Section  3.1.3  we  extend  this  in  the  present  work  to  allow  the  species  to  be 
partitioned  into  two  classes,  one  which  is  themally  equilibrated  and  the 
other  which  is  thermally  isolated.  This  extension  was  thought  to  be 
appropriate  for  the  treatment  of  a  mixture  which  contains  burning  particles. 

As  in  Reference  2  we  assume  that  the  flow  is  qua  si-one-dimensional  and 
we  recognize  the  possible  presence  of  reactive  sidewalls  attached  to  the 
tnbe  or  to  the  centerline.  Combustion  of  the  sidewalls  causes  variations  in 
flow  area  and  results  in  mass  addition  to  the  mixture  of  combustion 
products.  It  is  assumed,  however,  that  there  is  no  projectile  afterbody 
intrusion  to  be  considered  simnl  tane ously  with  the  traveling  charge  boundary 
condi tion. 

We  assume  that  we  have  a  total  of  N  chemical  species  which  may  be 
either  gas-  or  condensed-  phases.  A  total  of  E  chemical  reactions  takes 
place  in  the  mixture  of  combustion  products.  A  total  of  J  types  of 
propellant  are  present  in  a  given  cross-section  of  the  flow. 

We  take  A  to  be  the  area  of  the  cross-section.  We  assume  A  is  defined 
by  the  inner  surface  of  the  reactive  layer  on  the  tube  wall  and  the  outer 
surface  of  the  reactive  layer  on  the  centerline  and  that  A  excludes  the 
region  occupied  by  the  unburned  igniter.  The  mass  balance  for  the  mixture 
may  be  written  as 


d  .  .  3  . 

-  epA  +  —  epAu 

dt  dz 


J  N 

•  •  V  .  V  • 

9  +  “s  +  V  +  2  “j  -  L  wi 

*  j=l  i=l 
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where  e  is  the  porosity)  p,  the  density)  u,  the  velocity)  t,  time)  z,  axial 
distance;  <p,  rate  of  addition  of  igniter  products  per  nnit  volume*  m#e  and 
m s rate  of  addition  per  unit  volume  of  outer  and  inner  sidewall  products 
respectively)  mj,  rate  of  addition  per  unit  volume  of  type  j  propellant 
products;  and  w^  is  the  rate  of  deposition  of  species  i  on  the  surface  of 
the  solid  propellant.  We  note  that  when  the  time  dependen< s  of  A  is  due 
solely  to  the  combustion  of  the  igniter  we  have  3A/dt  =  pA,  p^t  where  pjg  is 
the  density  of  the  unburned  igniter. 

The  momentum  equation  may  be  written  as 


Du  dp 

ep  —  +  £go  -L. 

Dt  d  z 


-f 

s 


(3.1.2) 


where  D/Dt  is  the  convective  derivative  along  the  mixture  streamline;  g0,  a 
constant  to  reconcile  units;  f8,  the  interphase  drag)  and  Up,  the  velocity 
of  the  solid  propellant.  The  energy  balance  may  be  written  in  the  following 

f  orm  : 


1 


(3.1.3) 


Here  we  have  e,  the  thermal  part 
the  tube  walli  ,  heat  loss  to 
of  the  igniter;  e#e  and  e#i»  the 


of  the  internal  energy;  qw»  heat  loss  to 
propellant  type  ji  ejg,  the  chemical  energy 
chemical  energies  of  the  outer  and  inner 
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sidewall  products)  Cy^  and  Pc^,  the  constant  volnne  heat  capacity  and 
density  of  condensed-phase  species  i;  Q|.»  the  heat  release  per  unit  mass  of 
reaction  kj  and  tj.  is  the  rate  of  reaction  k.  Ve  note  that  the  presence  of 
the  apparent  heating  term  Q^r^  is  due  to  the  convention  adopted  here  of 
regarding  e  as  excluding  the  chemical  bonding  energy. 

Finally,  we  have  the  governing  equation  for  each  of  the  i  species 
which  constitute  the  mixture  of  combustion  products 


Dt 
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j=l 


N  K 


w .  +  Y.  5  w  + 
1  1  L  n 


n=l 


k=l 


(3  .1.4) 


where  is  the  mass  fraction  of  species  i>  Yjg^,  YSe^,  Ysii  are  the  mass 
fractions  of  species  i  in  the  products  of  combustion  of  the  igniter,  the 
outer  sidewall  and  the  inner  sidewall  respectively)  Y- .  i$  the  mass 
fraction  of  species  i  produced  by  the  near  field  (fizz)  reaction  of 
propellant  jj  and  rik  is  the  rate  of  production  of  species  i  by  reaction  k. 

For  computational  purposes  it  is  convenient  to  eliminate  the  derivative 
of  e  from  Equation  (3.1.3).  Let  £  j'  represent  the  right  hand  side  of 
Equrtion  (3.1.3)  and  let  represent  the  right  hand  side  of  Equation 
(3.1.4)  for  the  i-th  species.  Then  if  c  is  the  the  isentropic  sound  speed 

at  constant  composition  it  follows  that  the  energy  equation  may  be  restated 

as 
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Ve  note  that  certain  of  the  terns  in  Equations  (3.1.1)  -  (3.1.4)  will 
not  be  exercised  in  the  present  study.  Ve  will  have  $  =  ms  =  bs.  =  w,  =  0. 

3.1.2  Balance  Equations  for  the  Solid  Propellant 


The  velocity  and  temperature  of  the  solid  booster  propellant  are 
distinguished  from  those  of  the  mixture  of  combustion  products.  Ve  have  the 
balances  of  mass  and  momentum  for  the  solid  propellant  in  the  following 
f  orm  s : 
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where  pp  is  the  density  of  the  propellant,  up  is  the  velocity,  a  is  the 
intergranular  stress  and  D/Dtp  i5  the  convective  derivative  along  the  solid 
propellant  streamline. 

Since  the  solid  propellant  is  assumed  to  be  incompressible  we  do  not 
state  an  energy  balance.  The  thermal  property  of  interest  is  the  surface 
temperature  which  is  initially  deduced  from  the  interphase  heat  transfer  and 
a  solution  of  the  heat  conduction  equation  applied  to  the  interior  of  the 
solid  propellant.  When  the  surface  temperature  satisfies  an  ignition 
criterion,  the  heat  transfer  condition  is  replaced  by  a  steady-state 
combustion  law. 


3.1.3  Constitutive  Laws 


The  constitutive  laws  required  for  closure  include  the  mixture  equation 
of  state,  the  intergranular  stress  law,  the  interphase  drag  and  heat 
transfer  correlations,  the  wall  heat  loss  correlation,  the  ignition 
criterion,  the  booster  burn  rate  law  and  the  chemical  reaction  rate  lawa. 
Here  we  discuss  only  the  mixture  equation  of  state  and  the  reaction  rate 
laws  as  these  incorporate  some  modifications.  Reference  may  be  made  to 

earlier  work,  for  a  discussion  of  the  other  constitutive  laws.  1*2,6*20 

Ve  have  characterized  the  mixture  in  terms  of  density,  p,  pressure,  p, 

internal  energy,  e  and  species  mass  fractions  i  =  1 . . *N.  Ve  also 

introduce  the  temperature  T  and  we  assume  that  the  mixture  obeys  an 
effective  covolume  equation  of  state 


p(l  -  bp) 


(3.1.9) 


where  bis  the  effective  covolume,  Rg  is  the  universal  gas  constant  and 

is  the  effective  molecular  weight  of  the  mixture.  Moreover,  since  e  is 
understood  to  exclude  the  energy  of  chemical  bonding  we  have  the  caloric 
equation  of  state 


e  =  c  T 
v 


(3.1.10) 


where  cy  is  the  effective  specific  heat  at  constant  volume  for  the  mixture. 

We  assume  for  the  moment  that  thermal  equilibrium  prevails  among  the 
species  so  that  all  have  the  same  temperature  T,  Then  the  value*  of  b  and 
«w  follow  from  a  consideration  of  the  covolume  equation  of  state  for  the 

gas-phase  components  of  the  mixture.  Consider  a  unit  volume  of  the  mixture. 


20  Gough,  P.  S.  "Extensions  to  NOVA  Flamespread  Modeling  Capacity” 
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The  condensed  phases  occupy  a  fraction  of  the  volume  equal  to 
P*i 

-  where  we  assume  species  i  =  ng  +  l, . ,N  to  be  condensed 

"Ci 

phases  and  pc^  is  the  density  of  condensed-phase  species  i.  Within  the  jnit 
volume  are  pY^/M^  moles  of  gas-phase  species  i,  i  =  l,.....,ng  where  is 
the  molecular  weight  of  species  i.  The  gas-phase  molecules  occupy  a 


N 

5 

i=ng+! 


n 

e 

volume  )  pY  .b .  where  b.  is  the  covolume  of  species  i.  The  mixture 
L  11  i 

i»l 

consisting  of  the  gas-phases  alone  evidently  satisfies  the  covolume  equation 
of  state 


[■-  5  £-i.v,U.5£ 

L  i=n„+l  Ci  i=l  J  i=l  Wi 


Comparing  (3.1.11)  with  (3.1.9)  and  (3.1.10)  we  see  that 
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With  all  species  thermally  equilibrated  it  is  clear  that  the  specific  heats 
obey 
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=  5  V 


(3.1.14) 
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where  cp  is  the  specific  heat  at  constant  volume.  We  also  have  y,  the  ratio 
of  specific  heats 


Y  =  c  /c 
P  v 


(3 .1 .16) 


We  also  recall  that  for  the  covolume  equation  of  state 


c  -  c  = 

P  v 


M 


so  that  from  (3.1.9),  (3.1.10)  and  (3.1.16)  we  have 


e  =  - ? -  [1  -  bp) 

<Y  -  1)P 


(3.1.17) 


We  emphasize  the  importance  of  y  in  respect  to  the  thermal  response  of 
the  mixture.  Let  heat  AQ  be  added  to  the  mixture  at  constant  volume  and 
negligible  initial  pressure.  Then,  neglecting  the  covolume,  we  have 
p  =  (y  -  l)pAQ.  The  increase  in  pressure  is  proportional  to  y  -  1.  For 
gas-phase  species  we  will  have  typically  Cp^/cv^  ~  1.25.  However,  for 

condensed-phase  species  we  will  have  «p^/cv^  ~  !•  If  the  mixture  consists 


of  equal  mass  fractions  of  a  gas  and  a  solid  then  we  will  have  y  ~  1.17.5  and 
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the  pressure  increase  due  to  AQ  will  be  seen  to  be  one  half  the  increase 
that  would  occur  if  the  mixture  consisted  only  of  a  gas.  This  is  the 
intuitively  expected  result,  but  it  is  important  to  see  that  it  is  conveyed 
through  the  dependence  of  7  on  the  composition  of  the  mixture. 

We  have  assumed  thu?  far  that  all  species  have  the  same  temperature. 

N^w  suppose  that  certain  condensed  species  are  thermally  isolated  as  would 
be  the  case  for  large  particles  or  for  small  particles  surrounded  by  a  flame 
zone.  Let  T  be  the  temperature  of  the  thermally  equilibrated  species.  Then 
it  is  easy  to  see  that  (3.1.12)  and  (3.1.13)  still  apply  since  T  is  the 
temperature  of  the  gases.  Moreover,  (3,1.16)  and  3.1.17)  also  apply  but  it 
is  necessary  to  replace  (3.1.14)  and  (3.1.15)  with 
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where  E ;  =  1  if  species  i  is  thermally  equilibrated  and  Ej  =  0  if  species  i 
is  insulated.  We  emphasize  that  this  simple  modification  is  only 
appropriate  when  the  insulated  species  are  condensed-phases.  A  two- 
temperature  gas  mixture  is  not  considered  here. 

The  second  constitutive  lsw  of  interest  here  is  that  for  the  rate  of 
chemical  reaction  which  appears  in  Equation  (3.1.4).  We  have  previously 
assumed  r^  to  be  given  by  an  Arrhenius  law.**  Here  we  wish  to  model  the 
case  in  which  reaction  k  represents  the  combustion  of  condensed  species  i  by 
normal  surface  regression.  Thus  r^  i*  the  negative  value  of  the  rate  of 
decomposition  of  species  i  per  unit  volume. 

Let  species  i  consist  of  an  aggregate  of  droplets  or  particles  which 
are  locally  identical.  Let  Vp^  and  SPi  be  the  volume  and  surface  area  of 
one  particle.  Then  the  number  of  particles  per  unit  volume  is  given  by 
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If  i$  the  rate  of  surface  regression  it  follows  that 


r  =  -  p  S  n  d. 

ik  Cj.  pi  Pi  i 


which,  in  view  of  (3.1.20)  implies 
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where  we  have  introduced  Dpi  =  6VPi/Spi  as  the  effective  diameter  of  species 
i.  We  may  assume  that  d^  is  given  in  the  usual  form,  having  an  exponential 
dependence  on  pressure. 

In  the  present  study  we  characterize  the  combustion  of  an  aggregate  of 
particles  in  tems  of  a  characteristic  diameter  Dp^  and  the  burn  rate 
d  =  Bpn.  However,  we  do  not  follow  the  changes  in  particle  diameter  aa 
combustion  proceeds.  This  is  not  difficult  to  do  but  it  wss  not  thought  to 
be  worth  the  additional  computational  harden  in  the  present  context. 

However,  we  note  that  when  Dpi  is  allowed  to  vary,  we  expect  that  Y*  and  Dpi 
will  tend  to  zero  together,  maintaining  a  finite  value  of  rik  and  aasnring  a 
clean  burn  out  of  species  i.  If  the  value  of  Dp  is  kept  constant,  as  is 
done  here,  rik  tends  asymptotically  to  zero  with  Y^  and  the  particles  never 
quite  burn  out  completely. 

Accordingly,  we  use  the  following  law  to  describe  the  pressure 
dependent  rate  of  reaction  of  species  i 
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where  is  the  burn  rate  pre-exponent  and  nj  is  the  exponent.  We  note,  in 
comparing  (3,1.22)  with  (3.1.21)  that  we  have  simply  dropped  the  factor  of 
T^.  When  we  characterize  the  rate  of  reaction  by  the  inpnt  datum  Dp^  we  may 
interpret  Dp^  =  Dp^Y^  as  the  initial  particle  diameter. 

3.1.4  Traveling  Charge  Balance  Equations 

In  order  to  distinguish  them  from  the  corresponding  quantities  for  the 
booster  charge,  we  denote  the  traveling  charge  solid-phase  state  variables 
by  the  subscript  tc.  Thus  we  have  ptc,  u^c  and  crtc  which  denote  the 
density,  velocity  and  pressure  in  the  traveling  charge.  We  have  the 
bal*nces  of  mass  and  momentum 
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where  fw  is  the  frictional  force  exerted  on  the  traveling  charge  by  the  wall 
of  the  tube.  It  is  implicitly  assumed  that  the  traveling  charge  moves 
through  a  constant  area  section  of  the  tabe. 

3,1.5  Boundary  Conditions  at  the  Base  of  the  Traveling  Charge 

When  the  traveling  charge  has  not  ignited,  the  physical  boundary 
condition  at  its  base  expresses  the  contact  of  the  mixture  of  combustion 
products  and  the  non-penetration  of  the  solid  propellant.  Ignition  of  the 
traveling  charge  is  taken  to  occur  after  a  predetermined  delay.  Subsequent¬ 
ly,  the  physical  boundary  conditions  consist  of  the  finite  balances  of  mass, 
momentum  and  energy  at  the  regression  front  and  either  one  or  two  data  to 
determine  the  regression  rate.  (Hie  datum  is  required  to  determine  the 
regression  rate  if  the  products  are  subsonic  relative  to  the  front  and  two 
data  are  required  if  the  products  are  sonic  relative  to  the  front.  The 
condition  of  supersonic  products  —  the  strong  deflagration  wave  —  is  not 
admitted. 
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The  finite  balances  of  mass,  momentum  and  energy  may  be  stated  as 
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Here  we  have  us  as  the  velocity  of  the  regression  front  relative  to  the  gun 
tube  so  that 
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where  rtc  is  the  regression  rate  relative  to  the  unburned  traveling  charge. 

We  have  used  the  .subscript  tc  to  denote  properties  of  the  traveling  charge 
and  to  distinguish  them  from  the  booster  solid  propellant.  We  note  that  we 
have  introduced  Atc,  the  crOss-se ctio&al  area  of  the  traveling  charge 

increment,  which  is  not  necessarily  assumed  equal  to  that  of  the  tnbe,  in 
order  to  account  for  the  possible  presence  of  an  external  liner  used  to 
support  the  traveling  charge. 

It  remains  to  discuss  the  conditions  which  specify  rtc<  We  assume  that 
rtc  obeys  any  of  the  following  laws: 

(a)  Measured  burn  rate  r^c  is  given  as  s  function  of  the  pressure  on 
either  side  of  the  flame  in  either  exponential  or  tabular  form. 

(b)  Ideal  burn  rate — r^,,  is  chosen  so  as  to  yield  a  predetermined 

value  of  pressure  on  either  side  of  the  fine,  or  to  yield  a  predeterminsd 
acceleration  of  the  projectile,  or  to  yield  a  predetermined  value  of  the 
Mach  number  of  the  combustion  products  relative  to  the  flame. 

(c)  Composite  burn  rate — rtc  may  be  required  to  satisfy  a  measured  or 

ideal  burn  rate  law  subject  to  the  constraint  that  the  Mach  number  of  the 
products  be  equal  to  one. 
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It  is  always  assumed  that  the  Mach  number  of  the  products  of  combustion 
of  the  traveling  charge  relative  to  the  flame  is  less  than  or  equal  to  one. 
If  the  Mach  cumber  is  less  than  one  it  is  assumed  that  the  burning  process 
is  acoustically  coupled  to  the  state  of  mixture  of  combustion  products  and 
use  is  made  of  the  appropriate  characteristic  constraint.  If  the  Mach 
number  is  one.  acoustic  coupling  is  not  assumed  and  the  characteristic 
constraint  is  replaced  by  the  condition  of  choking.  In  the  latter  case  the 
traveling  charge  burns  as  a  nozzleless  rocket. 

Ignition  of  the  first  increment  is  assumed  to  occur  at  a  prespecified 
point  in  time.  Following  ignition,  the  full  born  rate  as  described  by  the 
appropriate  law  is  assumed  to  be  reduced  by  a  coefficient  which  increases 
from  zero  to  one  over  a  second  nser-selectabl e  interval.  When  the  fine 
passes  from  one  increment  to  the  next,  a  delay  can  be  specified  for  the 
ignition  of  the  new  increment.  The  burn  rate  achieves  the  value  given  by 
the  appropriate  law  for  the  new  increment  after  an  additional  delay  during 
which  the  actual  value  varies  linearly  in  time  from  the  final  value  for  the 
previous  increment  to  the  full  value  for  the  new  increment. 

A  strong  rarefaction  can  be  formed  when  the  traveling  charge  burns  out. 

As  in  BHLTC,^  we  use  a  simple  wave  solution  for  five  steps  after  burnout  to 
allow  the  state  variables  in  the  mixture  of  combustion  products  to  become 
reasonably  smooth. 

3.2  Numerical  Results 


The  numerical  results  presented  here  satisfy  two  different  objectives. 
First,  vs  seek  to  determine  the  effect  of  a  finite  length  reaction  zone  on 
the  interior  ballistics  of  a  traveling  charge.  A  second  objective,  which 
had  to  be  satisfied  prior  to  the  first,  was  to  complete  the  de"elopment  of 
XJQ'C,  including  in  particular,  the  traveling  charge  option  and  its  linkage 
to  the  chemistry  option.  To  satisfy  these  objectives  we  started  rith  a 
BHLTC  data  base  which  was  considered  to  exhibit  a  reasonable  level  of 
ballistic  benefit  from  a  traveling  charge  increment.  This  data  base  was 
then  made  compatible  with  NOVATC  and  necessary  algorithm  revisions  were 
identified  and  incorporated  to  achieve  a  satisfactory  numerical  solution. 
This  solution  then  served  as  a  benchmark  against  which  the  operability  of 
1KTC  could  be  verifed,  at  least  for  the  ideal  case  of  an  infinitesimal 
traveling  charge  combustion  zone.  Finally,  the  chemistry  option  was  invoked 
to  generate  a  finite  fine  zone  in  the  XKTC  solutions  and  the  effect  of  the 
flame  thickness  on  ballistic  performance  was  appraised. 
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Table  3.1  XKTC  Input  Data  for.  Ncr^laal  Simulation  of  Txaveiing  Charge  w ith 

Finite  Flame  Thicknesi 


CONTROL  BATA 

LOGICAL  VARIABLES: 

PRINT  T  DISK  WRITE  F  DISK  READ  F 

I.B.  TABLE  T  FLAME  TABLE  F  PRESSURE  TABLE(S)  F 

EROSIVE  EFFECT  0  WALL  TEMPERATURE  CALCULATION  0 

BED  PRECOMPRESSED  0 

HEAT  LOSS  CALCULATION  1 

BORE  RESISTANCE  FUNCTION  2 

TRAVELING  CHARGE  OPTION  <0=NO,  1=YES)  1 

CONSERVATIVE  SCHEME  TO  INTEGRATE  SOL  ID- PHASE 

CONTINUITY  HHJAIION  (0=NO,  0LD>  1=YES,NEW)  0 

KINETICS  MODE  (CNNONE,  1=GAS-PB>SE  CNLY,  2=B0TH  PHASES)  1 

TANK  GUN  OPTION  (0=NO,  1=YES)  0 

INPUT  ECHO  OPTION  0 


INTEGRATION  PARAMETERS 


NUVBEk  OF  STATIONS  AT  WHICH  DATA  ARE  STORED 

NU1CER  OF  STEPS  BEFORE  LOGOUT 

TIME  STEP  FOR  DISK  START 

NUMH21  OF  STEPS  FOR  TERMINATION 

TIME  INTERVAL  BEFORE  LOGOUT*  SEC) 

TINE  FOR  TERMINATION  (SEC) 

PROJECTILE  TRAVEL  FOR  TERMINATION  (IN) 

MAXIMUM  TIME  STEP  (SEC) 

STABILITY  SAFETY  FACTOR 

SOURCE  STABILITY  FACTOR 

SPATIAL  RESOLUTION  FACTOR 

TIME  INTERVAL  FOR  I.B.  TABLE  STORAGE(SEC) 

TIME  INTERVAL  FOR  PRESSURE  TABLE  STORAGE  (SEC) 


30 

5000 

0 

5000 

.100  X  10”* 
10.00 

157.48  (400  cm) 

.100  X  10-1 
2.00 
.200 
.010 

.100  X  10~* 
.100  X  10”* 
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FILE  OOUNTERS 


NUMBER  OF  STATIONS  TO  SPECIFY  TUBE  RADIOS  2 

NUMBER  OF  TIMES  TO  SPECIFY  PRIMER  DISCHARGE  0 

NUWER  OF  POSITIONS  TO  S  Hi  GEY  PRIMER  DISCHARGE  0 

NDJCER  OF  ENTRIES  IN  BORE  RESISTANCE  TABLE  0 

NUIBER  OF  ENTRIES  IN  WALL  TEMPERATURE  TABLE  0 

NUMBER  OF  ENTRIES  IN  FORWARD  FILLER  ELEMENT  TABLE  0 

NOIBER  OF  TYPES  OF  PROPELLANTS  1 

NUWER  OF  BORN  RATE  DATA  SETS  1 

NUMBER  OF  ENTRIES  IN  VOID  FRACTION  TABLE(S)  0  0  0 

NUMBER  OF  ENTRIES  IN  PRESSURE  HISTORY  TABLES  0 

NOIBER  OF  ENTRIES  IN  REAR  FILLER  ELEMENT  TABLE  0 


GENERAL  PROPERTIES  OF  INITIAL  AW  LENT  GAS 


INITIAL  TEMPERATURE  (R) 
INITIAL  PRESSURE  (PSI) 
MOLECULAR  WEIGHT  ( LB  M/ LB  MOL) 
RATIO  OF  SPEGFIC  HEATS 


540.0  (300K) 

14.7  (0.101  MPa) 

28.960 
1.400 


GBH&AL  PROPERTIES  OF  PROPELLANT  BED 


INITIAL  TEMPERATURE  (R) 


540.0  (300K) 
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ROPERIES  OF  R0PELLANT  1 


PROPELLANT  TYPE 

BOOSTER 

MASS  OF  PROPELLANT  (IBM) 

.5327 

(242  kb) 

DENSITY  OF  PROPELLANT  ( LBM/ IN* ) 

.06  00 

(1.66  gm/  am* ) 

FORM  FUNCTION  INDICATOR 

7 

OUTSIDE  DIAMETER  (IN) 

.2415 

(0.613  cb) 

INSIDE  DIAMETER  (IN) 

.0281 

<0.071  «) 

LENGTH  (IN) 

.57  95 

(1.472  cb) 

NUMBER  OF  PERFORATIONS 

7. 

S.OT  WIDTH  (NF0RM~11) 

OR  SCROLL  DIA.  <NF0RM=13)  (IN) 

0.0000 

PROPELLANT  STACKED  (0=NO,  1=YES)  0 

ATTACHMENT  CONDITION  ( 0=FREE,  1= ATTACH  ED  TO  TUBE, 

2=ATTACBED  TO  PROJECTILE)  0 

BOND  STRENGTH  (LBF) 

(N.B.  ZERO  DEFAULTS  TO  INFINITY)  0. 

UROLOGICAL  ROPERIES 


SPEED  OF  COMPRESSION  WAVE 
IN  SETTLED  BED  (IN/ SEC) 

SETTLING  POROSITY 

SPEED  OF  EXPANSION  WAVE  (IN/ SEC) 

POISSCN  RATIO  <-) 


17400  (44196  cn/iec) 

1.0 

30000.  (127000  on/ sc  c) 

0-0 


SOLID  PHASE 

MAXIMUM  PRESSURE  FOR  BUFN  RATE 
DATA  (LBF/IN* ) 

BURNING  RATE  PRE-EXPCN0JTIAL 
FACTOR  (IN/  SEC-PSIBN) 

BURN  I hT.  RATE  EXPONENT 

BURNING  RATE  CONSTANT  (IN/  SEC) 

IGNITION  TEMPERATURE  (R) 

THERMAL  CONDUCTIVITY  (LBF/SEC-R) 

fmwnut  r\Trmio  tirTTV  /  /  CT7P  \ 

iilLAC  UUj  U±  IT'UJ  A  W  LA.  A  \  Jil  /  wwv/ 


THBUOCHBMISIST 

100000.  (689.5  MPa) 

.1790  X  10~*  (0.3368  cm/*e  c-MPaBN) 

.8630 

0.0000 

539.0  (299. 4K) 

.2770  X  10_1  (2.22  X  10_*  J/cm--aec-I 

,1345  X  !Q~*  (8.677  X  10“4  cw*/*®o) 

.600 


EMISSIVITY  FACTOR 


GAS  FHALE  THERMOCHEMISTRY 


CHEMICAL  ENERGY  RELEASED  IN  BURNING(LBF-IN/LBM)  .17280  X  108  (4304  J/gm) 

MCLECULAR  WEIGHT  (IBM/ IBMOL)  19.4000 

RATIO  OF  SPECIFIC  HEATS  1.2500 

COT OLD ME  32.9000  (1.189  cn’/gm) 


LOCATION  OF  PACKAGE* S) 


PACKAGE 

LEFT  BDDY 

RIGHT  BDDY 

MASS 

INNER  RADIOS 

00 TER  RADIUS 

(IN) 

(IN) 

(IBM) 

(IN) 

(IN) 

1 

0.000 

7.963 

0.533 

0.000 

0.000 

(20.22  cm) 

(241.8 

gm) 

PARAMETERS  TO 

SPECIFY  TUBE  GEOMETRY 

DISTANCE(IN) 

RADIUS(IN) 

0.000 

.787 

200.000  (508 

cm) 

.787  (2. 

0  cm) 

THERMAL  PROPERTIES  OF  TUBE 

THERMAL 

CONDUCTIVITY 

(LBF/SEC-R) 

7.770  (0.662 

J/cm-  se  c-K) 

THERMAL 

DIFFUSIVm 

(IN* /SEC) 

.2280  X  10~* 

(0.147  an* / ne 

EMISSIVITY  FACTOR 

.700 

INITIAL 

TEMPERATURE 

(R) 

540.00  (300K) 

PROJECT  ILK  AND  RIFLING  DATA 


INITIAL  P06ITI(N  OF  BASE  OF  PROJECTILE*  IN)  0.000 

MASS  OF  PROJECTILE  (IBM)  0.000 

POLAR  MOMENT  OF  INERTIA  (IBM- IN1)  0.000 

ANGLE  OF  RIFLING  (DEG)  0.000 
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CHEMISTRY  OPTION  DATA 


NUMBER  OF  SPECIES  4 

NUMBER  OF  GAS-PHASE  REACTIONS  1 

NUMBER  OF  SOLID-PHASE  REACTIONS  0 

PROPERTIES  OF  SPECIES 


NAME 

PHASE 

CV 

IBF-1N7  LBM-R 

CP 

LBF-INy  LBM-R 

COVCLUMI3 

inVibm 

MOL.VGT 

LB/  LB  MOL 

DENSITY 

IBM/ IN* 

TRANSFER 

CQEF. 

KTEQL 

AIR 

G 

3  824.7 

4780.9 

32.900 

19.400 

0.00000 

0. 

0 

BOOSTER 

G 

3824.7 

47  80.9 

32.900 

19.400 

0.00000 

0. 

0 

TCI 

S 

3  824  .7 

3  824  .7 

0.000 

0.000 

<1 

.06000 
.66  gm/cm* 

0. 

) 

1 

TCF 

G 

3  824.7 

47  80 . 9 

32.900 

19.400 

0  .00000 

0. 

0 

(1.715  (2.144  (1.189 

J/gm-K)  J/gn-K)  cn'/gm) 

COMPOSITION  OF  LOCAL  OOJEUSTICN  PRODUCTS  OF  PROPELLANT  1 


ENERGY 

LB  F- IN/ IBM 

MASS  FRACTIONS  (-) 
Y0(  1)  Y0(  2) 

Y0(  3) 

Y0(  4) 

17280000. 

(4304  J/gra) 

0.00000  1.00000 

0.00000 

0.00000 

COMPOSITION 

OF  LOCAL  OOIBUSTIQN  PRODUCTS  OF  TRAVELING  CHARGE 

ENERGY 

LBF-IN/  IBM 

MASS  FRACTIONS  (-) 
YTC0(  1)  YTC0(  2) 

YTCO(  3) 

YTCO(  4) 

8640000. 

(2152  J/gm) 

0.00000  0.00000 

0.50000 

0.3G000 

COMPOSITION  OF  OOIBUSTIQN 

PRODUCTS  OF 

IGNITES 

ENERGY 

LBF-IN/LBM 

MASS  FRACTIONS  (-) 
Y0(  1)  Y0(  2) 

Y0(  3) 

Y0(  4) 

0. 

0.00000  0.00000 

0 .00000 

0.00000 
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composition  op  ajbiskt  hah 


ENERGY  MASS  FRACTIONS  (-) 

LB F- IN/'  LB M  YO(  1)  '  IO(  2)  TOC  3;  YOC  4) 


864727.  1.00000  0.00000  0 .0-0000  0.00000 

(  215  J/gm) 


GAS- PHASE  REACT  IWi  MXA 


REACTION 

1 

REACTANT  SPECIES  3 

0  0 

0  PRODUCT 

s  real's  4  o 

0  0 

STOICHIOMETRIC  COEFFICIENTS  (  LBM)  1. 

0.  0. 

0. 

1.  0. 

0.  0 

BEAT  OF  REACTION  (IBF-IN/LBM) 

17  280000 . 

(430 4  J/gm) 

PARTI aE  DIAMETER  (IN) 

0.1 

(0,264  c*) 

BURN  RATE  ADDITIVE  CONSTANT  (IN/ SEC) 

0. 

BURN  RATE  COEFFICIENT  ( IN/ SEC- PS  I8  N) 

.17  900  X  10"* 

BURN  RATE  EXPONENT  (-) 

.8650 

T.C.  CONTROL  DATA 


I  MEAL  BURN  RATE  LAI  2 

CONTINUUM  MODEL  OF  UNREACIED  PROPELLANT  1 

NUJBER  OF  PROPELLANTS  1 

PROPELLANT  WALL  FRICTHN  PARAMETER  0 

NUMBER  OF  ENTRIES  IN  PROJECTILE  BORE  RESISTANCE  TABLE  2 
INDICATOR  FOR  AIR  RESISTANCE  1 

NUMBER  OF  ENTRIES  IN  CBlU  EAT  OR  FRICTION  TABLE  0 


INTEGRATION  PARAMETERS 

MAXIMUM  NUMBER  OF  MESH  POINTS 
HINIMOM  MESH  SIZE  (IN) 


11 

.200  (0.508  ck) 


48 


rear,  and  tray,  charge  fiofbties 


T.  C.  DIAMETER  (IN)  1.575 

INITIAL  POSITION  OF  REAR  FACE  OF  PROPELLANT  (IN)  7.963 

PROJECTILE  MASS  (IBM)  .35270 

CHARGE  MASS  (IBM)  .542 

MAXIMUM  PRESSURE  IN  UNREACTED  PROPELLANT  (PSI), 

IF  IDEAL=2  100000. 

MAXI  HIM  MACH  NUMBS  OF  REACT  I  (N  PRODUCTS  .999 

MAXIMUM  ACCELERATION  OF  PROJECTILE  (G)  0. 

RATIO  OF  SPECIFIC  HEATS  OF  AIR  (-)  1.4000 

PRESSURE  OF  AIR  IN  BARREL  (PSI)  14.700 

TEMPERATURE  OF  AIR  IN  BARREL  (R)  540.0 

MOLECULAR  WEIGHT  OF  AIR  IN  BARREL  ( IBM/  LB  MOL)  28.9600 


RESISTIVE  PRESSURE  DUE  TO  OBTURATOR 


TRAVEL 

RESISTIVE  PRESSURE 

(IN) 

(PSI) 

0.000 

800. 

(5.52  MPa) 

.500  (0.27  cm) 

500. 

(3.45  MPa) 

(4.0  ca) 
(20.22  ca) 
(160  ga) 
(245.8  ga) 

(689.5  MPa) 


(0.101  MPa) 
(300K) 
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PROPERTIES  OF  PROHBLLANT  NUWEE  1 


RATIO  OF  SPECIFIC  HEATS  (-)  1.250 

COVOLUBE  (IN* /IBM)  32.900 

MOLECULAR  WEIGHT  ( LB  M/ LB  MOL)  19.400 

CHEMICAL  ENERGY  OF  PROPELLANT  (  LBF-IN/  IBM)  17280000. 

DENSITY  OF  PROPELLANT  (UBM/IN1)  .046 6 

INITIAL  MASS  (IBM)  .5415 

IGNITION  DELAY  (MSEC)  1.800 

DELAY  FOR  IRANS.  TO  FULL  BUSH  RATE  (MSEC)  .100 

BURNING  RATE  ADDITIVE  CONSTANT  (IN/ SEC)  -I 

BURNING  RATE  PRE-EXPONENTIAL  FACTOR  <IN/SEC-PS1®N)  -I 
BURNING  RATE  EXPONENT  (-)  -I 

TC  GRAIN  LENGTH  (IN)  5.966 

LENGTH  BREECH  TO  PROJECTILE  BASE  (IN)  13.929 

COMPRESSION  WAVE  SPEED  IN  PROPELLANT  (IN/ SEC)  118110. 
EXPANSION  WAVE  SPEED  IN  PROPELLANT  (IN/SEC)  0. 

BURN  RATE  FORMAT  (0=EXP>  l=TABULAk)  0 

BURN  RATE  DEPENDENCE  <0=FRES,  1-STRESS )  1 


(1.189  an‘/gm) 

(4304  J/gn) 
(1.290  gm/cm*) 
(160  gm) 


(15.15  cm) 

(35.38  cm) 

(300000.  cm/ sec) 
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A  repre  sentat  ive  XKTC  data  base,  including  the  chemistry  data,  is 
presented  in  Table  3.1.  The  problem  of  interest  involves  a  total  propellant 
mass  of  488  gm  and  a  projectile  mass  of  160  gm  so  that  C/M  is  approximately 
equal  to  3  The  gun  bore  diameter  is  4.0  cm  and  the  projectile  travel  is 
400  cm.  The  charge  is  divided  into  a  booster  component,  consisting  of 
seven-perforation  granular  propellant  and  having  a  mass  of  242  gm,  and  a 
traveling  charge  component  whose  mass  is  246  gm.  The  booster  granulation  is 
selected  to  achieve  a  maximum  breech  pressure  of  approximately  690  MPa  which 
occurs  at  about  1.2  msec  after  the  booster  is  ignited.  Ignition  of  the 
traveling  charge  occurs  at  1.8  msec  when  the  breech  pressure  has  fallen  to 
approximately  250  MPa.  The  rate  of  surface  regression  of  the  traveling 
charge  is  required  to  yield  a  stress,  on  the  unreacted  aide  of  the  flame, 
equal  to  690  MPa,  provided  that  the  Mach  number  of  the  products  does  not 
exceed  0.999.  If  the  stress  cannot  be  achieved  w ithout  violating  the  Mach 
number  constraint,  the  regression  rate  is  chosen  to  yield  a  Mach  number  of 
0.999.  The  solution  within  the  raveling  charge  is  assumed  always  to  be 
acoustically  coupled  to  that  in  the  mixture  of  combustion  products. 

The  d«  ta  of  Table  3.1  were  developed  from  a  BKLTC  data  base  which  we 
refer  to  as  40MTC3.  Certain  modifications  were  incorporated  in  order  to 
arrive  at  the  data  of  Table  3.1,  apart  from  considerations  of  a  chemical 
reaction  in  the  mixture  of  combustion  products.  BHLTC  models  the  booster 
increment  as  n  single-phase  substance  and  the  pressure  gradient  responds 
only  to  the  momentum  of  the  combustion  products.  Accordingly,  BHLTC  tends 
to  underestimate  the  pressure  gradient  when  compared  with  the  more  complete 
two-phase  analysis  of  NOVATC  or  XKTC.  The  NOVA  data  bases  therefore 
incorporated  a  somewhat  lower  burn  rate  coefficient  to  produce  the  same 
maximum  chamber  pressure  as  BHLTC,  It  also  became  necessary  in  the  NOVA 
runr  to  advance  slightly  the  time  of  ignition  of  the  traveling  charge 
increment  in  order  to  assure  burnout  prior  to  muzzle  exit.  In  Table  3.2  we 
compare  the  predictions  of  BHLTC,  NOVATC  and  XKTC  for  the  40MTC3  data  base. 


Table  3  .2 

Code  Dependence 

of  Nominal  Thin  Plane 

TC  Data  Base 

(40MTC3) 

Code 

Max  Press. 

Muzzle  Vel c-oi ty 

Anflb 

Ae% 

(MPa) 

(m/  sec) 

BKLTC* 

699 

287  9 

0.67 

0.36 

NOVATC** 

691 

2860 

-1.9 

0.04 

XKTC** 

6  87 

2854 

-2.1 

-0.95 

Booster  burn  rate  =  0 .00185 p°  •*  * 1 ,  TC  ignition  st  2.0  msec,  maximum  of 
31  points  for  booster  and  TC  combined. 

Booster  burn  rote  =  0 .00017  9p°  ■ *** ,  TC  ignition  at  1,8  msec,  16  points 
for  booster,  maximum  of  11  for  TC, 
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Good  agreement  it  seen  between  the  three  sets  of  cslcnl  stions.  Hie  so 
results  serve  to  demonstrate  the  operability  of  XSTC  for  the  thin  flame 
traveling  charge  calculations  and  to  confirm  the  data  of  Table  3.1  as 
reproducing  the  performance  of  40MTC3  in  BRLTC.  We  note  that  we  have  also 
tabulated  values  of  the  final  percent  mass  and  energy  defects.  As ft  and  Ae%, 
as  indicators  of  numerical  accuracy,  \alues  of  these  quantities  less  than 
1%  would  be  desirable,  but  we  have  accepted  results  in  this  study  for  which 
the  final  defects  were  as  large  as  2  or  3%. 

To  verify  that  these  mass  and  energy  defects  did  not  imply  excessive 
mesh  dependence  we  present  solutions  for  three  different  problems  obtained 
with  30  and  60  mesh  points  for  the  region  occupied  by  the  mixture  of 
combustion  products.  The  first  problem,  identified  as  the  thin  flame  TC, 


Table  3,3  Mesh  Dependence  of  XCTC  Solutions 


Probl em 

Number  of 
Booster 
Mesh  Points 

Max  Press. 
(MPa) 

Muzzle  Vel. 
(m/se  c) 

Am% 

Ae% 

Thin  Flame  TC 

30 

692 

2875 

-2.4 

-1.4 

60 

6  98 

2875 

-1.9 

-1.3 

Conventional  Equival 

ent  30 

689 

2476 

-1.16 

-0.54 

60 

6  9  6 

2483 

-0.3  8 

-0.11 

Fini  te  FI  ame  TC* 

30 

6  84 

2603 

-1  .58 

-2.22 

60 

6  89 

2589 

-1.23 

-1.82 

*  Burn  rate  =  0 .000179p° • * ‘ ,  TC-l/TC-F  =  50/50,  Dp  *  0.254  cm 


is  the  40MTC3  data  base  adapted  to  XKTC.  The  second  problem,  identified  as 
the  conventional  equivalent,  is  a  charge  consisting  entirely  of  granular 
propellant  and  having  the  same  total  mass  as  40MTC3.  The  third  problem, 
identified  as  the  finite  flame  TC,  corresponds  precisely  to  the  data  of 
Table  3.1.  Good  mesh  indifference  is  exhibited  for  all  three  problems  and 
it  is  assumed,  unless  otherwise  explicitly  noted,  that  the  aolntions 
presented  here  have  an  accuracy  of  1  or  2%. 

The  conventional  equivalent  data  base  was  formed  by  taking  the 
traveling  charge  increment  to  be  a  conventional  granular  increment  of  the 
same  mass  and  chemical  energy.  In  order  to  obtain  a  conventional  loading 
density  of  approximately  60%  it  was  necessary  to  move  the  projectile  forward 
a  distance  of  5  cm.  The  travel  was  still  assumed  to  be  400  cm  so  that  the 
conventional  equivalent  result  corresponds  to  a  slightly  longer  gun  than  the 
traveling  charge  results.  This  difference  is  not  believed  to  be 


significant.  The  granulation  of  the  conventional  equivalent  was  varied  to 
achieve  a  maximum  breech  pressure  of  approximately  690  MPa.  Table  3.3 
therefore  shows  a  ballistic  benefit  of  the  traveling  charge  in  the  sense 
that  the  muzzle  velocity  exceeds  that  of  the  conventional  equivalent  by 
400  m/sec  at  the  same  maximum  pressure. 

It  should  be  noted  that  the  maximim  pressure  for  the  traveling  charge 
is  achieved  not  only  at  the  breech  but  also  further  downbore  in  the  regions 
occupied  by  the  traveling  charge  following  ignition.  Ve  do  not  consider  the 
oxtent  to  which  the  ballistic  benefit  is  offset  by  the  potentially  increased 
structural  burden  on  the  tube. 

The  finite  flame  TC  calculation  shown  in  Table  3.3  assumes  that  the 
regression  of  the  base  of  the  traveling  charge  yields  a  mixture  of  final 
combustion  products  (TC-F)  and  intermediate  products  (TC-I)  in  the  ratio 
TC-F/TC-I  =  50/50.  The  intermediates  are  assumed  to  react  according  to 
Equation  (3.1.22)  with  Dp^  =  0.254  cm.  We  understand  this  value  to 

correspond  to  an  initial  particle  diameter  equal  to  0.127  cm. 

We  note  that  the  maximum  pressure  in  the  finite  flame  TC  calculation  is 
essentially  identical  to  that  in  the  two  preceding  cases.  This  will  be  true 
of  all  the  calculations  considered  here.  ‘ccordingly,  the  evaluation  of 
ballistic  performance  is  reduced  to  an  e.'xsination  of  the  muzzle  velocities. 
It  is  evident,  in  Table  3.3,  that  this  particular  finite  fine  TC  data  base 
exhibits  considerable  degradation  of  performance,  the  muzzle  velocity 
exceeding  that  of  the  conventional  equivalent  by  only  125  m/sec. 

In  order  to  relate  the  degradation  of  the  traveling  charge  performance 
benefit  to  the  thickness  of  the  flame  zone  we  .introduce  the  characteristic 
reaction  zone  length  L  defined  as 


„  /L  e PAxY*dx 

L  =  — - 1 —  .  (3.2.1) 

/0L  epAY.dx 


Here  x  is  a  coordinate  originating  at  the  base  of  the  traveling  charge  and  L 
is  the  distance  to  the  breechface.  We  have  e,  porosity)  p,  mixture  density) 
A,  cross-sectional  area  of  the  tube;  and  Yj  is  the  mass  fraction  of  the 
intermediate  product  species  TC-I. 

It  should  be  noted  that  we  have  Y|  in  the  numerator  and  in  the 

denominator.  Hence  L*  is  a  measure  not  only  of  the  length  of  the  reaction 
zone  but  also  the  concentration  of  TC-I. 
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It  i$  easy  to  verify  that  if 


Y 


i  ~ 


0  <  x  <  Lf 
x  >  Lf 


where  YQ  is  constant  and  L£  is  the  flane  thickness,  then 


(3.2,2) 


* 

L 


Y  L. 
o  f 


(3.2.3) 


In  such  a  case  Lj  =  2L*/YC  and  will  be  equal  to  4L*  if  YQ  =  0.5.  The 
distribution  (3.2.2)  corresponds  to  a  flame  sheet  model  with  TC-I  being 
consumed  very  rapidly  at  a  standoff  distance  L£.  If  ve  have  a  distribution 
more  appropriate  to  our  particle  burning  model  in  the  form 


Y 


i 


l  o 


then  we  see  that 


0  1  x  i  L{ 
x  >  Lf 


L 


Y  L. 
o  f 


(3.2.4) 


(3.2.5) 


We  may  therefore  interpret  values  of  L*  as  corresponding  to  an  ideal  flane 
length  L£  =  6L* /Y0  if  the  distribution  (3.2.4)  applies. 

Our  subsequent  results  will  introduce  a  non-diaeni ionsl  flame  length 
L^d  which  is  simply  related  to  L*  by 


(3.2.6) 


where  Dx  i  s  the  diameter  of  the  bore. 
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Table  3.4  presents  the  relationship  between  muzzle  velocity  and  flame 
thickness  for  a  series  of  data  bases  corresponding  to  Table  3.1  Only  the 
composition  of  the  products  of  surface  decomposition  TC-F/TC-I  and  the 
particle  diameter  Dp  are  varied. 


Table  3.4  Relation  Between  Muzzle  Velocity  and  Flame  Thickness 
(TC-I  Burn  Rate  =  0.337p8****  cm/sec) 


TC-F/TC-I 

DP 

cm 

^D,  MAX 

Muzzle  Vel. 
m/  sec 

f-TQ 

Asm 

Ae« 

50/  50 

0.127 

0.40 

2889 

0.000 

3.37 

2.89 

50/50 

0.254 

1.52 

2603 

0.050 

-1.58 

-2.22 

50/  50 

0.3  81 

3  .36 

2502 

0,123 

-2.17 

-2.77 

50/50 

0.508 

5.56 

2461 

0.188 

-2.25 

-2.86 

7  5/25 

0.127 

0.14 

2906 

0.000 

0.46 

1.61 

75/25 

0.254 

0.45 

27  93 

0.007 

-2.21 

-1.3  9 

75/25 

0.381 

0.89 

2732 

0.021 

-2.45 

-1.69 

75/25 

0.508 

1.45 

2703 

0.037 

-2.58 

-1.86 

25/75 

0.127 

0.79 

2830 

0.000 

2.58 

1.23 

25/75 

0.254 

2.54 

2434 

0.131 

1.60 

-1.51 

25/75 

0.3  81 

6.07 

2249 

0.301 

-0.18 

-3.3  8 

25/75 

0.508 

8.42 

217  5 

0.3  89 

-0.91 

-4.17 

For  each  run  we  tabulate  TC-F/TC-I,  the  ratio  of  the  mass  fractions  of 
the  final  and  intermediate  TC  products  at  the  base  of  the  graveling  charge* 
Dp.  the  particle  diameter)  MAX’  *ke  maximum  value  of  which  occurs  in 

the  solution*  the  muzzle  velocity*  f— ' TCI,  the  value  of  the  mass  of  the 
unburned  intermediates  at  muzzle  exit  divided  by  the  initial  mass  of  the 
traveling  charge*  and  An ft,  Ae%,  the  final  mass  and  energy  defects. 

It  is  evident  in  Table  3.4  that  many  of  these  runs  involved  an 
appreciable  quantity  of  unburned  traveling  charge  intermediates  at  muzzle 
exit.  In  order  to  remove  the  complicating  influence  of  unburned  propellant 
from  the  relation  between  muzzle  velocity  and  flsme  length,  we  repeated  the 
matrix  of  runs  with  a  particle  burn  rate  calculated  to  ensure  complete 
burnout  of  the  intermediates.  These  results  are  shown  in  Table  3.5. 
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Table  3.5  Relation  Between  Muzzle  Velocity  and  Flame  Thickness 
(TC-I  Burn  Rate  =  50.8  cm/ sec) 


TC-F/TC-I 

DP 

cm 

lnd,  MAX 

Muzzle  Vel. 

m  /  sec 

f-TCI 

AsA 

Ae« 

50/  50 

0.127 

0.54 

2832 

0.000 

1.10 

0.12 

50/50 

0.254 

1.27 

2726 

0.000 

-1.6  8 

-2.3  8 

50/50 

0.3  81 

1.90 

2659 

0.000 

-2.10 

-2.77 

50/50 

0.508 

2.50 

2604 

0.000 

-2.17 

-2.82 

7  5/25 

0.127 

0.14 

2883 

0.000 

-0.37 

0.54 

7  5/25 

0.254 

0.33 

2  824 

0.000 

-2.23 

-1.41 

7  5/25 

0.3  81 

0.48 

27  92 

0.000 

-2.47 

-1.68 

7  5/25 

0.508 

0.64 

2766 

0.000 

-2.51 

-1.77 

25/75 

0  .127 

1.12 

2774 

0.000 

1 .25 

-0.51 

25/75 

0.254 

2.04 

266  5 

0.000 

0.95 

-2.56 

25/75 

0.381 

2.92 

2556 

0.000 

-0.27 

-3.71 

25/75 

0.508 

3.70 

245  9 

0.000 

-0.77 

-4.23 

The  results  of  Tables  3.4  and  3.5  are  presented  graphically  in  Figure 
3.4.  It  is  evident  that,  for  those  data  corresponding  to  complete 
combustion,  there  is  quite  a  good  correlation  between  muzzle  velocity  and 

^ND,  MAX’  *1 »°  indicate  the  muzzle  velocities  obtained  with  the  ideal 

thin  flame  traveling  charge  and  its  conventional  equivalent.  It  is  further 
evident  from  Tables  3.4  and  3.5  and  from  Figure  3.4  that  there  are  solutions 
with  values  of  Lj^  MAX  the  order  of  0.5  to  1.0  for  which  the  degradation 


of  traveling  charge  performance  enefit  is  quite  small.  Thus  we  conclude 
that  a  characteristic  flame  length  of  one  caliber  is  quite  compatible  with 
the  theoretical  benefit  of  the  traveling  charge.  Since  Equation  (3.2.4)  is 
a  reasonable  representation  of  the  distribution  of  the  intermediate  species 
and  since  Y0  corresponds  to  a  value  between  0.25  and  0.75  in  these 


calculations,  we  may  say  that  a  physical  combustion  zone  thickness  of  five 
to  ten  calibers  appears  tolerable  without  substantial  degradation  of 
performance . 
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Figure  3.4  Relation  Between  Muzzle  Velocity  end  Non-Dimens  mal  TC  Fiane  Tliickne** 


Ve  conclude  our  discussion  with  one  other  consent.  In  Table  3.6  we 
examine  the  effect  on  performance  of  the  time  of  ignition  of  the  traveling 
charge  for  the  case  TC-F/TC-I  =  50/50,  Dp  =  0.508  cat  in  Table  3.5.  The 

nominal  delay  is  1.8  msec  and  we  consider  a  range  from  1.4  to  2.4  msec.  It 
is  evident  that  for  this  problem  there  is  virtually  no  dependence  of  muzzle 
velocity  on  ignition  delay.  Accordingly,  the  loss  of  performance  due  to 
flame  thickness  cannot  be  compensated  by  a  change  in  ignition  delay,  at 
least  in  this  case. 


Table  3.6 

Effect  of  TC-Ignition  Delay 
(TC-F/TC-I  =  50/50.  Dp  =  0. 

on  Finite 
508,  B.  R.  = 

Plane  TC  Performance 
=  50.8  cm/ sec) 

TC-Ignition 
mse  c 

Delay  Muzzle  Velocity 

m/se  c 

^D,  MAX 

f-TQ 

Am* 

Ae* 

1.4 

2575 

2.28 

0.000 

-2.21 

-3  .08 

1.6 

2590 

2.43 

0.000 

-2.14 

00 

cs 

t 

1  .8 

2604 

2.50 

0.000 

-2.17 

-2.82 

2.0 

2605 

2.46 

0.000 

-2.23 

2.2 

2607 

2.39 

0.009 

-1.79 

2.4 

2592 

2.30 

0.044 

-0.87 

-2  ._ 
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NOIENCUIOBB 


Cross-sectional  urea  of  flow.  Equals  area  of  tube  minus 
intrusions  of  reactive  sidewalls,  nnbnrned  igniter  and 
projectile  afterbody. 

Born  rate  coefficient. 

Covoiume, 

Isentropic  sound  speed  in  mixture  of  combustion  products  at 
constant  composition. 

Specific  heat  at  constant  pressure. 

Specific  heat  at  constant  volume. 

Effective  diameter  of  particle. 

Surface  regression. 

Initial  diameter  of  perforation. 

Internal  energy.  Thermal  component  only. 

Heat  released  per  unit  mass  of  igniter. 

Heat  released  per  unit  mass  of  propellant  type  j. 

Heat  released  per  unit  mass  of  sidewall. 

Steady-state  interphaae  drag. 

Friction  force  between  traveling  charge  and  tube  wall. 
Constant  used  to  reconcile  units  of  measurement. 

Total  number  of  types  of  solid  propellant. 

Total  number  of  reactions  in  mixture  of  combustion  prcducts. 
Mass  of  projectile. 

Molecul  ar  weight. 

Rate  of  decomposition  per  unit  volume  of  propellant  type  j. 
Rate  of  decomposition  per  unit  volume  of  reactive  sidewall. 
Total  number  of  species  in  mixture  of  combustion  products. 


p 

Ok 


<iw 


*k 

*ik 


r 


te 


t 


o 


u 


p 


He 


V 


P 


*i 


Y.  . 

ij,  o 


TIGi 


z 


Born  rate  exponent. 

Pressure. 

Heat  released  per  unit  mass  by  reaction  k. 

Heat  loss  per  unit  volume  to  tube  wall. 

Heat  loss  per  unit  volume  to  propellant  type  j. 

Universal  gas  constant. 

Rate  of  reaction  k  per  unit  volume. 

Rate  of  production  of  species  1  by  reaction  k. 

Rate  of  regression  of  rear  surface  of  traveling  charge. 

Surface  area  of  a  grain. 

Temperature. 

Time. 

Velocity  of  mixture  of  ccsobustion  products. 

Velocity  of  solid  propellant. 

Velocity  of  traveling  charge. 

Volume  of  a  grain. 

Projectile  velocity. 

Rate  of  deposition  of  condensed-pha se  species  i  onto  surface  of 
solid  propellant. 

Mass  fraction  of  species  i  in  mixture  of  combustion  products. 

Ma  ss  fraction  of  species  i  in  near  field  products  of  combustion  of 
propellant  type  j. 

Mass  fraction  of  species  i  in  products  of  combustion  of  igniter. 
Mass  fraction  of  reactive  sidewall. 

Axial  coordinate. 
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Greek  Symbol  s 


r 

e 

$ 

P 

PCi 

PIG 

Ps 

Ptc 

o 

c 


Ratio  of  specific  heats. 

Porosity  or  fraction  of  unit  volume  occupied  by  mixture  of 
combustion  products. 

No  n-homoge  neons  term  in  balance  equations. 

Density  of  mixture  of  combustion  products. 

Density  of  condensed  species  in  mixture  of  combustion  products. 
Density  of  unbuxned  igniter 
Density  of  sidewall. 

Density  of  traveling  charge. 

Intergranular  stress. 

Stress  in  traveling  charge 
Rate  of  decomposition  of  igniter. 


Special  Symbols 

—  Convective  derivative  along  streamline  of  mixture  of  combustion 

Dt  products. 


-  Convective  derivative  along  solid-propellant. 

Dt 

P 
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APPENDIX: 

XNOVAKTC  (XKTC)-STRUCTURE  AND  USE 


65 


INTENTIONALLY  LEFT  BLANK. 


66 


Oar  intention  in  this  Appendix  is  to  provide  the  user  of  the  code  with 
some  understanding  of  its  macro  structure  and  full  details  of  the  specifica- 
tion  of  input  data.  Ve  do  not  attempt  to  document  the  code  in  such  detail 
as  to  permit  revisions  by  the  user.  XXT C  is  an  extension  of  N OV AT C ,  XNOVAT 
and  XNOVAK  which  are,  in  turn,  extensions  of  the  express  version  of  the  NOVA 
Code  known  as  XNOVA.  XNOVA  is  a  subset  of  the  NOVA  Code  and  certain  of  the 
details  of  the  original  versions  of  the  NOVA  Code  which  were  deleted  in  the 
preparation  of  XNOVA  are  also  absent  from  its  successors  XNOVAK  and  XNOVAT. 
Discussions  of  the  structure  of  NOVA,  NOVATC,  XNOVA,  XNOVAK  and  XNOVAT  may 
be  found  in  earlier  reports.1 '3 '4.6,2  Although  we  provide  a  complete 
tabulation  of  the  present  set  of  code  subroutines  and  functions,  our  general 
discussion  is  confined  to  the  differences  between  XKTC  and  its  predecessors. 

The  discussion  of  structure  and  of  input  is  contained  in  the  two 
following  Sections.  Before  proceeding,  however,  we  wish  to  draw  the  user’s 
attention  to  certain  general  restrictions  on  the  use  of  the  code. 

First,  it  should  be  noted  that  existing  NOVA,  XN0VA  and  NOVATC  data 
bases  cannot  be  read  by  XNOVAK,  XNOVAT  or  XKTC.  Although  we  always  attempt 
to  maintain  data  base  compatibility  between  the  various  code  versions, 
XNOVAK,  XNOVAT  and  XKTC  do  require  a  minor  change  to  pre-existing  data 
bases.  File  [M4] ,  described  in  Table  A. 5,  always  consists  of  two  cards  in 
XNOVAK/ XNOV  AT /XKTC  data  sets.  Pre-existing  NOVA,  XNOV  A  and  NOVATC  data  sets 
in  which  File  [M4]  consisted  of  jnst  one  card  may  be  made  compatible  with 
XNOVAK/XNOVAT/XKTC  by  the  incorporation  of  a  single  blank  card  following  the 
pre-existing  File  [M4] .  No  change  is  required  to  pre-existing  data  bases 
for  which  [M4]  already  consisted  of  two  cards.  A  corresponding  revision  to 
XNOVAK/XNOVAT/XKTC  data  bases  is  required  if  they  are  to  be  run  on  NOVATC, 
XNOVA  or  NOVA.  However,  all  XNOVAK  data  bases  can  be  read  by  XNOVAT  and 
XKTC. 


Second,  in  regard  to  the  summarized  solution  histories  given  at  the 
conclusion  of  each  run,  it  should  be  noted  that  the  energy  defect  calcula¬ 
tion  does  not  presently  support  the  full  chemistry  option.  If  intemodiate 
combustion  products  are  present,  the  tabulated  energy  defect  may  become 
quite  large.  Only  the  mass  defect  can  be  depended  npon  to  gage  the  accuracy 
of  the  numerical  solution  in  such  cases.  The  single  exception  to  this 
statement  is  that  of  the  traveling  charge  with  a  finite  flame  thickness.  In 
this  special  case,  where  there  is  just  one  intermediate  species,  the  energy 
defect  calculation  reflects  the  chemical  energy  stored  in  the  intermediates. 
Moreover,  in  this  special  case,  the  mass  fraction  calculation  for  propellant 
type  three  in  the  summary  table  is  used  to  tabulate  the  ratio  of  the  mass  of 
the  unburned  intermediates  to  that  of  the  original  traveling  charge 
increment. 
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STRUCTURE  AND  LINKAGES 


A  complete  summary  of  the  routines  and  their  linkages  is  given  in 
Table  A.l.  The  overall  macro  structure  is  unchanged  from  that  of  XNOVA  os 
documented  in  Reference  4.  The  following  brief  comments  provide  a  summary 
of  the  routines  which  were  added  to  XNOVA  to  produce  XNOVAK,  the  revisions 
which  subsequently  transformed  XNWAK  into  XNOVAT,  and  the  final  asalgam  of 
XNOVAT  and  BRLTC^  to  produce  XETC.  It  should  be  noted  that  we  have 
previously  commented  on  the  fact  that  XNWAK  was  written  on  a  32-bit  word 
machine  in  double  precision, **  and  that  users  who  implement  the  Code  on  a  60- 
bit  machine  such  as  the  CYBQl  7600  might  wish  to  convert  the  Code  to  single 
precision.  The  version  of  XNOVAK  used  to  develop  XXTC  was  implemented  on 
the  CYBER  at  BRL  as  a  single  precision  code.  Accordingly,  XKTC  is  a  single 
precision  code  and  it  follows  that  implementation  of  XKTC  on  a  32-bit  word 
machine  should  incorporate  a  conversion  to  double  precision. 


Routines  Added  in  Preparation  of  XNOVAK 

XNOVAK  was  developed  as  an  extension  of  XNOVA.  Accordingly,  XNOVA  is  a 
subset  of  both  XNOVAK  and  NOVA.  Moreover,  in  developing  XNOVAK  it  was  found 
necessary  to  restore  certain  capabilities  which  had  been  removed  from  NOVA 
during  the  preparation  of  XNOVA.  Specifically,  we  restored  the  invariant, 
embedding  solution  for  the  thermal  response  of  the  solid  propellant.  Thus, 
in  addition  to  the  capabilities  of  XNOVA,  XNOVAK  and  NOVA  also  share  a 
transient  combustion  modeling  capability. 

The  NOVA  routines  restored  to  XNOVAK  were  CMBUST,  FBAK,  IKBED  and  ROOT 
together  with  necessary  I/O  linkages.  The  routines  were  then  modified  to 
support  the  XNOVAK  features  of  subsurface  reactivity  and  the  evaporative 
boundary  condition.  New  routines  added  to  support  the  kinetics  option 
included  QIEMR,  CBEMRS,  CONLOS,  COVC  and  GAMMOL.  The  new  routine  FILIBR 
implements  the  analysis  of  the  rear  filler  elements  and  was  created  by 
making  appropriate  editing  changes  to  the  pre-existing  routine  FILLER.  One 
additional  routine,  AVN,  was  included  for  the  purpose  of  documenting  Code 
revisions  and  printing  the  version  number  on  each  run. 


Routines  Added  in  Preparation  of  XN WAT 


Several  routines  were  required  to  produce  XNWAT,  Subroutine  INTO  PC 
determines  the  rate  of  combustion  of  the  reactive  layers  ascribed  to  the 
tube  wall,  the  centerline,  and  the  projectile  afterbody.  To  support  the  new 
form  functions  we  have  added  the  routines  BLSL,  HEX19,  and  SLIVER.  To  treat 
the  reactivity  and  flow  resistance  of  the  enthralls  of  the  bundles  of  pro¬ 
pellant  we  added  the  routines  CALFLO,  CALPRM,  FLOX,  FLDXER  AND  TDBCAL,  all 
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of  which  were  adapted  from  the  TDNOVA  Code.^  The  routine  TEBUF  acts  as  a 
link  between  the  existing  boundary  value  solver  BCAL  and  the  routine  TEBCAL. 
In  addition  we  added  the  routines  JOON,  SETSUB  and  RCASE  which  perform 
various  utility  functions. 

Structure  of  NOVATC 


XKTC  was  essentially  formed  by  substituting  IN0VAT  in  place  of  XNOV^  in 
the  previously  developed  NCVATC  Code.^  As  we  have  previously  discussed. 
NOVATC  is  an  amalgam  of  INOVA  and  BRLTC.*'® 

In  forming  the  amalgam  care  was  taken  to  minimize  cross-contamination 
of  the  coding.  The  approach  was  to  take  BELTC  and  delete  all  references  to 
the  flow  behind  the  base  of  the  traveling  charge  while  retaining  the 
storage,  data  reads  and  processing  pertinent  to  the  regression  of  the 
surface,  the  traveling  charge  and  the  projectile.  This  stripped  down  code 
was  linked  to  a  complete  version  of  XNOVA’  in  which  a  minimal  number  of 
linking  subroutine  calls  were  placed.  Although  a  more  tightly  written  code 
could  have  been  developed  it  would  have  been  more  difficult  to  maintain  such 
a  code  in  a  form  compatible  with  future  versions  of  INOVA  and  BELTC.  Also, 
it  is  the  case  that  INOVA  and  BELTC  have  many  variables  whose  Fortran  names 
are  identical,  but  whose  meanings  differ  and,  conversely,  some  variables 
with  different  names  but  identical  meanings. 

It  will  be  noted  in  Table  A.l  that  all  routines  pertinent  to  the 
traveling  charge  have  the  prefix  T C.  A  single  routine,  TCLINE,  has  the 
function  of  reconciling  the  storage  conventions  of  the  two  codes  and  of 
passing  data  between  them.  Processing  follows  the  usual  NOVA  conventions 
with  the  traveling  charge  logic  activated  at  a  number  of  places.  NWSDB 
calls  TCDATA  and  TCLINE  to  read  and  initialize  data.  INTEG  calls  TCLINK, 
TCdJT  to  print  the  solutions  and  TCBR4  to  determine  the  traveling  charge 
time  step  constraint.  INT/vL  calls  TCLINK  and  TCXC,  the  main  integration 
executive  for  the  traveling  charge.  TABLES  calls  TCLINK.  It  was  also 
necessary  to  restructure  slightly  the  routines  BCAL  and  BDMOV. 

Routines  Added  in  Preparation  of  IKTC 


The  approach  followed  in  the  development  of  NOVATC  made  the  replacement 
of  INOVA  by  INOVAT  fairly  straightforward.  It  should  be  noted,  however, 
that  the  present  effort  involved  substantial  exercising  of  the  algorithm  and 
several  revisions  were  made  to  the  analysis  of  the  boundary  values  at  the 
base  of  the  burning  traveling  charge.  Accordingly,  the  versions  of  TCLINK 
and  TQBASE  used  here  differ  fiom  those  in  earlier  versions  of  NOVATC. 
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Apart  from  linking  the  kinetics,  tank  gun  and  traveling  charge  options 
in  a  physically  complete  fashion,  it  was  necessary  to  add  several  new 
routines  in  the  development  of  XKTC.  Subroutine  GETLS  computes  the  values 
of  L*  for  the  intermediate  products  of  combustion  of  the  traveling  charge. 
ATACH  enforces  the  attachment  of  a  propellant  type  to  the  tube  or  to  the 
projectile.  Subroutine  LPBV  provides  an  update  of  the  combnstion  chamber  of 
a  control  charge.  It  is  supported  by  FORMCR  which  computes  the  form 
functions,  and  by  the  previously  developed  NOVA  routine  FRATE  which  computes 
mass  transfer  rates. 


DESCRIPTION  OF  INPUT  FILES 


Ve  preface  the  detailed  description  of  the  input  file  structure  with 
some  general  observations.  We  note,  in  particular,  the  nomenclature  for  the 
definition  of  the  geometry  of  the  propelling  charge. 

Figure  A.  1  illustrates  &  charge  configuration  which  does  not  make  use 
of  the  tank  gun  option  (MDDET  =0).  We  illustrate  th*  base  of  the  pro¬ 
jectile,  an  arbitrary  number  of  filler  elements  —  some  at  the  rear  and  some 
at  the  front  of  the  charge  —  and  two  types  of  propellant,  the  first  of 
which  is  in  two  increments,  the  latter  increment  overlapping  that  of  the 
second  propellant.  The  figure  also  illustrates  the  significance  of  the 
input  data  ZGR,  XEL,  XBL  and  ZBPR, 

Figure  A. 2  illustrates  a  charge  configuration  which  does  make  use  of 
the  tank  gun  option  (M3DET  «  1).  We  illustrate  the  intrusion  of  the 
projectile  afterbody  into  the  combustion  chamber,  the  reactive  layers  on  the 
tube  wall,  the  centerline  and  the  projectile  afterbody,  and  three  types  of 
propellant,  the  first  of  which  is  present  as  two  increments.  To  the  rear  of 
the  afterbody  we  represent  the  three  types  of  propellant  as  being  parallel 
packaged  and  we  illustrate  the  significance  of  the  input  data  RGRI  and  RGRO. 

Figures  A.l  and  A. 2  both  assume  that  the  propellant  is  of  the 
conventional  type  Figure  A. 3  illustrates  a  charge  configuration  in  which  a 
traveling  charge  is  present.  We  note  that  the  traveling  charge  may  consist 
of  several  increments,  each  having  distinct  mechanical  and  thexmochemical 
properties.  The  simulation  of  the  traveling  charge  may  be  performed  in 
conj unction  wi th  an  arbitrarily  configured  conventional  charge.  It  is 
assumed,  however,  that  when  the  traveling  charge  is  present  there  are  no 
compa ctibl e  filler  materials  at  the  front  of  the  chamber,  as  in  Figure  A.l, 
and  that  the  projectile  does  not  have  an  afterbody  as  in  Figure  A. 2. 

Reactive  sidewalls  may  nevertheless  be  present,  as  shown  in  Fignre  A. 3. 
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nomenclature  for  Definition  of  Charge  Configuration  in  .'AT*  with  MODET 


Input  to  XKTC  consists  of  XNOVAT  data  followed  by  a  subset  of  the  BHLTC 
data  base  when  the  traveling  charge  option  is  selected  in  the  control  data 
field  of  me.  Tables  A. 2,  A. 3  and  A. 4  enumerate  the  "Mandatory", 

Contngent  and  Traveling  Charge  data  files.  The  files  prefixed  with  an  M 
are  always  mandatory  in  XNOVAT,  Here,  File  [M13]  is  optional.  We  have 
retained  the  file  labels  of  XNOVAT,  in  spite  of  this  minor  inconsistency,  to 
promote  commonality  among  the  various  codes. 

The  enumeration  of  the  M  and  C  files  follows  that  for  the  NOVA  Code. 

We  have  attempted  to  maximize  the  compatibility  of  XKTC,  XNOVAT,  XNOVAK. 
XNOVATC,  XNOVA  and  NOVA  data  bases.  As  noted  in  the  introduction  to  this 
appendix,  compatibility  is  complete  with  the  possible  exception  of  File 
[M4]  .  In  XKTC,  XNOVAT  and  XNOVAK,  File  [M4]  always  consists  of  two  cards. 

In  XNOVATC,  XNOVA  and  NOVA  it  may  consist  of  either  one  or  two  cards. 
Existing  NOVATC,  XNOVA  and  NOVA  data  sets  for  which  File  [M4]  consists  of 
only  one  card  may  be  made  compatible  with  XKTC,  XNOVAT  and  XNOVAK  by  the 
insertion  of  a  single  blank  card  after  the  existing  [M4]  data  card.  NOVA 
files  not  supported  by  XKTC,  XNOVAT  and  XNOVAK  are  read  and  subsequently 
ignored.  The  tabulation  of  data  by  XKTC  only  includes  those  elements  which 
are  supported. 

All  mandatory  files.  Table  A. 2,  are  common  to  all  codes.  However,  not 
all  the  members  of  every  file  are  supported  by  XKTC,  XNOVAT  and  XNOVAK.  The 
detailed  definition  of  the  input  files  given  in  Table  A. 5  describes  all  the 
elements  which  are  supported  and  simply  denotes  the  others  as  "Inactive." 
Those  contingent  files.  Table  A. 3,  which  are  supported  are  defined  while  the 
others  are  denoted  as  "not  supported."  Only  those  contingent  files  which 
are  supported  by  XKTC  appear  in  the  detailed  discussion  of  Table  A. 5. 

The  p.ogrrai  input  files  are  structured  to  permit  the  running  of  only 
one  problem  at  u  time.  A  problem  may  be  restarted  from  disc,  provided 
suitable  options  have  been  selected  in  a  prior  run.  Each  problem  may  be 
terminated  according  to  a  criterion  of  number  of  integration  steps,  pro¬ 
jectile  displacement  or  problem  time.  Termination  occurs  at  whichever  of 
the  three  criteria  is  first  satisfied. 


It  should  be  noted  that  file  counters  are  edited 
to  establish  their  conformality  with  internal  storage 
other  data  are  not  edited,  and  it  is  recommended  that 
perform  a  first  run  on  u  new  problem  with  termination 
set-up,  so  that  all  data  may  be  validated. 


by  subroutine  NOV  SUB 
dimensions.  However, 
the  program  nser 
set  to  follow  initial 


The  following  notes  (1)  -  (10)  apply  whether  or  not  the  traveling 
charge  option  is  exercised.  Subsequent  notes  (11)  -  (16)  pertain  to  the 
traveling  charge  option, 

( 1)  For  problems  involving  only  granular  propellant,  the  representations  of 
the  grain  distributions  are  arbitrary.  The  granular  propellants  may  doubly 
or  triply  overlap.  However,  a  region  occupied  by  stick  propellant  may  not 
be  shared  by  any  other  type  unless  MODET  =  1.  IF  HOIST  =  1,  it  is 
permissible  to  have  an  increment  of  stick  propellant  parallel  packaged  with 
another  increment  of  stick  or  granular  propellant.  The  monolithic  charge  is 
not  permitted  to  share  a  cross-section  of  the  tube  with  any  other  type.  If 
a  given  type  is  specified  as  bonded  to  the  tube  or  projectile,  the  bonding 
will  be  presumed  to  apply  to  any  other  propellant  which  shares  a  region  with 
the  given  type. 

(2)  The  mode  of  representation  of  one  type  of  propellant  is  independent  of 
the  mode  of  representation  of  any  other.  If  tables  of  volume  fractions  are 
used,  it  should  be  recalled  that  the  volume  fraction  is  the  complement  of 
the  porosity,  being  zero  when  no  propellant  is  present  and  increasing  to 
unity  as  the  pore  volume  vanishes.  The  volume  traction  defaults  to  zero 
outside  the  range  of  positions  defined  by  the  tables. 

(3)  An  internal  boundary  is  identified  at  any  location  in  which  any  volume 
fraction  is  discontinuous.  Thus  an  internal  boundary  may  be  established  at 
a  point  at  which  the  porosity  is  continuous. 

(4)  Absolutely  no  restrictions,  other  than  that  on  note  (1),  are  placed  on 
the  relative  positions  of  Jie  bags  of  propellant.  After  the  data  are  read, 
they  are  all  converted  internally  to  volume  fraction  tables.  These  are  all 
then  scanned  to  create  a  list  of  internal  boundary  positions.  This  list  is 
then  sorted  to  put  them  into  order  and  to  define  the  computational  regions. 

(5)  The  right-hand  boundary  is  always  set  equal  to  XEL(  1) .  If  NEL  =  0, 

XEL(  1)  defaults  to  ZBFR.  Similarly,  the  left-hand  boundary  is  always  set 
equal  to  XBL(l).  If  NBI.  =  0,  XBL(l)  defaults  to  zero.  It  should  be  noted 
that  the  configuration  of  the  forward  elements  is  specified  by  reference  to 
the  position  of  the  left-hand  or  rear  boundary  of  each  element.  Conversely, 
the  conf igure tion  of  the  rear  elements  is  specified  by  reference  to  toe 
position  of  the  right-hard  or  forward  boundary  of  each  clement.  In  both 
cases  the  elements  are  ordered  so  that  the  first  is  closest  to  the 
propelling  charge, 

(6)  Although  filler  elements  may  be  included  if  MDPET  =  1,  the  coding  does 
not  consider  the  influence  of  the  afterbody  on  the  behavior  of  the  filler 
el  cme  nt  s. 
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(7)  It  is  emphasized  that  every  case  closure  element  must  be  accounted  for, 
including  spaces.  Any  element,  other  than  the  first,  may  be  entered  as  a 
space  by  setting  its  mass  equal  to  zero. 

(8)  It  is  noted  that  one  may  use  a  single  element  to  represent  more  than 
one  physical  component  and  that  conversely,  one  may  use  several  elements  to 
characterize  a  single  component. 

(9)  The  present  storage  allocation  in  the  code  imposes  certain  restrictions 
on  the  modeling  of  subsurface  reactions  in  the  propellant.  It  is  assumed 
that  only  one  type  of  propellant  is  present  if  an  invariant  embedding 
solution  of  the  thermal  response  is  desired  (DIODE  =  2  in  File  [M2]).  At 
most  21  stations  are  permitted  for  the  analysis  of  the  thermal  response.  At 
most  two  chemical  species  can  be  associated  with  the  subsurface  chemistry. 
While  these  limits  are  thought  to  be  adequate  for  present  needs,  they  may 
readily  be  relaxed  by  means  of  changes  in  certain  dimension  statements. 

(10)  It  will  be  noted  that  the  data  required  when  KMODE  ^  0  Files  ([C23.1]  - 
[C23.6])  are  to  some  extent  already  provided  by  the  data  of  Files  [M7]  - 

[C7] .  The  latter  data  are  superseded  in  such  a  case  and  consistency  between 
the  two  sets  of  data  is  not  essential.  However,  inconsistencies  may  yield 
an  apparent  energy  defect  in  the  summarized  interior  ballistics  table. 

(11)  Onit  9  is  reserved  for  storage  of  the  traveling  charge  data  if  disk 
logout  is  specified. 

(12)  The  code  does  not  support  the  traveling  charge  option  simultaneously 
with  the  compactible  case  closure  option.  It  is  also  assumed  that  the 
projectile  has  no  afterbody  if  the  traveling  charge  option  is  exercised. 

(13)  The  code  termination  data  are  always  specified  according  to  the  NOVA 
format.  The  projectile  data  follow  the  NOVA  format  only  if  the  traveling 
charge  option  is  not  in  effect. 

(14)  Due  to  the  cross-linking  of  the  codes,  no  more  than  98  mesh  points  and 
no  more  than  9  increments  can  be  specified  for  the  booster  region  when  the 
traveling  charge  option  is  in  effect. 

(15)  A  composition  dependent  covolume  is  not  considered.  The  covolume  for 
the  calculation  is  presently  based  on  the  booster  charge  data  according  to 
the  NOVA  convention. 

(16)  The  gun  tube  is  assumed  to  have  a  constant  area  in  the  region  occupied 
by  the  unreacted  traveling  charge. 
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TABLE  A,  1  SUMMARY  OF  ROUTINES  AND  LINKAGES 


XKTC 

Purpose  : 

XKTC  is  a  dummy  maiu  routine  which,  immediately 
transfers  control  to  NOVSUB. 

Calls: 

NOV  SUB. 

Called  by: 

None . 

ATACH 

Purpo  se  : 

Subroutine  ATACH  enforces  the  attachment  of  a  charge 
increment  to  either  the  tube  or  the  projectile. 

Calls  : 

None . 

Called  by: 

INTAL. 

AVN 

Purpose  : 

Subroutine  AVN  is  used  to  document  code  revisions 
and  print  the  version  number. 

Calls: 

None  . 

Called  by: 

NOV  SUB. 

“V 

BCAL 

Purpose  : 

Subroutine  BCAL  accepts  trial  boundary  values  for 
the  gas-  and  solid-phases  and  adjusts  them  so  as  to 
satisfy  simultaneously  the  physical  and  character¬ 
istic  boundary  conditions. 

Calls  : 

COVC,  EPTOR,  GAMMOL,  LUMPS,  OUTPUT,  PRTDC, 

PR  TOE,  TABLES,  TDBUF. 

Called  by: 

INTAL. 

BDMOV 

Purpo  se : 

Subroutine  BDMOV  updates  the  position  of  the  region 
boundaries.  It  also  updates  the  projectile 
displacement  and  velocity,  making  use  of  FILLED 
and/or  FILLER  when  esse  closure  elecontr.  are  present. 

Calls : 

FILLER,  FILLER,  RCASE,  RES  FUN,  EBEt. 

Called  by: 

INTAL. 

BLED  AT 

Purpose  : 

BLKDAT  performs  block  dsts  ini tisl izstion. 

Calls: 

No  ne , 

Called  by: 

None . 
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BLSL 

Purpo  se  : 

Subroutine  BLSL  computes  the  volume  and  surface 
area  for  the  blind  slit  form  function  during  the 
slivering  phase 

Cells: 

None. 

Called  by: 

FORM. 

RRNOUT 

Purpo  se  : 

Subroutine  BRNOUT  detects  the  instant  when  all  the 
propellant  has  been  consumed  and  then  consolidates 
all  the  computational  regions  into  a  single  region 
with  ur.ifonn  mesh  spacing. 

Calls: 

FITSP,  KBG21 . 

Called  by: 

None . 

CAL  FLO 

Purpo  se : 

Subroutine  CAL  FLO  computes  the  rate  of  reaction  of 
the  surface  fluxes  associated  vith  the  endwalls  of 
the  charge  increments. 

Calls: 

None . 

Called  by: 

TDBUF. 

CAL  PR  M 

Purpo  se : 

Subroutine  CALPRM  computes  the  flow  resistance  co¬ 
efficient  for  the  endwalls  of  the  charge  increments. 

Calls  : 

None. 

Called  by: 

TDBUF. 

CHEMR 

Purpo  se  : 

Subroutine  CHEMR  computes  the  effects  of  chemical 
reactions  on  the  species  mass  fractions  in  the 
combustion  product  mixture  and  returns  the  value  of 
the  heat  liberated  by  chemical  reaction. 

Call s: 

None . 

Called  by : 

INTAL,  LUMPS. 
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CHEMRS 


CMBUST 


CONLOS 


COVC 


DINT 


DPDER 


Purpo  se : 

Subroutine  CHEMRS  performs  a  function  similar  to  ihat 
of  CHEMR)  however  it  addresses  the  reactions  within 
the  solid  propellant. 

Calls: 

None . 

Called  by: 

IMBED. 

Purpo  se  : 

Subroutine  CMBUST  interfaces  the  invariant  embedding 
solution  of  the  solid  propellant  thermal  response 
with  the  macroscopic  two-phase  flow. 

Calls: 

FBAK,  IMBED. 

Called  by: 

INTEG ,  NOV SUB . 

Purpo  se : 

Subroutine  CONLOS  computes  the  rate  of  transfer  of 
condensed  species  from  the  combustion  product  mixture 
to  tne  surface  of  the  solid  propellant. 

Call s: 

None . 

Called  by: 

HEATP. 

Purpose  : 

Function  COVC  computes  the  covolume  of  the  combustion 
product  mixture. 

Calls: 

None. 

Called  by: 

BCAL,  FRICT,  HEATP,  HEATW,  INTAL,  INTEG.  INTRPC, 
LUMPS.  REG12 ,  REG32,  STATES,  TCL INK,  TDBCAL,  TDBUF. 

Purpose : 

Function  DINT  perfoncs  a  truncation  of  a  floating 
point  variable  to  its  integer  part. 

Calls: 

None . 

Called  by: 

TABLES. 

Purpo  se : 

Function  DPDER  computes  the  partial  derivative 
(dp/3e)p  for  the  covolume  equation  of  state. 

Calls  : 

None . 

Called  by: 

LUMPS. 
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DFDRE 

Pur  po  se  : 

Function  DPDRE  computes  the  partial  derivative 
Op/dp)  e  for  the  covo)  mce  equation  of  state. 

Call s : 

None . 

Called  by: 

LUMPS. 

EPTOR 

Fur po se : 

Function  EPTOR  evaluates  the  density  as  a  function 
of  internal  energy  and  pressure  for  the  covolume 
equation  of  state. 

Calls: 

None . 

Call ed  by : 

BCAL,  FRICT.  HEATP,  HEATS,  INTRPC,  NOV  SOB,  TIB  CAL. 

ERTOP 

Pur  po  se  : 

Function  ERTOP  evaluates  the  pressure  as  a  function 
of  internal  energy  and  density  for  the  covolume 
equation  of  state. 

Calls: 

None . 

Call ed  by : 

INTAL,  LPVB,  LUMPS,  NOVSUB,  REG12,  RBG32. 

FBAK 

Purpose : 

Function  FBAK  computes  the  temperature  gradient 
at  the  surface  of  the  solid  propellant  taking  into 
account  heat  transfer  from  the  igniter  and  heat 
feedback  from  the  flame. 

Calls: 

None. 

Called  by: 

CMBUST,  NOVSUB. 

FXLLBR 

Purpose : 

Subroutine  FILLBR  updates  the  motion  of  the  rear  ca: 
closure  elements. 

Calls: 

OUTPUT,  TABLES. 

Called  by: 

BDMOV,  NOVSUB. 

FILLER 

Purpose  : 

Subroutine  FILLHl  updates  the  motion  of  the  forward 
case  closure  elements  and  the  projectile. 

Calls: 

OUTPUT,  TABLES. 

Called  by: 

BDMOV,  NOVSUB. 
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FITSP 


FLUX 


FLUXDR 


FORM 


FORMCR 


FRATE 


Purpose  ;  Subroutine  FITSP  sets  values  of  the  state  variables, 
following  re-allocation  of  the  mesh,  by  means  of  a 
cubic  spline  interpolation  scheme. 

Calls:  SETSUB. 

Called  by:  BRNOUT,  1CSH,  NOV  SUB. 


Purpo se  :  Subroutine  FLUX  computes  the  state  of  gas  at  the 

boundaries  of  the  reactive  endwalls. 

Calls :  None. 

Called  by:  FLUXDR,  TDBCAL. 


Pm po se  :  Subroutine  FLUXDR  computes  the  derivatives  of  the 

state  variables  in  a  lumped  parameter  region  with 
respect  to  the  fluxes  to  that  region  from  the 
contiguous  continuum  regions. 

Calls:  FLUX. 

Called  by:  TEBUF. 


Purpose  :  Subroutine  FORM  computes  the  surface  area  and  volume 

of  the  propellant  grains.  It  is  supported  by  PERF19 
which  treats  the  slivering  phase  of  nineteen 
perforation  propellant. 

Calls:  BLSL,  HEX19,  PERF19. 

Called  by:  IN  TEG,  NOV SUB. 


Purpose  : 

Calls : 
Called  by: 

Purpose  : 

Calls: 
Called  by: 


Subroutine  FCRMCR  computes  the  form  functions  for 
grains  used  to  define  a  control  charge. 


None . 
LPVB. 


Function  Frate  computes  the  rate  of  mass  transfer 
through  a  nozzle  subject  to  the  assumption  of 
isentropic  flow. 

None  . 

LPVB. 
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2RICT 

Purpose  : 

Subroutine  FRICT  computes  the  interphase  dreg  for 
granular  propellent. 

Cells: 

COVC,  EPTOR,  VIS. 

Celled  by: 

INTEG. 

GAMMOL 

Purpose  : 

Subroutine  GAMMOL  computes  the  ratio  of  specific 
heats  and  the  molecular  weight  of  the  combustion 
product  mixture. 

Calls: 

None . 

Called  by: 

BCAL,  LUMPS,  STATES,  TCLINK,  TUB CAL,  TDBUF. 

GEILS 

Purpose  : 

Subroutine  GETLS  computes  L*  values  for  the 
intermediate  combustion  products  of  the  traveling 
charge. 

Calls: 

None . 

Called  by: 

INTAL,  INTEG. 

HEATP 

Purpoi^ : 

Subroutine  HEATP  determines  the  interphase  heat 
transfer,  updates  the  surface  temperature  of  the 
solid-phase  prior  to  ignition,  and  also  determine 
the  interphase  drag  in  stick  propellant. 

Calls: 

COJtOS,  COVC,  EPTOR,  VIS. 

Called  by: 

INTEG. 

HEATW 

Purpose  : 

Subroutine  HEA1R  determines  the  rate  of  heat  tran 
to  the  tube  wall  and  updates  the  tube  surface 
temperature. 

Cells: 

CWC,  EPTOR,  VIS. 

Called  by: 

INTEG,  NOV SUB, 

HEX19 

Purpose : 

Subroutine  HRX19  computes  the  volume  and  surface 
area  of  nineteen-perforation  hexagonal  grains. 

Call s : 

SLIVER. 

Ca  1 1  e  d  by  : 

FORM. 
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IGNITE 

Pur  po  se : 

Subroutine  IGNITE  determines  the  time  of  ignition 
of  the  propellant. 

Call s : 

None . 

Called  by: 

INTEG. 

IMBED 

Purpose : 

Subroutine  IMBED  integrates  the  thermal  profile  in 
the  solid  propellant  by  means  of  the  method  of 
invariant  embedding. 

Calls: 

CHEMRS,  ROOT,  ROOTS. 

Called  by: 

CMBUST,  NOV  SOB. 

INTAL 

Pur po se  : 

Subroutine  INTAL  updates  the  principal  state 
variables,  p,  p,  u,  a,  u^,  or  at  all  interior 
mesh  points.  It  also  determines  trial  update 
boundary  valnes  together  with  characteristic 
coefficients  which  are  transmitted  to  BCAL  for 
imposition  of  the  physical  boundary  conditions. 

Calls: 

ATACH,  BCAL,  BDMOV,  CHEMR,  COVC,  ERTOP,  GETLS, 

INTR.P1 ,  LPVB,  FRTDC,  PRIDE.  RCASE,  RDEDPR,  IBEO. 
TCLINK,  TCXC. 

Called  by: 

INTEG. 

INTEC 

Pur  po  se : 

Subroutine  INTEG  is  the  principal  integration 
executive.  It  is  cycled  twice  per  time  step,  once 
for  the  predictor  level  and  once  for  the  corrector 
level.  It  establishes  the  constitutive  data  by  calls 
to  individual  subroutines  and  effects  the  update  of 
the  principal  state  variables  by  a  call  to  INTAL. 
INTEG  itself  updates  the  surface  regression  and 
thermal  parameter  H  for  the  solid-phase  and,  if 

EMODE  =0  (File  [M2]  of  Table  A. 4) ,  it  also  updates 
the  molecular  weight  and  ratio  of  specific  heats  for 
the  gas  phase. 

Calls  : 

ClfiUST,  COVC,  FORM,  FRICT,  GETLS,  HEATP,  BEATW, 
IGNITE,  INTAL,  INTRPC,  IN'PIPl,  INTRP2,  Z  CON,  LESS, 
OUTPUT,  PR  TOC,  REGRES,  STAr-^S,  TABLES,  TCBR4,  TCLINK, 
TCOUT,  TERMIN,  VOIDS,  WALTEH. 

Call ed  by : 

NOV  SUB. 
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INIRPC 

Pur  po  se  : 

Subroutine  INIRPC  computes  the  local  mass  addition 
due  to  combustion  of  the  reactive  sidewalls. 

Calls: 

COVC,  EPTOR,  JCON,  RCASE,  VIS. 

Called  by: 

INTEG,  NOV  SUB. 

INTO  PI 

Purpo  se : 

Subroutine  INTOP1  determines  the  cross-sectional  area 
of  the  tnbe  at  each  mesh  point,  including  an 
allowance  for  the  unburned  portion  of  the  primer. 

• 

Calls: 

JCON. 

- 

Called  by: 

INTAL,  INTEG,  NWSUB,  REG12,  REG21 ,  REG32,  TABLES. 

INTOP2 

Pur  po  se  : 

Subroutine  INTOP2  computes  the  rate  of  discharge  of 
the  igniter  at  each  mesh  point. 

Calls : 

None. 

Called  by: 

INTEG,  NOV  SUB,  REG12,  REG32. 

J  CON 

Pur  po  se  : 

Function  JCON  locates  the  nearest  continuum  mesh 
point . 

Calls: 

None  , 

Called  by: 

INTEG,  INIRPC,  INTO PI,  LUMPS. 

LPDATA 

Purpo  se  : 

Subroutine  LPDATA  converts  the  XNOVA  data  base 
into  an  equivalent  data  base  for  a  lumped  parameter 
interior  ballistic  model. 

Calls: 

N  one . 

Call ed  by: 

NOV  SUB. 

LPVB 

Purpo  se  : 

Subroutine  LPVB  updates  the  state  of  the  control 
charge  combustion  chamber  and  prepares  data  to  permit 
the  mass  transfer  to  the  chamber  of  the  gua  to  be 
represented  as  a  source  term. 

Calls  : 

ERTOP,  FORMCR,  FRATE. 

Called  by: 

INTAL,  NOV  SOB. 

VA 

LDMPS 

Purpo se  : 

Subroutiim  LUMPS  provides  a  trial  update  of  the  state 
of  all  lumped  parameter  regions.  It  also  computes 
derivatives  of  the  state  variables  with  respect  to 
the  mass  fluxes  to  the  lumped  parameter  region. 

These  are  used  by  B CAL  to  enforce  the  physical 
boundary  conditions. 

Calls : 

CEE  MR,  COVC.  DPDER,  DPDRE,  ERTOP,  G  A.  MMOL,  JCON, 

RCASE. 

Call ed  by : 

BCAL,  TTBUF. 

MESH 

Purpo  se  : 

Subroutine  MESH  determines  the  level  of  modeling  of 
each  region  and  assigns  mesh  points  to  the  continuum 
regions  on  the  basis  of  their  relative  sizes. 

Calls: 

FITSP,  OUTPUT,  REG12,  REG 21 ,  RBG32,  STATES. 

Called  by: 

INTEG,  NOV  SUB. 

NOV  SOD 

Purpo  se : 

Subroutine  NOVSUB  reads  and  prints  the  input  data 
used  to  define  the  problem.  NOVSUB  also  performs  all 
the  initialization  of  variables  vhich  is  not  accom¬ 
plished  through  BLKDAT.  Some  of  the  initialization 
is  performed  by  executing  a  call  to  certain 
subroutine  s  with  an  appropriately  set  switch. 
Following  initialization,  NOVSUB  transfers  program 
control  to  INTEG. 

Calls: 

AVN,  CUBUST,  EPTOR,  ERTOP,  FBAK,  FILLBR,  FILLER, 
FITSP,  FORM,  HEATW,  I  WED,  INTEG,  INTRPC,  INTRP1 , 
INTRP2 ,  LPDATA,  LFVB,  1ESB,  SETSUB,  TABLES.  TCDATA, 
TCINIT,  TaiNK, 

Called  by: 

XKTC, 

OUTPUT 

Purpo  se : 

Subroutine  OUTHJT  is  responsible  for  logout  of  the 
solution  to  the  printer  and/or  direct  access  device 
unit  8. 

Calls: 

TABLES. 

Called  by: 

BCAL,  FILLBR,  FILLER,  INTEG,  WSH, 
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PER*’ 19 

Pur  po  se  : 

Subroutine  PERF19  computes  the  surface  area  and 
covolume  of  nineteen-perforation  propellant  daring 
the  slivering  phase. 

Call  s : 

None . 

Call ed  by : 

FORM. 

PRTOC 

Purpose : 

Function  PRTOC  computes  the  square  of  the  speed  of 
sound  according  to  the  covoluae  equation  of  state. 

Calls: 

None. 

Called  by: 

BCAL,  INTAL,  IN  TEG,  TDBCF. 

PR  TOE 

Purpose  : 

Function  FRTOE  computes  the  internal  energy  from  th< 
pressure  and  density  according  to  the  covolume 
equation  of  state. 

Calls: 

None. 

Called  by: 

BCAL,  INTAL,  STATES,  T  CL  INK,  TDBCAL,  TDBUF. 

RCASE 

Pur  po  se : 

Subroutine  RCASE  determines  the  density  of  the  case 
as  a  function  of  pressure. 

Calls: 

None . 

Called  by: 

BDMOV,  INTAL,  INTRPC,  LUMPS. 

RDEDPR 

Purpose  : 

Function  RDEDPR  computes  the  quantity  p(de/dp)^ 
according  to  the  covolume  equation  of  state. 

Calls: 

None. 

Called  by: 

INTAL. 

REGRES 

Pur  po  se : 

Subroutine  RBGRES  computes  the  rate  of  regression  o 
the  surface  of  the  propellant. 

Calls: 

None . 

Ca  1 1  e  d  by  : 

INTEG. 
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RBG12 


REG  21 


REG  3  2 


RES  FUN 


RBEO 


ROOT 


Pnrpo  se  : 

Subroutine  REG12  transforms  data  for  a  region  of 
continuum  ullage  into  data  for  a  region  of  lumped 
parameter  ullage. 

Calls: 

COVC,  ERTOP,  INTO  PI,  INTRP2 . 

Called  by: 

MESH. 

Purpo  se  : 

Subroutine  REG21  transforms  data  for  a  region  of 
lumped  parameter  ullage  into  data  for  a  region  of 
continuum  ullage. 

Call s : 

INTRPl. 

Called  by: 

BRNOCT,  MESH. 

Purpo  se : 

Subroutine  REG32  sets  initial  data  for  a  lumped 
parameter  region  of  ullage  which  is  newly  opened. 

Calls : 

COVC,  ERTOP,  1NTRP1,  INTRP2 . 

Called  by: 

MESH. 

Purpose  : 

Function  RESFUN  computes  the  resistance  to  project 
motion  due  to  interference  of  the  rotating  band  wi 
the  gun  tube. 

Calls: 

None. 

Called  by: 

BDMOV. 

Pur  po  se  : 

Subroutine  RBEO  computes  the  rate  of  propagation  o 
intergranular  disturbances. 

Calls: 

None . 

Called  by: 

BDMOV,  INTAL. 

Purpo se : 

Subroutine  ROOT  selects  a  new  value  for  the  soluti 
of  an  arbitrary  equation  by  the  aiethod  of  regula 
falsi. 

Call s: 

None . 

Called  by: 

I M3  ED. 
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ROOTR 

Purpose : 

Subroutine  ROOIR  is  similar  to  ROOT,  but  contains 
additional  logic  for  the  situation  in  which  the 
islope  of  the  test  function  vanishes. 

Calls: 

None . 

Called  by: 

IMBED. 

SETSOB 

Purpose  : 

Subroutine  SETSUB  sets  the  reactivity  and  perme¬ 
ability  pointers. 

Calls: 

None. 

Called  by; 

F  ITS  P,  NOV  SOB, 

SLIVER 

Purpose  : 

Subroutine  SLIVER  computes  t>e  volume  and  surface 
area  of  nineteen-perforation  hexagonal  grains  during 
the  slivering  phase. 

Calls: 

None. 

Called  by: 

HEX19 . 

STATES 

Purpose  : 

Subroutine  STATES  computes  certain  of  the  dependent 
state  variables  from  the  principal  storage  arrays  for 
the  gas-phase. 

Calls: 

COVC,  u  A  MMOL,  PR  TOE. 

Called  by: 

IN  TEG,  MESH,  TERMIN. 

TABLES 

Pur  po  se  : 

Subroutine  TABLES  compiles  the  summary  data  which 
describe  the  conventional  interior  ballistic 
quantities  such  as  histories  of  breech  and  base 
pressure  and  of  the  projectile  trajectory. 

Ca  1 1  s : 

DINT,  INTRP1,  TaiNK. 

Called  by: 

BCAL,  FILLS  R,  FILLER.  IN  TEG,  NOV SOB,  OUTPUT. 
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TCARAC 

Purpo  se  : 

Subroutine  TCARAC  defines  characteristic 
data  in  conj  unction  with  the  determination  of 
boundary  values  for  the  traveling  charge 
when  modeled  as  a  continuum. 

Call s; 

TQ3R2 ,  T CDS  DR. 

Called  by: 

TCBASE,  TCBPR,  TCCNTC. 

TCBASE 

Purpo se : 

Subroutine  TCBASE  determines  the  regression 
rate  of  the  traveling  charge  and  provides 
trial  boundary  values  for  the  unreacted  side. 

Calls: 

TCARAC,  TCBOPN,  TCDSDR,  TCIG,  T CL  INK,  TCPROP, 
TCRESP,  TCRNOM. 

Called  by: 

TCXC. 

TCBMDC 

Purpo  se  : 

Subroutine  TCBMDC  computes  mechanical 
properties  of  the  traveling  charge  when  its 
rheology  is  specified  in  tabular  format. 

Calls: 

ISC1C0  (An  IMSL  routine). 

Called  by: 

TCDATA,  TCDSDR,  TCRNOM,  TCSNOW. 

TCBPR 

Purpose  : 

Subroutine  TCBPR  npdates  the  boundary  valnes 
of  the  vxaveling  charge  at  the  base  of  the 
proj  e  ctil  e 

Calls ; 

TCARAC,  TCDSDR,  TCRNOM. 

Called  by: 

TCXC. 

TCBR1 

Pur  po  se  : 

Subroutine  TC8R1  initializes  the  computational 
arrays  associated  wi  th  the  traveling  charge. 

Calls : 

TCRFIT,  TCRNOM. 

Call ed  by : 

TCINIT. 
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TGBR2 


Purpose:  Subroutine  TCSR2  moves  the  traveling  charge 


computational  arrays  into  working  arrays.  This 
routine  performs  a  trivial  function  in  the 
present  code  but  it  is  retained  in  order  to 
maintain  structural  cosmonality  with  BRLTC  where 
its  function  is  non- trivial. 

Calls: 

None . 

Called  by: 

TCARAC,  T(BR3,  TCBR4,  TCOUT. 

TCBR3 

Purpose  : 

Subroutine  TCBR3  updates  the  state  variables  for 
the  traveling  charge  at  all  internal  mesh  points. 

Calls: 

TCBR2,  TCDSDR,  TCWFR. 

Called  by: 

TCXC. 

TCBR4 

Purpose  : 

Subroutine  TCBR4  determines  the  CFL  time  step 
restriction  for  the  traveling  charge. 

Calls: 

TCBR2 ,  TCDSDR,  TO.  INK. 

Called  by: 

IN  TEG. 

TCBURN 

Purpo  se : 

Subroutine  TCBURN  calculates  the  regression  rate 
of  the  traveling  charge  when  specified  as  an 
empirical  function  of  pressure  (on  the  reacted 
side)  or  stress  (on  the  unxeacted  side). 

Calls: 

TCDVDI . 

Called  by: 

TCBASE, 

TCDATA 

Purpose : 

Subroutine  TCDATA  reads  and  prints  the  input  data 
associated  w ith  the  traveling  charge. 

Calls: 

TCBMDC,  TaNIT. 

Called  by: 

NOV  SUB. 

TCDSDR 

Purpose  : 

Function  TCDSDR  calculate*  (he  speo  1  of  sound  in 
the  traveling  charge. 

Calls: 

TCBMDC,  TCSNOM. 

Called  by: 

TCARAC,  TCBASE,  TCBPR,  TCBR3 ,  TCBR4,  TCCNTC. 
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TCDVD1 

Purpo  se : 

Subroutine  TCDVDI  performs  interpolation  by  the 
method  of  divided  differences. 

Call s : 

None . 

Called  by: 

TCBURN. 

TOGETK 

Purpo  se  : 

Subroutine  TOGETK  updates  the  ordinary 
differential  equations  associated  with  the 
traveling  charge. 

Calls : 

TCRESP. 

Called  by: 

TCXC. 

TCIG 

Purpose  : 

Subroutine  TCIG  determines  the  time  at  which 
ignition  occurs  for  the  traveling  charge. 

Calls  : 

None. 

Called  by: 

TCBASE. 

TCINIT 

Pur  po  se : 

Subroutine  TCINIT  performs  the  initialization  of 
constants  and  pointers  for  the  traveling  charge. 

Calls: 

TCBR1,  TCRFIT,  TCVUPD. 

Called  by: 

NOV  SOB.  TCDATA. 

TCLINK 

Pur  po  se  : 

Subroutine  TCLINK  transfers  data  between  the  NOVA 
routines  and  the  BRLTC  routines. 

Calls: 

COVC,  GAMMOL,  F&TOE,  TCFROP. 

Called  by: 

INTAL,  INTEG,  NOV  SUB,  TABLES,  TCBASE.  TCBR4 , 

TCOUT,  TCXC. 

TCONTC 

Purpose  : 

Subroutine  TCONTC  updates  the  boundary  values  for 
the  traveling  charge  at  the  interfaces  between 
increment  s. 

Calls : 

TCARAC,  TCDSDR,  TCRNOM. 

Called  by • 

TCXC. 
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TCOUT 

Purpo  se : 

Subroutine  TOOUT  performs  the  logout  of  the 
traveling  charge  state  variables. 

Calls: 

TCBR2,  TCLINK. 

Called  by: 

INTEG. 

TCPROP 

Purpose  : 

Subroutine  TCPROP  moves  the  traveling  charge 
property  data  from  vector  into  scalar  storage. 

Calls: 

None . 

Called  by: 

TCBASE,  TCLINK. 

TCRESP 

Purpo  se  : 

Subroutine  TCRESP  computes  the  resistance  on 
the  projectile  when  the  traveling  charge  option 
is  in  effect. 

Calls : 

None. 

Called  by: 

TCBASE.  TCGETK. 

TCRFIT 

Pur  po  se  : 

Subroutine  TCRFIT  performs  the  mesh  allocation 
for  the  traveling  charge  and  the  interpolation  i 
data  associated  with  revisions  to  the  number  of 
points. 

Calls: 

None. 

Call ed  by : 

TCBR1,  TaNIT,  TCXC. 

TCRNOM 

Purpo  se  : 

Function  TCRNOM  computes  the  density  of  the 
traveling  charge  as  a  function  of  pressure  on  tl 
nominal  loading  curve. 

Calls: 

TC8MDC . 

Called  by: 

TCBASE,  TCBPR,  TCBR1,  TCCNTC. 

TCSCHK 

Purpo  se  ; 

Subroutine  TCSCHK  enforces  the  requirement  that 
the  pressure  in  the  traveling  charge  not  exceed 
the  nominal  loading  value  for  the  current  value 
of  density. 

Calls: 

TCSNOM. 

Called  by: 

TCXC. 
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TCSNOM 

Pur po  se  : 

Function  TCSNOM  coaputes  the  pressure  in  the 
traveling  charge  as  a  function  of  density  on 
the  nominal  loading  curve. 

Calls : 

TCBMDC. 

Called  by : 

T CDS  DR,  TCSCHK. 

TCVCHK 

Pur  po  se : 

Subroutine  TCVCHK  checks  that  the  velocity  in 
the  traveling  charge  has  not  been  reversed  as  a 
consequence  of  the  friction  tera. 

Calls: 

None . 

Called  by: 

TCXC. 

TCVIIPD 

Purpose  : 

Subroutine  TCVUPD  updates  the  continuum  boundary 
velocities  for  the  traveling  charge. 

Calls: 

None . 

Called  by: 

TCIN  IT,  TCXC. 

TCTFK 

Pur  po  se  : 

Subroutine  TCTFR  computes  the  friction  between 
the  traveling  charge  and  the  tube  wall. 

Calls: 

None . 

Called  by: 

TCBR3. 

TCXC 

Purpo  se  : 

Subroutine  TCXC  is  the  integration  executive  for 
the  traveling  charge  state  variables. 

Calls: 

T CHASE,  TCBPR,  TC8R3 ,  TCGETK,  TCLINK,  TCONTC. 
TCRFIT,  TCSCHK,  TCVCHK,  TCVUPD,  TCZDPD. 

Called  by: 

INTAL. 

TCZDPD 

Purpose  : 

Subroutine  TCZDPD  updates  the  positions  of  the 
continuum  boundaries  associated  with  the 
traveling  charge. 

Calls: 

None . 

Called  by: 

TCXC. 
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TDBCAL 

Purpo  se : 

Subroutine  TDBCAL  uses  the  physical  boundary 
conditions  and  the  characteristic  foons  to  determine 
the  boundary  values  at  each  time  step. 

Calls: 

CCVC,  EPTOR,  FUJI,  GAMMOL.  PRTCE. 

Called  by: 

TDBUF. 

TI»DF 

Purpo  se : 

Subroutine  TDBUF  acts  as  an  interface  from  the  NOVA 
storage  conventions  to  those  of  the  solver  TDBCAL. 

Calls: 

CAL  FLO,  CAL  PR  M,  COVC,  FLUXDR,  GAMMOL,  LUMPS.  PR  TOC, 
PRTOE,  TDBCAL. 

Called  by: 

BCAL. 

TER  MIN 

Purpo  se  : 

Subroutine  TERMIN  tests  for  the  termination  of  the 
calculation. 

Calls: 

STATES. 

Called  by: 

INTEG. 

VIS 

Purpose  : 

Function  VIS  computes  the  viscosity  of  the  gas-phase 
as  a  function  of  temperature. 

Calls : 

None. 

Called  by: 

FRICT.  HEATP,  HEATW,  INTRPC. 

VOIDS 

Purpo  se  : 

Subroutine  VOIDS  updates  the  volume  fractions  of  ea 
of  the  constituent  solid-phase  species  from  the 
history  of  the  porosity. 

Calls: 

None . 

Called  by: 

INTEG. 

WALTEM 

Purpo  se : 

Subroutine  WALTEM  updates  the  thermal  equation  for 
the  tube  vail  temperature  determined  according  to  a 
cubic  profile  approximation. 

Call s: 

None . 

Called  by: 

INTEG. 
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TABLE  A.  2  SUMMARY  OF  XETC  "MANDATORY"  IN  POT  FILES 


FILE  NUM3ER 

FILE  NAME 

Ml 

Problem  Name 

M2 

Control  Data 

M3 

Integration  Parameters 

M4 

File  Counters 

M5 

Properties  of  Ambient  Gas 

M6 

General  Properties  of  Propellant  Bed 

M7 

Ith  Propellant  Type,  Location  and  Mas* 

M8 

Ith  Propellant  Form  Function  Data 

M9 

Ith  Propellant  Rheology 

M10 

Ith  Propellant  Solid-Phase  Thermochemistry 

Mil 

Ith  Propellant  Gas-Phase  Thermochemistry 

M12 

Tube  Geometry 

*M13 

Projectile  and  Rifling  Characteristics 

*  File  M13 

is  not  mandatory  in  XKTC.  It  is  required  only  if 

traveling 

charge  option  is  not  exercised. 
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TABLE  A. 3  SUMMARY  OF  XKTC  CONTINGENT  INPUT  FILES 


FILE  NUMBER 

Cl 

C2 

C3 

C4 

C5 


C5.1 
C5  .2 
C6 
C7 
C8 
C9 
CIO 
Cll 
C12 
C13 
C14 
C15 
C16 
C17 
C18 
C19 
C19.1 
C19.2 
C20 
C21 
C22 


FILE  NAME 

Logout  Counters 

Erosive  Burn  Rate  Parameters 

Not  Supported. 

Position  and  Mass  of  Additional  Bags  of 
Ith  Propellant 

Distribution  of  Volume  Fraction  of 
Ith  Propellant 

Ith  Propellant  Port  Diameter  File  Counter 
Ith  Propellant  Port  Diameter 
Not  Supported. 

Ith  Propellant  Deterred  Layer  Properties 

Igniter  Themochemi stry 

Igniter  Discharge  Times 

Igniter  Discharge  Positions 

Igniter  Rate  of  Discharge 

Bore  Resistance 

Tube  Thermal  Properties 

Tube  Initial  Temperature  Profile 

Not  Supported. 

Not  Supported. 

Not  Supported. 

Ith  Forward  Compactible  Filler  Element  Properties 
Constitutive  Data  for  Ith  Forward  Element 
Ith  Rear  Compactible  Filler  Element  Properties 
Constitutive  Data  for  Ith  Rear  Element 
Not  Supported. 

Not  Supported. 

Not  Supported. 
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TABLE  A.3  continued 


FILE  NOMBTt 


FILE  NAME 


C23 

C23 . 01 
C23  .1 
C23.2 
C23  .3 

C23  .4 
C23.5 
C23  .6 
C23  .61 

C23  .62 
C23  .63 
C23.7 

C23  .71 

C23  .8 
C23  .9 
C24 

C25 

■’.5.1 

25  .2 
C25  .3 
Ci 5  .31 
C25.4 
C25  .5 
C25  .6 


Positions  for  Pressure  Table  Storage 
Projectile  Afterbody  Pressure  History  Option 
Kinetics  Option  Counters 
Properties  of  Species 

Composition  of  Propellant  Near  Field  Combustion 
Products 

Composition  of  Igniter  Discharge  Combustion  Products 

Composition  of  Initial  Ambient  Gas 

Initial  Composition  of  Solid  Propellant 

Composition  of  Traveling  Charge  Near  Field 
Combustion  Products 

Types  of  Reactions  in  Miztnre  of  Combustion  Products 

Species  Thermal  Equilibration  Switches 

Description  of  Arrhenius  Reactions  in  Mixture  of 
Combustion  Products 

Description  of  Pressure-Dependent  Reactions  in 
Mixture  of  Combustion  Product s 

Description  of  Reactions  in  Solid  Propellant 

Solid  Propellant  Thermal  Response  Parameters 

Parameters  to  Define  Mesh  in  Invariant  Embedding 
Analy  si s 

Tank  Gun  Option  Control  Data 
Geometry  of  Projectile  Afterbody 

Thickness  and  Density  Counters  for  Reactive  Layer  I 
Thickness  of  Reactive  Layer  I 
Segment  Pointers  for  Reactive  Layer  I 
Density  of  Reactive  Layer  I 

Counter  to  Describe  Discharge  of  Reactive  Layer  I 
Positions  to  Describe  Discharge  of  Reactive  Layer  I 
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TABLE  A. 3  continued 


FILE  NDffiER 


FILE  NAME 


C25.7 
C25.8 
C25.9 
C2 6 
C26  .1 
CIS  .2 
C26.21 


C26.3 
C26.4 
C26.5 
C26  .6 
C26  .7 
C26.71 

C26.8 

C26.9 

C27 

C28 

C28.1 

C28.2 

C28.3 

C28.4 

C28.3 

C28.6 

C28.7 


Tines  to  Describe  Discharge  of  Reactive  Layer  I 

Rate  of  Discharge  of  Reactive  Layer  I 

Bure  Rate  Counters  for  Reactive  Layer  I 

Burn  Rate  Data  for  Reactive  Layer  I 

Ignition  Data  for  Reactive  Layer  I 

The  mo  chemical  Data  for  Reactive  Layer  I 

Composition  of  Reactive  Layer  Near  Field  Combustion 
Product  s 

Properties  of  Deterred  Layer  of  Reactive  Layer  1 

Enthrall  Property  Pointers 

Permeability  Model  Data 

File  Counters  for  1-th  Reactivity  Model 

Thermochemical  Data  for  I-th  Reactivity  Model 

Composition  of  Near  Field  Combustion  Products  of 
I-th  Reactivity  Model 

Ignition  Value  for  I-th  Reactivity  Model 

Burn  Rate  Data  for  I-th  Reactivity  Model 

Tabular  Description  of  I-th  Reactivity  Model 

Internal  Properties  of  Control  Charge  Combustion 
Chamber 

External  Properties  of  Control  Charge  Combustion 
Chamber 

Control  Charge  Type 
Control  Charge  Form  Function 
Control  Charge  Burn  Rate  Counter 
Control  Charge  Burn  Rate  Description 
Control  Charge  Thermochemistry 
Properties  of  Deterred  Layer 
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TABLE  A. 4  SUMMARY  OF  XKTC  TRAVELING  CHARGE  INPUT  FILES 


File  Noaber 


File  Naae 


TCI 

TC2 

TC3 

TC4 

TC5 

TC6 

TC7 

TC8 

TC9 

TC10 

TC11 

TC12 

TC13 

TC14 

TC15 

TC16 

TC17 


Traveling  Charge  Control  Data 
Mesh  Parameters 
General  Properties 

Propellant  Ihemocheaical  Properties 

Burn  Rate  Switches 

Tabular  Burn  Rate  Data 

Tabnler  Born  Rate  Data  (cont'd) 

Exponential  Born  Rate  Data 

Propellant  Analytical  Rheology  Data 

Propellant  Tabular  Rheology  Data 

Propellant  Tabular  Rheology  Data  (cont'd) 

Ablative  Fils  Data 

Propellant  Friction  Data 

Tabular  Bore  Resistance  Data 

Barrel  Shock  Resistance  Data 

Obturator  Setback  Resistance  Data 

Obturator  Friction  Data 
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TABLE  A. 5  DESCRIPTION  OF  INPUT  FILES 


File  [Ml] 

"Problem  Name" 

( 15  AO  One 

Card 

TITLE 

Problem  name. 

May  be  up  to 

£0  alphanumeric  characters. 

File  [M2] 

"Control  Data* 

(711 ,3311) 

One  Card 

NPRINT  1  if  print  on  logout,  0  otherwise. 

NGRAJPB  Inactive. 

NDISK  1  if  disc  write  on  logout,  0  otherwise. 

DSKRD  1  if  disc  start,  0  otherwise. 

1BTABL  1  if  summarized  interior  ballistic  data  are  required, 

0  otherwise. 

NFL  AM  1  if  summary  of  ignition  delay  data  required, 

0  otherwise. 

NPTABL  1  if  summarized  pressure,  pressure  difference  data  are 

required, 

0  otherwise, 

NEROS  1  if  erosive  effect  to  be  included  in  propellant  combustion, 

0  otherwise.  If  NEROS  *  1,  File  [C2]  is  required. 

NDTN  Inactive. 

MITW  0  if  wall  temperature  is  not  updated. 

1  if  cubic  profile  update  is  used. 

NBC(l)  Inactive. 

N}C(2)  Inactive. 

NRES(l)  Inactive. 

NRES(2)  Inactive. 

LDBED  0  if  propellant  beds  initially  uncompacted, 

1  otherwise.  See  connent  in  File  [M91  . 
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JHTW 


LYER 

1BRES 


NTC 


IfBIB(l) 

ir«IB(2) 
INBIB(  3 ) 
NXCW 
NBLDWN 
NSWSOL 

KMODE 


0  if  heat  lots  to  tube  is  neglected. 

1  if  heat  loss  is  calculated  by  empirical  correlation  for 
fully  developed  tnrbnlent  flow.  File  [C13]  is  required. 

Inactive. 

1  if  linear  interpolation  of  tabular  data  from  File  [C12]  is 
used  to  define  bore  resistance.  No  velocity  dependence  is 
considered. 

2  if  the  interpolated  value  is  multiplied  by  7.2V~9**,  where 
V  is  the  projectile  velocity  (ft/sec),  provided  V  1  27  fps. 

3  if  the  interpolated  value  is  multiplied  by 

(1  +  0.0004414V) /(I  +  0.005046V)  where  V  is  the  projectile 
velocity  (in/sec). 

If  any  other  value  is  used,  IBRES  will  default  to  the  value  2. 

0  if  traveling  charge  option  not  required. 

1  if  traveling  charge  option  is  in  effect. 

When  NTC  «  1,  NONIN  and  NUN00T  (File  [M3])  must  not  be  entered 
as  nine  (default  values  are  eight)  $  NEL  (File  [M4] )  must  be 
entered  as  zeroi  NDIM  (File  [M3])  must  not  exceed  ninety-eight 
for  a  granular  booster  or  forty-nine  for  a  stick  booster. 

Files  [TCI]  through  [TC17]  are  required. 

0  if  propellant  type  1  does  not  have  a  deterred  layer. 

1  if  propellant  type  1  does  have  a  deterred  layer. 

File  [C7]  is  required  in  this  case. 

As  per  INRIB(l),  but  for  second  propellant  type. 

As  per  INBIB(l),  but  for  third  propellant  type. 

Inactive. 

Inactive. 

0  if  non-conservative  form  of  solid-phase  continuity  equation 
is  to  be  used. 

1  if  conservative  fosa  is  to  be  used.  It  is  recommended  that 
this  option  be  selected  only  if  excessive  mass  defects  result 
from  the  use  of  the  non-conservative  form. 

0  if  kinetics  option  not  desired, 

1  if  gas-phase  is  to  be  treated  as  a  reacting,  homogeneous 
mixture.  See  Files  [C23.1]  through  [C23.9]. 

2  if,  in  addition  to  treatment  of  gas-phase  as  a  mixture, 
subsurface  reactivity  and/or  finite  difference  solution  of 
solid-phase  thermal  response  is  desired.  See  Files  [C23.1] 
through  [C23.9]. 
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MDDET  0  if  tank  gun  option  not  desired. 

1  if  tank  gun  option  is  desired.  See  Files  [C25.1]  through 
[C27] . 

NECHO  0  if  interpretation  of  input  and  ooml  run  desired. 

1  if  input  data  not  interpreted  but  normal  run  desired. 

2  if  interpretation  of  input  data  with  no  subsequent  run 
desired. 

Note:  Independently  of  the  value  of  NECBO,  the  Code  prints  an  image  of  the 
card  input  file. 


File  [M3]  "Integration  Parameters”  (6I5,F10.0/8F10 .0)  Two  Cards 


NDIM  Number  of  grid  points  to  be  used  in  calculation.  In  general, 

NDIM  must  not  exceed  99.  If  perforated  stick  propellant  is 
present,  NDIM  must  not  exceed  49. 

NS1EP  Number  of  integration  steps  before  logout.  (Each  step 

consists  of  a  predictor  followed  by  a  corrector.)  If 
NSTEP  =  0,  logout  follows  every  predictor  as  well  as  every 
corrector.  If  NSTEP  <  0,  File  [Cl]  is  required  with  iNSTEPl 
pairs  of  data.  NSTEP  1  -10. 

NDTST  Step  number  for  disc  start. 

NSTOP*  Number  of  steps  for  termination. 

NONIN  Unit  number  for  disc  read.  If  NDNIN  *  0,  it  is  defaulted 

to  8. 

N0N00T  Onit  number  for  disc  logout.  If  NUN00T  =  0,  it  is  defaulted 

to  8. 

DTPRT  Time  interval  before  logout  (sec).  If  entered  as  zero,  it  is 

defaulted  to  10lo» 

TSTOP*  Time  for  termination  (sec). 

ZSTOP*  Projectile  travel  for  termination  (in). 

TINT  Maximum  time  step  (sec). 

SAFE  Stability  safety  factor.  Recommended  value  is  1.1  although 

larger  values  may  sometimes  be  required. 

CRIT  Stability  safety  factor  for  source  terms.  Recommended  value 

is  0.05. 
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RZOLV 


Spatial  resolution  factor.  Recomoended  range  of  values  is 
0.01  i  RZOLV  i  0.05.  If  RZOLV  is  entered  as  zero,  it  will 
default  internally  to  1/(NDIM  -  1). 

TABLIB**  Time  interval  for  storage  of  summarized  interior  ballistic 

data. 

TABLP**  Time  interval  for  pressure  table  storage. 


Termination  occurs  at  whichever  of  NSTOP,  TSTOP  or  ZSTOP  is  first 
satisfied, 

Table  sizes  are  dimensioned  to  100.  When  overflow  is  imminent,  tables 
are  condensed  by  deleting  every  second  datum.  Subsequently,  storage 
time  is  doubled. 


File  [Cl]  "Logout  Counters*  (1615)  One  or  TVro  Cards 


N.B.  File  required  if  and  only  if  NSTEP  <  0.  (See  File  [M3].) 

NDTEP(l)  Maximum  value  of  integration  step  number  for  which 

MSTEP(l)  is  to  be  used. 

MSTEP(l)  Number  of  integration  steps  (predictor  plus  corrector 

count  as  one  unit)  between  logout  cycles. 


NDTE?(  iNSTEPi )*  Maximum  value  of  integration  step  number  for  which 
MSTEP (INSTEP I )  is  to  be  used. 

MSTEP( I NSTEP I )  Corresponding  value  of  number  of  steps  between  logout 
cycl es. 


If  NHTEP( INSTEP] )  is  exceeded  during  the  calculation,  the  value 
MSTEP (INSTEPi)  will  be  used  as  a  default  quantity. 
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File  [M4] 


"File  Counters 


(1615)  Two  Cards 


NSTA 

JJ 

II 

NBRES 

NTEM 

NEL 

NPROP 

NBRDS 

NEPS(l) 

NEPS(2) 
NEPS( 3) 
NZFT 

NHS(l) 


Number  of  entries  in  tube  geometry  table.  2  £  NSTA  £  10. 

Number  of  times  at  which  primer  discharge  is  specified. 

0  <  JJ  <  40. 

If  II  and  JJ  1  2,  Files  [C8]  -  [Cll]  are  required. 

Number  of  positions  at  which  primer  discharge  is  specified. 
0  <  II  <  40. 

Number  of  entries  in  bore  resistance  table. 

0  £  Mi  RES  <  20. 

If  NBRES  2  2,  File  [C12]  is  required. 

Number  of  entries  in  tube  initial  temperature  profile. 

If  NTEM  =  0,  tube  temperature  defaults  to  TEMST  (see 
File  [M5] ) .  If  NTEM  t  0,  File  [C14]  is  required. 

0  <  NTEM  £  10. 

Number  of  elements  to  characterize  compactible  filler 
between  bed  and  projectile.  If  NEL  =  0,  projectile  base 
is  taken  as  right-hand  boundary  for  computational  domain. 

If  NEL  it  0,  File  [C18]  is  required.  0  <  NEL  <  10. 

Number  of  types  of  propellant  grains.  1  £  NPROP  £  3. 

Number  of  burn  rate  data  sets  to  describe  pressure 
dependence  of  exponent  and  pre-exponential  factors. 

See  File  [M10]  .  Note  that  the  same  value  is  assumed  for 
all  propellant  species,  if  more  than  one  are  defined. 

Number  of  entries  in  initial  distribution  of  volume 
fraction  of  propellant  type  1.  If  NEPS  =  0,  porosity  is 
calculated  from  mass  and  position  data  given  in  File  [M7 3 . 
If  NEPS  /=  0,  File  IC5]  is  required.  Propellant  types  may 
be  entered  in  either  mode,  independently  of  the 
representation  of  other  types.  0  £  NEPS(  1)  £  10. 

Number  of  entries  for  type  2. 

Number  of  entries  for  type  3. 

Number  of  entries  in  table  of  positions  for  summarized 
pressure  data.  If  hEFT  t  0,  File  [C23]  is  required. 

0  £  NZPT  <  8. 

Inactive. 
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NMS(  2) 

Inactive. 

Wf  IB 

Inactive. 

M0RE(1) 

Nuaber  of  additional  bags  of  propellant  of  type  1. 

If  MORE(l)  >  0,  File  [C4]  ia  required.  0  <  M0RE(1)  <  9. 

MORE ( 2) 

Nuaber  of  additional  baga  of  propellant  of  type  2.  Only 
required  if  NlfcOP  1  2.  Starts  a  new  card. 

MORE ( 3 ) 

Nuaber  of  additional  bags  of  propellant  of  type  3.  Only 
required  if  NPROP  -  3. 

NBL 

Number  of  eleaents  to  characterize  oonpactible  filler 
between  bed  and  breech.  If  MIL  i  0,  File  [C19.1]  is 
required.  0  i  NBL  1  10. 

File  [MS] 

"Properties  of  Ambient  Gas"  (8F10.0)  One  Card 

TRUST 

Initial  temperature  (R)  . 

PST 

Initial  pressure  (psi). 

GMST 

Molecular  weight  ( lbm/ lbnol ) . 

GAMST 

Ratio  of  specific  heats  (-) . 

File  [M6] 

"General  Properties  of  Propellant  Bed"  (8F10.0)  One  Card 

TPO 

Initial  temperature  (R)  . 

File  [C2] 

"Erosive  Burn  Rate  Paraaeters*  (8F10.0)  Tie  Card 

N.B.  File  required  if  and  only  if  NEROS  =  1.  (See  File  [M2].) 

CEROS  Erosive  burning  pre-exponent ial  factor  (in*  #R/lbf). 

BEROS  Erosive  burning  exponential  factor  (-). 
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Note:  The  subsequent  files  [M7] ,  [C4] .  [C5] ,  [M8],  [M9] ,  [M10] ,  [C6] . 

[Mil]  end  [C7]  are  repeated,  as  a  group  for  each  of  the  NPR0P  (see 
(File  [M4] )  types  of  propellant  which  constitute  the  solid-phase. 


File  [M7]  "Ith  Propellant  Type,  Location  and  Mass" 

(5A4,6F10.0)  One  Card 


GRNAM(I) 
ZGR( 1, I) 

ZGR(2, I) 

MTCRA(1, I) 


RHOP(I) 

RGRI(l.I) 

RGRO(l.I) 


Name  of  propellant.  Up  to  20  alphanumeric  characters. 

Left-hand  boundary  of  first  bag  of  propellant  of 
type  I  (  in) . 

Right-hand  boundary  of  first  bag  of  propellant  of 
type  I  ( in) . 

Mass  of  first  bag  of  propellant  of  type  I  (lbm). 

A  negative  value  of  WIGRA  may  be  entered.  In  such  a  case, 
the  absolute  value  is  used  to  define  the  mass.  The 
negative  sign  is  used  to  set  an  internal  switch  which 
causes  the  computation  of  initial  porosity  to  reflect  the 
assumption  that  the  propellant  is  packaged  as  a  cylinder 
rather  than  filling  uniformly  each  cross-section  of  the 
tube . 

Density  of  propellant  type  I  ( lbm/ in’ ) . 

Inner  radius  of  rear  of  bag  (in).  Only  required  if 
MODET  =  1.  (See  File  [M4].)  If  RGRI(1,I)  is  input  as  zero 
it  is  assumed  to  coincide  with  the  inner  radius  of  the 
available  flow  cross-section. 

Outer  radius  of  rear  of  bag  (in).  Only  required  if 
MODET  =1.  (See  File  [M4].)  If  RGR0(1,I)  is  input  as  zero 
it  is  assumed  to  coincide  with  the  outer  radius  of  the 
available  flow  cross-section. 
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File  [C4] 


'Position  and  Mats  of  Additional  Bags  of  Ith  Propellant” 
(8F10.0)  MORE (I )  Cards 


N.B.  File  required  if  and  only  if  MORE(I)  =  0.  (See  File  [M4].) 

ZGR(3,I)  Left-hand  boundary  of  second  bag  of  propellant  o 

type  I  ( in) . 


ZGR(4, I) 

WH}RA(2, I ) 
BGRI(2, I) 
RGR0(2, 1) 
ZGR(5, I)* 


Right-hand  boundary  of  second  bag  of  propellant 
type  I  (in). 

Mass  of  second  bag  of  type  1  (lbm). 

Inner  radios  of  second  bag  (in). 

Outer  radius  of  second  bag  (in). 

Left-hand  boundary  of  third  bag. 


f 

of 


ZGR(2*MORE(I)+2, I)  Right-hand  boundary  of  bag  MORE(l)+l  (last  bag). 
WTCRA(MORE(I)+l)  Mass  of  last  bag. 

RGRI(MORE(I)+l)  Inner  radius  of  last  bag  (in). 

/ 

RGRO(MORE(I)+  1)  Outer  radius  of  last  bag  (in). 

*  Begin  new  card  for  each  bag. 


File  [C5]  "Distribution  of  Volume  Fraction  of  Ith  Propellant” 

(8F10.0)  One  to  Three  Cards 

N.B,  File  required  if  and  only  if  NEPS(I)  =  0.  (See  File  [M4] .) 
ZEPS(1,I)  First  position  relative  to  breech  (in). 

EPSO(l.I)  Corresponding  volume  fraction  of  propellant  type  I. 


ZEPS(NEPS(I) , I)  Last  position. 

EPSO(NEPSd) ,  I)  Corresponding  volume  fraction. 
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File  [M8]  " I th.  Propellant  Fora  Function  Data”  (13 , 5F10  .0, 15  ,F10 .0) 

One  Card 


NFORM(I)  Fora  fraction  indicator.  1  i  |NF0RN(1)|  <  16.  N FORM  ■ ay 

be  entered  ae  a  negative  number.  In  snob  a  cate,  the 
absolute  value  is  used  to  determine  tbe  fora  fraction,  but 
the  grains  are  tale.)  be  stacked.  Ibis  option  allows 
granular  propellant  J^o  respond  as  though  it  had  the  flow 
resistance  of  stick  propellant. 

(H)(1)  Grain  dimension  (in).  (See  table  below.) 

BPEB.F  ( I )  Grain  dimension  (in).  (See  table  below.) 

GLIft(I)  Grain  length  (in). 

NH&Fd)  Number  of  perforations  (-)  except  for  NFORM  =  14  in  which 

case  NPERF  is  the  number  of  slots. 

SLW(I)  Grain  dimension  (in).  (See  table  below.) 

NFIX(I)  If  0,  grains  are  assumed  to  be  free  to  move. 

xi  1,  grains  are  assmed  to  be  attached  to  tbe  tube. 

If  2,  grains  are  assmed  to  be  attached  to  the  projectile. 

Strength  of  bond  to  tube  or  projectile  (lbf).  If  BQNDl(I) 
is  entered  as  zero,  it  is  taken  to  be  infinite.  Separation 
from  tube  or  projectile  occurs  when  absolute  value  of  force 
required  to  maintain  attachment  exceeds  BGNDX(I). 


BCNDX(I) 


FUNCTION  PARAMETERS 


FORM 


NFOEM 

GRAIN  TYPE 

OP 

DPERF 

SLV 

1 

Sphere 

Diameter  (in) 

- 

- 

2 

Cylinder 

Diameter  (in) 

- 

- 

3 

Stick* 

Diameter  (in) 

Perforation 
Diameter  (in) 

— 

4 

Strip* 

Vidth  (in) 

Thickness  (in) 

- 

5 

Mo  nop  erf  orated 

Diameter  (in) 

Perforation 
Diameter  (in) 

“ 

6 

Nonoperf orated  with 
outside  inhibition 

Diameter  (in) 

Perforation 
Diameter  (in) 

— 

7 

Seven  Perforations 

Diameter  (in) 

Perforation 
Diameter  (in) 

— 

9 

Nineteen  Perforations 

Diameter  (in) 

Perforation 
Diameter  (in) 

— 

10 

Hexagonal  Seven- 
Perforation  Stick 

Distance  between 
flats  (in) 

Perforation 
Diameter  (in) 

— 

11 

Slotted  stick* 

( singl e-voidage ) 

Diameter  (in) 

Perforation 
Diameter  (in) 

Slot 

Vidth(in) 

12 

Scroll* 

(singl e-voida ge ) 

Vidth  (in) 

Thickness  (in) 

— 

13 

Scroll* 

( dual-voids ge ) 

Vidth  (in) 

Thickness  (in) 

External 
Diameter  of 
Soroll  (in) 

14 

Blind  slit 

Diameter  (in) 

Diameter  of 
slotted  core  (in) 

15 

Hexagonal  Nineteen- 
Perforation 

Diameter  (in)** 

Perforation 
Diameter  (in) 

— 

16 

Monolithic  Grain 
with  central  port*** 

Port  diameter 
at  rear  (in) 

Port  diameter 
at  front  (in) 

• 

End  burning  is  neglected 

for  these  forms. 

Hence  GLEN  (File 

[M8] ) 

may  be  entered  arbitrarily. 

Grain-  assumed  to  have  rounded  corners.  Diameter  is  between  opposite 
corners. 

Burning  is  assoeed  to  take  place  only  on  the  snrfaces  of  the  internal 
port.  Projectile  afterbody  say  intrude  into  port.  NFIX(I)  and 
BCNDX(I)  default  to  1  and  0  respectively.  If  OD(I)  *  0,  or 
TPERFd)  =  0,  then  Files  [CS.l],  [CS.2]  are  required  to  define  the 
4i  sSvtsr  cf  the  pc?£  •  •  •  f  sc  ties  of  po  si  ti on. 


09 


1 


File  [C5.1]  * Ith  Propellant  Port  Dimeter  File  Counter"  (15)  One  Card 

N.B.  File  required  if  and  only  if  NFCRM  (I)  =  16  end  either  OD(I)  -  0  or 
DPERF(I)  =  0.  (See  File  [M8] . ) 

MNOP(I)  Number  of  entries  in  file  to  define  dimeter  of  port  of 

monolithic  grain.  2  i  WOP  (I)  &  8. 


File  [CS.2]  "Ith  Propellent  Port  Dimeter"  (8F10.0)  One  to  TVo  Cards 

N.B.  File  required  if  and  only  if  NFCRM(I)  =  16  and  either  OD(I)  -  0  or 
EPERF(I)  =  0.  (See  File  [M8]  .) 

ZfrNO(l.I)  First  axial  position  relative  to  rear  surface  of 

grain  ( in)  . 

DMNO(l.I)  Corresponding  dimeter  of  port  (in). 


Z*NO(MNOP(I),I)  Last  axial  position. 
EMNO(MNOP(I) ,  I)  Corresponding  dimeter. 


File  [M9]  "Ith  Propellant  Rheology"  (8F10.0)  One  Card 


GAP(I)  Rate  of  propagation  of  intergranular  stress  in 

settled,  loading  granular  bed  of  type  I  (in/sec).  If 
propellant  I  consists  of  stick  then  GAP(I)  is  independent 
of  the  settling  porosity. 

GEPO(I)  Porosity  of  settled  bed  for  type  I  propellant.  If 

LDBED  =  0  (see  File  [M2])  and  the  porosity  of  the  bed  is 
less  than  GEPO(I)  at  any  point  where  type  I  is  present, 
GEPO(I)  is  automatically  defaulted  to  the  smallest  value  of 
porosity  which  occurs  in  the  initial  distribution.  This 
property  may  be  used  to  avoid  calculating  GEPO(I). 
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GCAP(I) 

Rate  of  propagation  of  intergrannl ar  stress  during 
unloading  and  reloading  of  propellant  type  I  (in/sec). 

Not  required  for  stick  propellant. 

PRPYLD 

Inactive. 

ANO(I) 

Poisson  ratio  of  propellant  (-).  Must  be  specified  for 
stick  propellant.  AKU(I)  defaults  to  0.5  internally  if 
propellant  I  is  granular  or  if  ANTJ(I)  is  entered  as  zero. 

File  [M10] 

"Ith  Propellant  Solid-Phase  Themocheai stry"  (8F10.0) 

One  to  Three  Cards 

DPPR(1,I) 

Maximum  value  of  pressure  for  corresponding  values  of  burn 
rate  pre-exponent ial  and  exponential  factors  (psi)  . 

B22(1,I) 

Burning  rate  pre- exponent ial  factor  ( in/aec-ps ' n) . 

BNN( 1, I) 

* 

Burning  rate  exponent  (-) . 

• 

e 

D  PPR  (  NBRDS,  I )  * 

Maximum  value  for  last  set  of  burn  rate  data. 

B22(  NBRDS,  1 ) 

Corresponding  pre-exponential  factor. 

BNN(  NBRDS,  I ) 

Corresponding  exponent. 

Bid) 

Burning  rate  additive  constant  (in/  sec  . 

lEMPIG(I) 

Ignition  temperature  (R) . 

KP(I) 

Thermal  conductivity  (lbf/sec/R). 

ALPBAP(I) 

Thermal  diffusivity  (in*/sec). 

ALPHA  ( I ) 

Emissivity  factor  (-). 

If  pressure  exceeds  OPPR(MJRDS,  I)  ,  the  corresponding  values  of 
pre-exponential  and  exponential  factors  are  used  as  defanlt  valnes. 


Ill 


File  [Mil]  *Ith  Propellant  Gas-Phase  Thexmocheai stry*  (8F10.0) 

One  Card 

ECHEM(I)  Internal  energy  released  in  combustion  (lbf-in/lbm)  . 

GMW(I)  Molecular  weight  (lbm/lbmol). 

GAMMA(I)  Ratio  of  specific  heat  (-). 

BV(I )  Covolume  (in*/lbm). 

File  [C7]  "Ith  Propellant  Deterred  Layer  Properties*  (8F10.0) 

One  Card 

N.B.  File  required  if  and  only  if  INEIB(I)  +  0.  (See  File  [M2].) 

ECBIB(l)*  Internet  energy  released  in  combustion  at  start  of 

deterred  layer  (lbf-in/lbm). 

ECHIB2(I)  Internal  energy  released  in  combustion  at  end  of 

deterred  layer  (lbf-in/lbm). 

RGFAC(I)*  Factor  by  which  burn  rate  is  multiplied  at  start  of 

de terred  1 ayer  (-) . 

RGFAC2 ( I )  Factor  by  which  burn  rate  is  multiplied  at  end  of 

deterred  layer  (-)  . 

HIBX(I)  Depth  of  inhibited  layer  (in). 

Values  within  deterred  layer  deduced  by  linear  spacewise 
interpolation.  Final  values  need  not  be  the  same  as  those  of  the 
undeterred  propellant. 

File  [C8]  "Igniter  Thermochemistry"  (8F10.0)  One  Card 

N.B.  File  required  if  and  only  if  II,  JJ  1  2.  (See  File  [M4]  .) 

EIG  Internal  energy  released  in  combustion  (lbf-in/lbm). 

GMIG  Molecular  weight  (lbm/lbmol). 

GAMIG  Ratio  of  specific  heats(-). 

VIG  Specific  volume  of  igniter  solid-phase  (in’/lbo). 
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File  [C9l  'Igniter  Di (charge  Tinea"  (8F10.0)  One  to  Five  Cerda 

N.B.  File  required  if  and  only  if  II,  JJ  1  2.  (See  File  [K4] .) 

TFBI(l)  Fir  at  tine  in  tabular  repreaentation  of  igniter  diacharge 

.  (sec). 

TPHI(JJ)  Last  tiae. 


File  [CIO]  'Igniter  Discharge  Positions*  (8F10.0)  One  to  Five  Cards 

N.B.  File  required  if  and  only  if  II,  JJ  >  2.  (See  File  [M4].) 

ZPHI(l)  First  position  in  tabnlar  representation  of  igniter 

.  discharge  (in). 

e 

ZPBI(II)  Last  position. 


File  [Cll]  'Igniter  Rate  of  Discharge*  (8P10.0)  Two  to  Two  Hundred 

Cards 


N.B.  File  required  if  and  only  if  II,  JJ  1  2 .  (See  File  [M4].) 

FBI( 1,1)  Igniter  discharge  per  nnit  length  per  unit  time  at  first 

position  and  first  tine  (lba/ in/tec) . 

PHI ( 2 , 1 )  Discharge  at  second  position  and  first  tine. 


PHI(II.l)  Discharge  at  last  position  and  first  tiae. 

FHI(1, 2)*  Discharge  at  first  position  and  second  tiae. 


rai(II.JJ)  Discharge  at  last  position  and  last  tiae. 


PHKl.J)  begins  a  new  card,  J  =  1 
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JJ. 


File  [M12] 


Tube  Geometry"  (8F10.0)  One  to  Three  Cards 


ZA(  1) 

First  axial  location  relative  to  breech  (ins). 

RA(  1) 

• 

* 

Corresponding  radios  of  bore 

(ins) . 

• 

* 

ZA(NSTA) 

Last  axial  location. 

RA(NSTA) 

Corresponding  radius. 

File  [C12] 

•Bore  Resistance*  (8F10.0) 

One  to  Five  Cards 

N.B  File 

required  if  and  only  if  NBRES  1  2. 

(See  File  [M4]  .) 

BRZ(l) 

First  axial  location.  If  BRZ(l)  >  0  it  is  understood  to 
be  relative  to  breech  (in).  If  BRZ(l)  is  entered  as  zero, 
however,  all  values  of  BIZ  are  understood  to  be  relative  to 
the  initial  position  of  the  projectile. 

BR(1) 

• 

* 

Corresponding  resistance  exerted  on  projectile  when  base  is 
at  BRZ(l)  (psi). 

• 

« 

BRZ(NSRES) 

Last  location. 

BR(N)RES) 

Corresponding  resistance. 

File  [C13] 

"Tube  Thermal  Properties" 

(8F10.0)  One  Card 

N.B.  File  required  if  and  only  if  JH1W  t  0.  (See  File  [M2].) 
KW  Thermal  conductivity  (lbf/sec/R). 

ALPHAS?  Thermal  diffuaivity  (in*/sec). 


Emissivity  factor  (-) . 


ALFAW 


File  [C 14] 


(8F10.0) 


"Tube  Initial  Temperature  Profile* 

One  to  Three  Cards 

N.B.  File  required  if  and  only  if  N1EM  jh  0.  (See  File  [M4]  .) 

ZW(1)  First  axial  location  relative  to  breech  (in). 

TEMW(l)  Corresponding  wall  temperature  (R) . 


e 

ZV(NTEM)  Last  location. 

lEMW(NTEM)  Corresponding  temperature. 


File  [M13]  "Projectile  and  Rifling  Characteristics"  (8F10.0) 

One  Card 


ZBPR  Initial  axial  location  of  ba^e  of  projectile  relative  to 

breech  (in) . 

WTPR  Mass  of  projectile  (lbm). 

PRIN  Polar  moment  of  inertia  of  projectile  (lbm-in1). 

RIF  Angle  of  rifling(°). 


File  [C18j  "Ith  Forward  Coopactible  Filler  Element  Properties* 

(3F10. 0,215)  One  Card 

N.B.  File  required  if  and  only  if  NEL  t  0.  (See  File  [M4] .)  File  is 
repeated  NEL  tunes. 

XEL(I)  Position  of  left-hand  boundary  of  Ith  element  (in).  The 

array  must  be  well  ordered  with  respect  to  I.  Ve  require 
XEL(I)  1  XEL(I-l).  Forward  elements  are  ordered  so  that 
the  first  is  to  the  left. 

HJL(l)  Mass  of  Ith  element  (lbm).  If  <  10-10  element  is 

interpreted  as  a  space,  MEL(I)  must  be  greater  than  10~1#. 


115 


FEL(I) 


Initial  resistance  to  motion  of  Ith  element  (lbf). 


NTYPEd)  If  NTYIE(I)  =  0,  element  is  treated  as  perfectly  plastic 

(deformation  under  loading  only) . 

If  NTYFE(I)  =  1,  element  is  treated  as  elastic. 

If  NTYPE(I)  =  2,  element  is  treated  as  rigid. 

If  NTYFE(I)  =  3,  element  is  treated  as  incompressible. 

0  i  NTYPE(I)  ^  3. 

NDATAd)  Number  of  pairs  of  entries  in  stress-strain  tatle  of 

Ith  element.  0  <  NDATA(I )  <  10. 


File  [C19]  "Constitutive  Data  for  Ith  Forward  Element*  (8F10.0) 

One  to  Three  Cards 


N.B.  File  required  if  and  only  if  MEL(I)  >  10-10  and  NTYPEd)  <  2.  (See 
File  [C18].)  File  is  repeated  for  each  element  for  which  NDATA  >  0. 
Files  [C19]  ,  as  a  group,  follow  Files  [Cl. 8]  ,  as  a  group. 

YEL(I,1)*  First  engineering  strain  taken  positive  in  compression 

(dimensionless).  Must  be  zero. 

RESEL(I,1)**  Corresponding  stress  taken  positive  in  compression  (psi) . 


RES EL (I, NDATA(I ))  Last  engineering  stress.  Should  exceed  maximum  pressure 

in  gun. 


The  array  YEL  must  be  well  ordered.  All  entries  must  be  in  the 
interval  [0,1] . 

The  array  RESEL  must  have  non-zero  entries  and  must  be  non-decreasing 
for  each  element. 
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File  [C19.1]  "Ith  Rear  Compactible  Filler  Eleaient  Properties" 

(3F10 .0,215)  One  Card 

N.B.  File  required  if  and  only  if  NBL  jfc  0.  (See  File  [M4].)  File  it 
repeated  NBL  times. 

XBL(I)  Position  of  right-hand  boundary  of  Ith  element  (in).  The 

array  must  be  well  ordered  with  respect  to  I.  XBL(I) 
observes  the  opposite  convention  from  that  of  XEL(I),  that 
is  XBL(I)  1  XBL(I-l)  for  all  I.  Rear  elements  are 
tabulated  so  that  the  first  is  to  the  right. 

MBL(I)  Mass  of  Ith  element  (lbo).  If  <  10_l®  element  is 

interpreted  as  a  space.  MBL(I)  must  be  greater  than  10'10 

FBL(I)  Initial  resistance  to  motion  of  Ith  element  (lbf). 

NTYPB(I)  If  NTYPB(I)  -  0,  element  is  treated  as  perfectly  plastic 

(deformation  under  loading  only)  . 

If  NTYPB(I)  =  1,  element  is  treated  as  elastic. 

If  NTTfPB(I)  =  2,  element  is  treated  as  rigid. 

If  NTYPB(I)  -  3,  element  is  treated  as  incompressible. 

0  <  mPB(I)  i  3. 

NDATB(I)  Number  of  pairs  of  entries  in  stress-strain  table  of 

Ith  element.  0  <  NDATB(I)  <  10. 


File  [C19.2]  'Constitutive  Data  for  Ith  Rear  Element"  (8F10.0) 

One  to  Three  Cards 

N.B.  File  required  if  and  only  if  »EL(I)  >  10_1°  and  NTYPE(I)  <2.  (See 

File  [C19.ll.)  File  is  repeated  for  each  element  for  which  N0ATB  >  0. 
Files  [C19.2] ,  as  a  group,  follow  Files  [C19.1],  as  a  group. 

YBL(I.l)*  First  engineering  strain  taken  positive  in  compression 

(dimensionless).  Must  be  zero. 

RESBL(I,1)**  Corresponding  stress  taken  positive  in  compression  (psi). 


RESBL(I,  NDATB(I))  Last  engineering  stress.  Should  exceed  maximum  pressure 

in  gun. 

*  The  array  YBL  must  be  well  ordered.  All  entries  must  be  in  the 
interval  [0,1]  . 

The  array  RESBL  must  have  nuirzero  entries  end  must  be  scn-dc creasing 
for  each  element. 
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File  [C23]  "Positions  for  Pressure  Table  Storage"  (8F10.0)  One  Card 

N.B.  File  required  if  and  only  if  KZPT  ^  0.  (See  File  [M4]  .) 

ZFT(l)  First  position,  relative  to  breecb  (in). 

« 

« 

« 

ZFT(NZFT)*  Last  position. 

*  Output  consists  of  tables  of  pressure  at  each  of  the  NZPT  positions 
together  with  value  less  that  at  the  NZPTth  position. 

File  [C23.01]  "Projectile  Afterbody  Pressure  History  Option"  (815) 

One  Card 

N.B.  File  required  if  and  only  if  NZPT  jt  0  and  MODET  =  1.  (See  File  [M4].) 

IZFT(  1)  0  if  ZPT(l)  is  fixed  in  the  tube  reference  frame. 

.  1  if  ZPT(l)  is  fixed  in  the  projectile  afterbody 

.  reference  frame. 

• 

IZPT(NZPr) 

File  [C23.1]  "Kinetics  Option  Counters"  (315)  One  Card 

N.B.  File  required  if  and  only  if  KMODE  ^  0.  (See  File  [M2].) 

NSPEC  Number  of  species.  1  i  NSFEC  1  10. 

NGASR  Number  of  reactions  occurring  in  the  mixture  of  combustion 

products.  0  1  NGASR  i  10.  If  NGASR  >  0,  File  [C23.7]  is 
required. 

N39LR  Number  of  reactions  occurring  in  the  s '1 id-propell ant . 

0  <  NS0LR  <  10.  If  NSCLR  >  0,  File  [C23.8]  is  required. 
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File  [C23.2] 


"Properties  of  Species 


(2A4.1X.A1.6F10.0) 


NS PEC  Cards 


N.B.  File  required  if  and  only  if  KMODE  0.  (See  File  [M2].) 

SPCNAM(I)  Name  of  species,  up  to  8  alphanumeric  characters. 

FAZE(l)  Phase  of  species:  G  =  gas  or  vapor*  L  -  liquid*  S  =  solid. 

One  alphanumeric  character. 

SPCV(I)  Specific  heat  at  constant  volume  ( lbf- in/ lto-R) . 

SPCP(I)  Specific  heat  at  constant  pressure  ( lbf-in/ lhn-R) . 

SIBV(I)  Covolume  (in*/lbm).  Only  required  if  FAZE  =  G. 

SPHOL(I)  Molecular  weight  ( lbei/  Ibmol)  .  Only  required  if  FAZE  =  G. 

SPDEN(I)  Density  ( lbm/ in* ) .  Only  required  if  FAZE  =  L  or  S. 

GLOS(I)  Transfer  coefficient  in  correlation  for  rate  of  deposition 

of  condensed  species  on  surface  of  solid  propellant 
(sec/in1).  Only  required  if  FAZE  =  L  or  S.  If  CL06  =  0, 
the  rate  of  deposition  is  zero. 

N.B.  The  following  File  [C23.3]  is  repeated  NPROP  times,  once  for  each  type 

of  solid  propellant. 

File  [C23.3]  "Composition  of  Propellant  Near  Field  Combustion  Products" 
(8F10.0)  NPROP  to  2*NPR0P  Cards 

N.B.  File  required  if  and  only  if  KMODE  £  0.  (See  File  [M2].) 

ECHO(I)  Chemical  energy  released  in  near  field  combustion 

( lbf- in/ lbm)  . 

70(1,1)  Mass  fraction  of  near  field  combustion  products 

.  corresponding  to  species  type  1  (-). 


I0(I,NSPEC)  Mass  fraction  of  near  field  combustion  products 

corresponding  to  species  type  NSH£C  (-)  . 
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File  [C23 .4]  "Composition  of  Igniter  Discharge  Combustion  Products” 
(EF10.0)  One  to  TVo  Cards. 

N.B.  File  required  if  and  only  if  OODE  jfc  0.  (See  File  [M2].) 

ECHO(I)  Chemical  energy  released  by  combustion  of  igniter  material 

prior  to  injection  (lbf-in/lbm) . 

10(1,1)  Mass  fraction  of  igniter  products  corresponding  to  species 

.  type  1  (-). 


10(1, NS PEC)  Mass  fraction  of  igniter  products  corresponding  to  species 

type  NSFEC  (-)  . 


File  [C23 .5]  "Composition  of  Initial  Ambient  Gas"  (8F10.0) 

One  to  TVo  Cards 

N.B.  File  required  if  and  only  if  KMODE  £  0.  (See  File  [M21.) 

ECHO(I)  Internal  energy  (thermal  component)  of  ambient  gas  (or 

mixture  if  condensed  species  are  initially  present) 
(lbf-in/lbm) . 

Y0  (1,1)  Mass  fraction  of  species  1  (-) . 


Y0(I, NSFEC)  Mass  fraction  of  species  NSHSC  (-) . 
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N. B.  The  following  File  is  reposted  NFROP  tines,  once  for  each  type  of 
solid  propellant.  Present  storage  limits  only  support  the  case 
NFROP  =  1  with  RHODE  =2. 


File  (C23  .6] 

"Initial  Composition  of  Solid  Propellant”  (8F10.0) 

NPROP  to  2*NPROP  Cards 

N, P.  File  required  if  and  only  if  RHODE  =  2.  (See  File  [M2] .) 

ECHO(I) 

Inactive. 

Y0(I,1) 

* 

« 

Mass  fraction  of  solid  propellant  corresponding  to 
species  1  (-). 

« 

• 

YO(I.NSreC) 

Mass  fraction  of  solid  propellant  corresponding  to  species 
NSPEC  (-). 

Fi.'e  [C23.61] 

"Composition  of  Traveling  Charge  Near  Field  Combustion 
Products"  (8F10.0)  One  to  TSio  Cards 

N. B.  File  required  if  and  only  if  RHODE  t  0  and  NTC  =  1.  (See  File  [M2].) 

ECHTCO 

Chemical  energy  released  in  near  field  combustion  of 
traveling  charge  (lbf-in/lbm) . 

YTCO(I) 

• 

« 

Mass  fraction  of  species  1  (-) . 

« 

YTCO(NSPEC) 

Mass  fraction  of  species  NSH&C  (-)  . 

File  [C23  .  62] 

"Types  of  Reactions  in  Mixture  of  r  -bastion  Products" 

(1015)  One  Card 

N.B.  File  required  if  and  only  if  KMODE  t  0,  NGASR  t  0  and  NTC  =  1 .  (See 
Files  [M2]  and  [C23.1].) 

KRCTYP(l) 

If  0,  reaction  1  is  of  the  Arrhenius  type  and  is  described 
by  File  [C23.7], 

If  1,  reaction  1  is  of  the  pressure  dependent  type  and  is 
described  by  [C23 .71] . 

If  NTC  =  0,  RRCTYP(I)  defaults  to  zero  for  all  I. 

SECTTP(NGASR) 

Type  of  reaction  NGASR 
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File  [C23.63]  'Species  Thermal  Equilibration  Switches*  (1015)  One  Card 

N.B.  File  required  if  and  only  if  DIODE  t  0.  NGASR  #0  and  NTC  -  1. 

(See  Files  [M2]  and  [C23.1].) 

XTEQL(l)  If  0,  species  1  is  assumed  to  be  in  constant  thermal 

.  equilibrium  with  mixture  of  combustion  products. 

If  1,  species  1  is  assumed  to  be  thermally  insulated  from 
.  mixture  of  combustion  products.  If  NTC  *=  0,  ETBOL(I) 

.  defaults  to  zero  for  all  I. 


KTBQL(NSFEC)  Thermal  equilibration  switch  for  species  NSHJC. 


File  [C23.7]  'Description  of  Arrhenius  Reactions  in  Mixture  of 
Combustion  Products' 

(815/ 8F10 .0/F10 .0,D10 .4 ,F10 .0,D10 .4 .4F10 .0/F10 .0) 

Four  Cards  For  Each  Such  Reaction 

N.B.  File  required  if  end  only  if  DIODE  *  0,  NGASR  r  0  and  DtCTIP(I)  «  0. 
(See  Files  [M2],  [C23.1]  and  [C23.62].) 

KRCNTB(  1,1)*  Pointer  to  first  species  acting  as  a  reactant  in 
.  reaction  I.  0  1  KRCN1B  *  NSPEC. 


KRCNTB(4,I)  Pointer  to  fourth  species  acting  as  a  reactant  in 

reaction  I. 

KPRODBd.I)*  Pointer  to  first  species  acting  as  a  product  in  reaction  I. 
0  <  KFRODB  i  NSPEC. 


EPR0DB(4,I)  Pointer  to  fourth  species  acting  as  a  product  in 

reaction  I. 

STOIBd.I)  Stoichiometric  coefficient  corresponding  to  first  reactant 

.  species  (-) .  Starts  a  new  card. 


STO1B(4,I)  Stoichiometric  coefficient  corresponding  to  fourth  reactant 

species  (-). 
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ST0IB(5,I)  Stoichiometric  coefficient  corresponding  to  first  product 

.  species  (-). 


STOIB(8,I)  Stoichiometric  coefficient  corresponding  to  fourth  product 

species  (-) . 

ECHB(I)  Chemical  energy  released  by  reaction  ( lbf-in/lbm) .  Starts 

a  new  card. 


ARCB(I) 


Pre-exponential  factor  in  Arrhenius  rate  law  (units  yield 
lbm/  in' -sec) . 


ARXB(I)  Temperature  exponent  in  Arrhenius  rate  law  (-) . 

AREB(I)  Activation  energy  in  Arrhenius  rate  law  ( Tbf-in/lbool; „ 


ARCB(1,I) 


Order  of  reaction  with  respect  to  concentration  of  first 
reactant  species  (-). 


AR0B(4,I)  Order  of  reaction  with  respect  to  concentration  of  fourth 

reactant  species  (-) . 

ARCB(5,I)  Order  of  reaction  with  respect  to  concentration  of  gas- 

phase  (-).  Starts  a  new  card. 


The  most  general  reaction  supported  involves  four  reactant  and  four 
product  species.  At  least  one  reactant  pointer  and  one  product 
pointer  must  be  different  from  zero.  A  zero  entry  simply  implies  a 
reduction  in  generality  of  the  reaction. 


123 


File  [C23.7-1]  "Description  of  Pressure-Dependent  Reactions  in  Mixtore  of 
Combustion  Products"  (8I5/8F10 .0/8F10 .0) 

Three  Cards  For  Each  Such  Reaction 

N.B.  File  required  if  and  only  if  ERODE  /  0.  NGASR  £  0  and  KRCTYP(I)  =  1. 
(See  Files  [M2],  [C23.1]  and  [C23.62].) 

KRCNTB(1,I)*  Pointer  to  first  species  acting  as  a  reactant  in 
.  reaction  I.  0  i  KRCNTB  £  NSKC. 


EROJTB(4,I)  Pointer  to  fourth  species  acting  as  a  reactant  in 

reaction  I. 

KPRODB(l.I)*  Pointer  to  first  species  acting  as  a  product  in  reaction  I. 
0  <  KPRODB  i  NS  PEC. 


EPR0DB(4.I)  Pointer  to  fourth  species  acting  as  a  product  in 

reaction  I. 

STDIB( 1,1)  Stoichiome tri 3  coefficient  corresponding  to  first  reactant 

.  species  (-) .  btart»  a  new  card. 


ST0IB(4,I)  Stoichiometric  coefficient  corresponding  to  fourth  reactant 

species  (-). 

STOIB(5,I)  Stoichiometric  coefficient  corresponding  to  first  product 

.  species  (-) . 


ST0IB(8,I)  Stoichiometric  coefficient  corresponding  to  fourth  product 

species  (-). 

ECHB(I)  Chemical  energy  released  by  reaction  (lbf-in/lbm) .  Starts 

a  new  card. 
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D1AB(I) 

Par ticl e 

diameter. 

assumed  constant  (in). 

B1B(I) 

Burn  rate 

additive 

constant  (in/sec). 

B2B(  I ) 

Burn  rate 

pre-exponential  factor  ( in/se c-psi^*®) 

BNB(I) 

Burn  rate 

exponent 

(-) 

The  most  general  reaction  supported  involves  four  reactant  and  four 
product  species.  At  least  one  reactant  pointer  and  one  product 
pointer  must  be  different  from  zero.  A  zero  entry  simply  implies  a 
reduction  in  generality  of  the  reaction. 


File  [C23.8]  "Description  of  Reactions  in  Solid  Propellant” 

( 815/  8F10 .0/F10.0,  D10.4 ,F10 .0, DIO .4 ,4F10 .0)  3*NS0LR  Ckrds 


N.B,  File  required  if  and  only  if  RMODE  ^  0  and  NSGLR  b  0.  (See  Files  [M2] 
and  [C23.1].) 

KRCNTS(l.I)*  Pointer  to  first  species  acting  as  a  reactant  in 
.  reaction  I,  0  <  KRCNTS  <  NSFEC. 


KRCMTS(4,I)  Pointer  to  fourth  species  acting  as  a  reactant  in 

reaction  I. 

KPRODS(l.I)*  Pointer  to  first  species  acting  as  a  product  in  reaction  I. 
0  <  KPRODS  <  NS  PEC. 


KPRODS (4,1)  Pointer  to  fourth  species  acting  as  a  product  in 

reaction  I. 

STOIS(l.I)  Stoichiometric  coefficient  corresponding  to  first  reactant 

species  (-) .  Starts  a  new  card. 


STOIS( 4,1) 


Stoichiometric  coefficient  corresponding  to  fourth  reactant 
species  (-). 
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ST0IS(5,I)  Stoi chi one  trie  coefficient  corresponding  to  first  product 

,  species  {-) . 


ST0IS(8,I)  Stoichiometric  coefficient  corresponding  to  fourth  product 

species  (-). 

ECHS(I)  Chemical  energy  relessed  by  reaction  (lbf-in/lbm)  .  Starts 

a  new  card. 


ittCS(I) 

ARXS(I) 

ARES(I) 


Pre-exponential  factor  in  Arrhenins  rate  law  (units  yield 
lb«/in*-sec) . 

Temperature  exponent  in  Arrhenius  rate  law  (-). 

Activation  energy  in  Arrhenius  rate  law  (lbf-in/lbmol) . 


AROS(l,I) 


Order  of  reaction  with  respect  to  concentration  of  first 
reactant  species  (-)  . 


AR0S(4,I) 


Order  of  reaction  with  respect  to  concentration  of  fourth 
reactant  species  (-). 


The  most  general  reaction  supported  inwolves  four  reactant  and  four 
product  species.  At  least  one  reactant  pointer  and  one  product 
pointer  must  be  different  from  zero.  A  zero  entry  simply  implies  a 
reduction  in  generality  of  the  reaction. 
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File  [C23  -9]  "Solid  Propellent  Thermal  Response  Parameters” 

(3I5.5F10.0)  One  Card 

N.B.  File  required  if  and  only  if  KMODE  «  2.  (See  File  [M2].) 

ICC  Number  of  stations  for  analysis  of  thermal  response  of 

solid  propellant.  0  i  ICC  1  21. 

ICQBC  If  0,  pyrolysis  lav  is  assumed  to  govei  i  regression  of 

surface  of  solid  propellant. 

If  1,  evaporative  (Clausius  Clapeyron)  lav  is  assumed  to 
govern  regression  of  surface  of  solid  propellant. 

NOFLAM  If  0,  heat  feedback  from  the  gas-phase  is  considered 

according  to  the  Zel’dovitch  formulism. 

If  1,  heat  feedback  from  the  gas-phase  is  neglected. 

AS  Pre-exponential  factor  in  solid  propellant  surface 

regression  lav.  If  ICCBC  =  0,  units  are  in/sec.  If 
NZCBC  ■  1,  units  are  psi. 

ES  Activation  energy  (NZCBC  =  0)  or  heat  of  vaporization 

(NZCBC  i=  1)  in  solid  propellant  surface  regression  lav 
(lbf-in/lbmol). 

TSEN  Temperature  sensitivity  of  propellant  steady  state  burn 

rate  (l/R). 

TIRANS  Time  interval  over  vhioh  external  stimulus  drops  to  zero 

following  ignition  of  the  surface  of  the  propellant  (sec). 

PTRANS  Pressure  above  vhlch  thermal  analysis  is  discontinued  and 

steady  state  combustion  is  assumed  to  occur  (psi). 
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File 

[C24] 

'Parameters  to  Define  Mesh  in  Invariant  Embedding  Analysis* 
(F10. 0,315)  One  Card 

N.B. 

File  required  if  and  only  if  KMODE  =  2.  (See  File  [M2].) 

DELS* 

Mesh  spacing  in  group  closest  to  heated  surface  (sec1/*) 

DELI  >0. 

N” 

Number  of  intervals  per  group. 

MDELX 

Integer  multiple  by  which  mesh  spacing  increases  from  group 
to  group  as  we  move  away  from  heated  surface. 

NGRP 

Number  of  groups.  1  i  NGRP  i  10. 

• 

Refers 

heated 

to  computational  coordinate  x/>|  a  where  x  is  distance  from 
surface  and  a  is  thermal  diffnsivity. 

«• 

N*NGRP 

<  200 . 

File 

[C25] 

"Tank  Gun  Option  Control  Data*  (1615)  One  Card 

N.B.  File  required  if  and  only  if  MODET  *  1,  (See  File  [M2].) 


NAFT 

JIS(l) 


J IS( 2) 


Number  of  pairs  of  data  to  describe  projectile  afterbody. 
NAFT  may  be  aero.  Otherwise  2  £  NAFT  ^  10. 

Switch  to  define  level  of  modeling  of  surface  source  term 
attributed  to  tube  wall. 

JIS(l)  -  Model 

0  No  source  term 

1  Tabular  source  term 

2  Rate  of  source  determined  according 

to  ignition  and  combustion  submodels. 

Switch  to  define  level  of  modeling  of  surface  source  term 
attributed  to  centerline  of  tube. 

Switch  to  define  level  of  modeling  of  surface  source  term 
attributed  to  afterbody  of  projectile. 


JIS( 3) 


NENDL 


NPR.M 


NRCT 

NSBGS(l) 


Number  of  endwalls  defined  by  packaging  of  charge. 

0  <  NENDL  <  21. 

Number  of  pemeability  data  seta.  0  1  NPRM  1  10. 

Number  of  reactivity  data  sett.  0  1  NRCT  1  9. 

Number  of  segments  having  different  properties  in  tube  vail 
surface  source.  0  1  NSBGS(l)  i  3.  Default  value  is  1. 


NSEG  S(  2) 


Number  of  segments  in  centerline  surface  ixmrce. 


NSBGS(3)  Number  of  segments  in  afterbody  surface  source. 

KCTRL  If  0.  control  charge  not  present. 

If  1.  control  charge  is  present.  Files  [C28]  -  [C28.7] 
required. 

NZCRE  Number  of  data  to  describe  external  geometry  of  ballistic 

control  tube. 


INHBCR(l)  If  0,  control  charge  not  deterred. 

If  1,  control  charge  is  deterred. 


File  [C25.1]  "Geometry  of  Projectile  Afterbody"  (8F10.0) 
One  to  Three  Cards. 

N.B.  File  required  if  and  only  if  MODET  =  1  and  NAFT  >  2. 
(See  Files  [M2]  and  [C25]  .) 

ZAFT(l)*  First  axial  position  on  afterbody  (in). 

RAFT(l)  Corresponding  radius  of  afterbody  (in). 


ZAFT(NAFT)  Last  axial  position. 

RAFT(NAFT)  Corresponding  radius  of  afterbody. 


The  origin  of  the  coordinate  ZAFT  is  independent  of  that  used  to  describe 
the  gun  tube.  XNOVAT  internally  reconciles  the  coordinate  frames  by 
assuming  that  ZAFT(NAFT)  coincides  with  ZBPR.  (See  File  [M13] .) 


N.B.  Hie  following  filet,  [C25.2]  through  [C26.3]  are  repeated,  at  needed, 
at  a  group  for  each  of  the  three  types  of  reactive  layer t  whose 
presence  is  indicated  by  a  non-zero  value  of  JIS(I). 

(See  File  [C25].) 


File  [C25.2]  "Thickness  and  Density  Counters  for  Reactive  Layer  I"  (415) 
One  Card 


N.B.  File  required  if  and  only  if  MODET  =  1  and  JIS(I)  #  0. 
(See  Filet  [M2]  and  [C25].) 


NSTAC(I) 


Number  of  data  to  describe  thickness  of  reactive  layer  I. 
2  <  NSTAC(I)  i  10. 


NSHOS(I.l)  Number  of  data  to  describe  density  of  first  segment  of 
reactive  layer  I  as  a  function  of  pressure. 

1  <  NRHOS(I)  <  10. 


NH0S(I,2) 
NRH0S(  1 , 3 ) 


Number  of  data  for  second  segment  of  reactive  layer  I. 
Number  of  data  for  third  segment. 


File  [C25.3]  "Thickness  of  Reactive  Layer  I”  (8F10.C)  One  to  Three  Cards 

N.B.  File  required  if  and  only  if  MODET  *  1  and  JIS(I)  4  0. 

(See  Files  [M2]  and  [C25].) 

ZAC(l.I)*  First  axial  position  (in). 

IBC(l.I)  Corresponding  thickness  of  reactive  layer  (in). 


ZAC(NSTAC(I) ,1)  Last  axial  position. 
THC(NSTAC(I ) , I)  Corresponding  thickness. 


Then  I  *  1  or  2,  ZAC  is  assumed  to  refer  to  the  coordinate  system  used  to 
define  the  tube  geometry.  Thru  1=3,  ZAC  is  assumed  to  correspond  to 
the  coordinate  system  used  to  define  the  afterbody  of  the  projectile. 

The  thickness  of  layer  I  is  taken  to  be  zero  outside  the  table  range. 
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File  [C25.31]  "Segment  Pointers  for  Reactive  Layer  I  (8IS)  One  Card 

N.B.  File  required  if  and  only  if  MODET  =  1.  JIS(I)  *  0  and  NSBBS(I)  >  1. 
( See  Files  [M2]  and  [C25]  .) 

NSG(1,I)  Segment  property  type  for  axial  positions  between  ZAC(l.I) 

and  ZAC(2,1) . 


NSG ( NSTAC ( I ) -1 , I)  Segment  property  type  for  axial  positions  between 

ZAC ( NSTAC (I ) -1 , I )  and  ZAC(NSTACU)  ,1) . 

NSG(NSTAC(I) ,1)  Not  used. 

N.B.  If  Reactive  Layer  I  consists  of  more  than  one  segment,  the  following 
Files  [C25.4]  through  [C26.3]  are  repeated  as  a  gronp  for  each  of  the 
segments  of  the  layer.  The  snbscript  pertaining  to  the  segment  is 
suppressed  in  the  following  discussion. 


File  [C25 .4]  "Density  of  Reactive  Layer  I”  (8F10.0)  One  to  Three  Cards 

N.B.  File  required  if  and  only  if  MDDET  -  1  and  JIS(I)  ^  0. 

(See  Files  [M2]  and  [C25] .) 

RH0S(1,I)*  Value  of  density  at  first  pressure  (lbai/ in* )  . 

PRH0S(1,I)  First  value  of  pressure  (psi) . 


jffl0S(NRH0S(I) , I)  Last  value  of  density. 
PBH0S(NRH0S(I) ,1)  Last  value  of  pressure. 


*  Outside  the  table  range,  the  first  or  the  last  value  of  RBOS  applies.  If 
a  single  value  is  specified,  the  density  is  taken  to  be  constant  and 
therefore  independent  of  pressure. 

At  present,  case  compressibility  is  only  modeled  if  JIS(I)  *  2.  (See  File 
[C25]  .) 
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File  [C25.5]  "Counter  to  Describe  Discbirge  of  Reactive  Layer  I”  (215) 
One  Card 

N.B.  File  required  if  and  only  if  MDDET  =  1  and  JIS(I)  =  1. 

(See  Files  [M2]  and  [C25].) 

IIS(I)  Number  of  axial  positions  in  discharge  table. 

2  <  IIS(I)  <  10. 

JJS(I)  Number  of  times  in  discharge  table.  2  i  JJS(I)  1  10. 


File  [C25.6]  "Positions  to  Describe  Discharge  of  Reactive  Layer  1* 
(8F10.0)  One  to  Two  Cards 

N.B.  File  required  if  and  only  if  MODET  -  1  and  JIS(I)  =  1. 

(See  Files  [M2]  and  [C25] .) 

ZPHIS(1,I)*  First  position  (in). 


ZPHIS(IIS(I) ,1)  Last  position  (in). 


The  coordinate  frame  for  ZPHIS  is  assumed  to  accord  with  that  for  ZAC. 
(See  File  [C25  .3]  .) 


File  [C25.7]  "Times  to  Describe  Discharge  of  Reactive  Layer  I"  (8F10.0) 
One  to  Two  Cards 

N.B.  File  required  if  and  only  if  MDDET  *  1  and  JIS(I)  =1. 

(See  Files  [M2]  aud  [C25]  .) 

TPBIS(l.I)  First  time  (msec). 


TFHIS(  JJS(I ) ,  I)  Last  time  (msec). 
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File  [C25.8]  "Rate  of  Discharge  of  Reactive  Layer  I"  (8F10.0) 

One  to  Thirteen  Cards 

N.  B.  File  required  if  and  only  if  SDDET  =  1  and  JIS(I)  =  1. 

(See  Files  [M2]  and  [C25]  .) 

PHIS(l.l.I)  Rate  of  discharge  at  first  position  and  first  time 
(lbm/ in-  sec) . 

FHIS(2,1,I)  Rate  of  discharge  at  second  position  and  first  time. 

FHISUIS(I)  .1,1)  Rate  of  discharge  at  last  position  and  first  time. 

a 

FHIS(  IIS(I) ,  JJS(I)  ,1)  Rate  of  discharge  at  last  position  and  last  time. 

File  [C25.9]  "Burn  Rate  Counters  for  Reactive  Layer  I"  (215)  One  Card 

N.B.  File  required  if  and  only  if  MDDET  =  1  and  JIS(I)  *  2. 

(See  Files  [M2]  and  [C25].) 

KBRDS(I)  Number  of  burn  rate  data.  1  1  EBRDS  1  10. 

KNHIB(I)  If  0,  layer  I  is  not  deterred. 

If  1,  layer  I  is  deterred. 

File  [C26]  "Burn  Rate  Data  for  Reactive  Layer  I"  (8F10.0) 

One  to  Three  Cards 

N.B.  File  required  if  and  only  if  MODE!  >=  1  and  JIS(I)  -  2. 

(See  Files  [M2]  and  [C25]  .) 

DPPRS(l.I)*  First  pressure  limit  (psi). 

B22S(1,I)  Value  of  pre-exponent  for  pressures  less  than  DPPRS(  1,1) 

( in/sec-p*iBNNs) . 


BNNS( 1, I) 


Value  of  exponent  for  pressure*  less  than  OPPRS( 1, I)  (-) . 


* 

OPPRS(KBRDSd) ,  I) *  Last  pressure  limit  (psi)  . 

B22S(KBRDS(I) , I)  Value  of  pre-exponent  for  pressures  less  than 

UPPRS(EBRDS(I) , I)  but  greater  that  UPPRS(KBRDSd)-l.I) . 

BNNS(KBRDSd) ,  I)  Corresponding  exponent  (-)  . 

Outside  the  table  range  the  first  and  last  values  are  used. 


File  [C26.1]  "Ignition  Data  for  Reactive  Layer  I"  (8F10.0)  One  Card 

N. B.  File  required  if  and  only  if  MODET  =  1  and  JIS(I)  =  2. 

(See  Files  [M2]  and  [C253  .) 

BISd)  Burn  rate  additive  constant  (in/  sec). 

TMPIBS(I)  Ignition  temperature  (R) . 

KPSd)  Thermal  conductivity  (lbf/sec-R). 

ALPBAS(I)  Thermal  diffusivity  (in*/sec). 


File  [C26.2]  "Thennochemical  Data  for  Reactive  Layer  I"  (8F10.0)  One  Card 

N.B.  File  required  if  and  only  if  MODET  =  1  and  JIS(I)  -  2. 

(See  Files  [M2]  and  [C25].) 

EIGS(I)  Chemical  energy  released  during  combustion  (lbf-in/lbm)  . 

GMS(I)  Molecular  weight  of  products  of  combustion  (lbm/lbmol). 

GAMAS d)  Ratio  of  specific  heats  (-)  . 
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File  [C26.21]  "Composition  of  Reactive  Layer  Near  Field  Combustion 
Product s"  (8F10.0)  One  Card 

N.  B.  File  required  if  and  only  if  MODET  =  1,  KM)OE  =  1  and  JIS(I)  =  2. 
(See  Files  [M2]  and  [C25] .) 

YSO(I.l)  Mass  fraction  of  species  1  (-) . 


« 

YS0( I, NSFEC)  Mass  fraction  of  species  NSPEC  (-) . 


File  [C26.3]  'Properties  of  Deterred  Layer  of  Reactive  Layer  I"  (8F10.0) 
One  Card 

N.B.  File  required  if  and  only  if  MODET  *  1,  JIS(I)  =  2  and  INHIB(I)  *=  1. 

(See  Files  [M2],  [C25]  and  [C25.9].) 

ECHIS(I)*  Internal  energy  released  in  combustion  at  start  of  deterred 
layer  ( lbf- in/ lbm) . 

EGBIS2(1)  Internal  energy  released  in  combustion  at  end  of  deterred  layer 
(lbf-in/lbm) . 

HGFAS(I) *  Factor  by  which  burn  rate  is  multiplied  at  start  of  deterred 
layer  (-)  . 

KGFAS2(I)  Factor  by  which  burn  rate  is  multiplied  at  end  of  deterred 
1 ayer  (-) . 

HIBS(I)  Depth  of  deterred  layer  ( in)  <■ 

*  Values  within  deterred  layer  deduced  by  linear  spacewise  interpolation 
with  an  allowance  for  compression  of  the  reactive  layer.  Final  values 
need  not  be  the  same  as  those  of  the  undeterred  layer. 


N.B.  File  required  if  end  only  if  MODET  =  1  end  NPRM  +  0. 
(See  Files  [M2]  end  [C25]  .) 


PRM(I)  Initial  flow  resistance  coefficient  for  I th  model  (-) , 

HDPSTR(l)  Piessure  differential  at  which  rupture  of  endwell  commences 
(psi) . 

HDPINT(l)  Tine  interval  over  which  rupture  of  enthrall  is  completed 
(msec) , 
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N.B.  Hie  following  Files,  [C26.6]  through  [C27],  are  repeated  as  a  group 
for  each  of  the  NRCT  reactivity  models. 


File  [C26.6]  "File  Counters  for  I- th  Reactivity  Xbdel*  (315)  One  Card 


N.B.  File  required  if  and  only  if  M3DET  =  1  and  NRCT  P  0, 

(See  Files  [M2]  and  [C25].) 

KBRDE(I)  Number  of  data  ia  modeled  burn  rate  description  (-) .  If 

KBRDE(I)  =  0,  it  is  assumed  that  tabular  discharge  data  are 
specified.  0  1  KBRDE(I)  £  10. 

IGCRIT(I)  Ignition  criterion  for  modeled  burn  rate  description. 

Ignition  is  determined  by  reference  to  VALIG(I) 

(File  [C26.8])  as  follows. 

IGCRIT(I)  VALIG(I) 

1  Time  Delay  (w sec) 

2  Ambient  gas  temperature  (R) 

Z  Neighboring  propellant  temperature  (R) 


JRCT(I)  Number  of  data  in  tabular  burn  rate  description  (-)  . 

JRCT(I)  =  0  or  2  £  JRCT(I)  £  8. 


File  [C26.7]  "Thermochemical  Data  for  I-th  Reactivity  Model"  (8F10.0) 
One  Card 

N.B.  File  required  if  and  only  if  M9DET  =  1  and  NRCT  P  0. 

(See  Files  [M2]  and  [C25] .) 

RHOE(I)  Density  ( lbm/ in* ) . 

EIGE  (I)  Chemical  energy  released  during  combustion  (lbf-in/lbo) . 

G*E(I)  Molecular  weight  of  combustion  products  (lbm/ Ibmol) . 

GAMAE(I)  Ratio  of  specific  heats  of  combustion  products  (-) . 
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File  [C26.71]  'Composition  of  Near  Field  Combustion  Products  of  I-th 
Reactivity  Model'  (8F10.0)  One  Card 

N.B.  File  required  if  and  only  if  MODET  =  1,  MODE  =  1  and  NRCT  4  0. 
(See  Files  [M2]  and  [C25] .) 

TEE (1,1)  Mass  fraction  of  species  1  (-)  . 


lEE(I.NSPEC) 


Mass  fraction  of  species  NSPEC  (-) . 


File  [C26.8]  "Ignition  Value  for  I-th  Reactivity  Model"  (8F10.0)  One  Card 


N.B.  File  required  if  and  only  if  MODET  =  1,  NRCT  4  0  and  KBRDE(I)  4  0. 

(See  Files  [M2],  [C25]  and  [C26.6].) 

VALIG(I)  Value  of  time  delay  (msec),  gas  temperature  (R)  or  propellant 

temperature  (R)  in  accordance  with  IGCRIT(J)  as  defined  above, 
(See  File  [C26.6] .) 


File  [C26.9]  "Burn  Rate  Data  for  I-th  Reactivity  Model*  (8F10.0) 

One  to  Three  Cards 

N.B.  File  required  if  and  only  if  MODET  =  1,  NRCT  4  0  and  KBRDE(I)  4  0. 
(See  Files  [M2],  [C25]  and  [C26.6].) 

FPPRE(1,I)*  First  pressure  limit  (psi), 

B22E(1,I)  Value  of  pre-exponent  for  pressures  less  than  UPPRE(1, I) 

(in/.ec-p.iBNNE). 

BNNE(1,I)  Value  of  exponent  for  pressures  less  than  UPPSiE(l,I)  (-) . 

e 

UPPRE(KBRDE (I) ,1) *  Last  pressure  limit  (psi). 

B22E(KBRDE(I) ,  I)  Value  of  pre-exponent  for  pressures  less  than 

UPPRE(KBRDEU).I)  but  greater  than  UPPRE(IBRDE(I)-1 , 1) . 

BNNT(XBRDE(I ) , I)  Corresponding  exponent  (-)  . 

Outside  the  table  ranee  the  first  and  last  valuet  are  used. 
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File  [C27] 


"Tabular  Description  of  I-tb  Reactivity  Model"  (8F10.0) 
One  or  Two  Cards 


N.B.  File  required  if  and  only  if  JCDET  *  1,  NR CT  +  0,  KBRDB(I)  “  0  and 
2  <  JRCT(I)  <8.  (See  Files  [M2],  [C25]  and  [C26.«].) 

KRCr(l,I)  First  value  of  tine  (nsec). 

FL0RCI(1,I)  Corresponding  rate  of  combustion  (lbm/in*- sec) . 

♦ 

* 

* 

TRCT( JRCl'(I) ,  1)  Last  value  of  tine. 

FLORCT(JRCT(.l) ,  I)  Corresponding  rate  of  coabustion. 


File  [C28]  "Internal  Properties  of  Control  Charge  Coabustion  Chanber" 
(81710.0)  One  Card 

N.B.  File  required  if  and  only  if  1DDET  =  1  and  KCTRL  £  0. 

(See  Files  [M2]  and  [C2S]  .) 

CDCR  Dischatg*.  coefficient  for  venting  to  gun  chanber  (-) . 

RCRI  Internal  radius  of  bore  (in). 

VCRO  Initial  chamber  voluae  corresponding  to  zero  afterbody 

offset  (in*). 

ZCRO  Initial  afterbody  offset  (in). 
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File  [C28.1]  'External  Properties  cf  Control  Charge  Combustion  Qi  amber" 
(3F10.0)  NLCBE  Cards 

N.B.  File  required  if  and  only  if  MO!>ET  -  1  and  ICTRL  t  0. 

(See  Files  [M2]  and  [C251J 

ZCRE(l)  First  axial  position  relative  to  breechface  of  gnn  (in). 

RCKE(l)  Corresponding  external  radius  of  chamber  (in). 

AVEKT(l)  Total  sidewall  vent  area  exposed  when  base  of  afterbody  is  at 
ZCEE(l)  (in*). 

ZCRE(2)  Second  axial  position.  (New  Card) 

ZCRE(NZCX£)  Last  axial  position.  (New  Card) 

ICBE(NCCRE) 

AVENT<  NZCRE) 

File  [C28.23  "Control  Charge  Type*  (5M.2F10.0)  One  Card 

N,  B.  File  required  if  and  only  if  WDET  -  1  and  KCTRL  i*  0. 

(See  Files  [M2]  and  [C25].) 

6BNMCS  Nave  of  propellant.  Up  to  20  alphanumeric  characters. 

TOTCE  Mass  of  control  chaxge  (lb*). 

Density  of  propellaut  (lbn/in* ). 


SBOKR 


File  [C28.3] 

"Control  Charge  Form  Function"  (I5.4F10.0)  One 

Card 

N.B.  File  required  if  and  only  if  MODET  •=  1  and  KCTRL  f  O. 
(See  Files  [M2]  and  [C25] .) 

NFRJf~R 

Form  function  indicator.  See  diaousaion  of  File 
Allowable  value*  of  NFRMCR  are  1.  2,  5,  6  and  7. 

[K8]  . 

ODCR 

Grain  dimension  (in).  See  File  [M8] . 

DPRFCR 

Grain  dimension  (in).  See  File  [M8] . 

GLENCR 

Grain  length  (In) . 

NPRFCR 

Number  of  perforations  (-). 

File  [C28.4] 

"Control  Charge  Burn  Rate  Connter"  (15)  One  Card 

N.B.  File  required  if  end  only  if  MODET  -  i  and  KCTRL  *  0. 

(See  Filet  [M2]  and  [C25]  .) 

NBRDCR  Number  of  sett  of  values  used  to  describe  burn  rate. 

1<  NBRDCX  <  10. 


File  [C28.5]  "Control  Charge  Born  Kate  Description"  (8F10.0) 

One  to  Four  Cards 

N.B.  File  required  if  and  only  if  MODET  ■  1  and  ECTKL  +  0. 

(See  Files  [M2]  and  [C25]  .) 

UPPGR(l)  Maximum  value  of  pressure  for  corresponding  values  of  burn  rate 
pre-exponential  and  exponential  factors  (psi) 

B22CR(1)  Burning  rate  pre-ezponential  factor  ( in/aeo-psin) , 

BWCX(l)  Burning  rate  exponent  (-) . 

* 

• 

UPPCR(MJRDCR)  Maximum  value  for  last  set  of  burn  rate  data. 

B22(NBKDCK)  Corresponding  pre-exponent. 

BNN(NBRDCR)  Corresponding  exponent. 

B1CR  Burning  rate  additive  constant  (in/sec). 

DELCR  Ignition  delay  (msec). 
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file  [C28.6]  "Control  Charge  Thermochemistry"  (8F10.0)  One  Card 


l.  B.  File  reqoired  if  and  only  if  MODET  =  1  and  KCTRL  0. 

(See  Files  [M2]  and  [C25] .) 

CECR  Internal  energy  released  in  combnstion  (lbf-in/lbm). 

MCR  Molecular  weight  of  combustion  products  (lbm/lbmol). 

AMA.CR  Ratio  of  specific  heats  (-). 


file  [C28.7]  "Properties  of  Deterred  Layer"  (8F10.0)  One  Card 


LB.  File  required  if  and  only  if  M5DET  -  1  and  KCTRL  *  0  and  INHBCR  =  1. 

(See  Files  [M2]  and  [025].) 

£CRIB*  Internal  energy  released  in  combustion  at  start  of 

deterred  layer  (lbf-in/lbm). 

1CRIB2  Internal  energy  released  in  combustion  at  end  of 

deterred  layer  (lbf-in/lbm). 

RGFCR  Factor  by  which  bu*:n  rate  is  multiplied  at  start  of 

deterred  1 ayer  (-) . 

BGFCR2  Factor  by  which  burn  rate  is  multiplied  at  end  of 

deterred  1 ayer  (-) . 

9IBXCR  Depth  of  inhibited  layer  (in). 


Values  within  deterred  layer  deduced  by  linear  spacewise 
interpolation.  Final  values  need  not  be  the  same  as  those  of  the 
undeterred  propellant. 
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File  [TCI] 


"Traveling  Charge  Control  Data"  (715)  One  Card  TC-Xandatory 


IDEAL 


NPRC 


NPROP 


Propellant  burn  rate  indicator. 

0  -  Measured  burn  rate  data.  See  Files  [TC5]-[TC8]. 

1  -  Not  supported. 

2  -  Ideal  burning  with  prespecif  ied  value  of  pressure  on 

either  side  of  gas/ traveling  charge  interface  or  of 
projectile  acceleration.  Note  the  discussion  of  SIBMAX, 
MACH,  APMAX  in  File  [TC2] . 

0  -  Traveling  charge  treated  aa  rigid. 

1  -  Traveling  charge  treated  as  a  continues  with  analytical 

description  of  rheology.  File  [TC9]  required. 

2  -  Traveling  charge  treated  as  continuum  with  tabular 

description  of  rheology.  Files  [TC10]  and  [TCll] 
required. 

Number  of  traveling  charge  increments  (maximum  of  20). 


Mf  FR 


Mi  RES 1 


M1SES2 


NBRES3 


Propellant  Vail  Friction  Parameter. 

0  -  Friction  between  propellant  and  tube  not  considered. 

1  -  Friction  due  to  gas  film.  File  [TC12]  required. 

>0  -  Number  of  entries  in  velocity  dependent  coefficient  of 
friction  table  (maximum  of  10).  File  [TC13]  required. 

0  -  Obturator  resistance  not  given  as  table, 

>0  -  Number  of  entries  in  table  of  resistive  pressure  versus 
travel  (maximum  of  10).  File  [TC14]  required. 

0  -  Resistance  due  to  shocked  air  cot  considered. 

1  -  Resistance  due  to  shocked  air  considered.  File  [TC15] 
required. 

0  -  Obturator  resistance  not  proportional  to  setback 
pressure. 

>0  -  Number  of  entries  in  table  of  velocity-dependent 

coefficient  of  friction  of  obturator  (maximum  of  10). 
Files  [TC16]  and  [TC17J  required. 


File  [TC2]  "Mesh  Parameters"  (I5,F10,0)  One  Card  TC-Mandstory 

MAXDIM  Maximan  number  of  mesh  points  to  be  used  in  continuum 

representation  of  traveling  charge  (1  100). 

DXMIN  Minimum  mesh  size  for  continuum  representation  (in). 
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File  [TC3]  "General  Properties"  (6F10.0)  One  Card  TC-Mandatory 


DB 

XEB 

PRM 

SIGMAX 


MACE 


APMAX 


Diameter  of  tube  (in). 

Initial  length  of  gas  col  man  (in). 

Mass  of  projectile  (lbm). 

Maximum  value  of  pressure  at  gas/propcllant  interface  (psi). 
If  SIGMAX  =  0.  no  restriction  is  considered.  SIBMAX  pertains 
to  the  reacted  or  the  unreacted  aide  of  the  flame  according 
as  NFORS(KX)  =  0  or  1  respectively.  (See  File  [TC5].) 

Maximum  value  of  Mach  number  of  combustion  products  relative 
to  regressing  surface,  If  MACH  =  0,  no  restriction  is 
consider  td. 

Maximum  value  of  acceleration  of  projectile  (gravities).  If 
APMAX  =  0,  no  restriction  is  considered. 


N.B.  Files  ITC4]  through  [TC133  pertain  to  a  specific  type  of 

propellant.  The  sequence  [TC4]  through  [TC13] ,  subject  to 
relevant  contingencies,  is  repeated  for  each  type  of  traveling 
charge  propellant  in  the  problem.  The  index  KK  used  in  the 
following  file  descriptions  runs  through  the  values 
1,2,  .  .  .,  NPROP,  where  NPROP  is  defined  in  File  [TCI]. 


File  [TC4]  "Propellant  Thermo  chemical  Properties"  (8F10.0)  One  Card 

TC -Mandatory 


XuAM(KK) 

XHV(KK) 

XMOL(KK) 

XECHEM(KK) 

XRiOP(KK) 

XCM(KK) 

XCRIT(KK) 


Ratio  of  specific  heats  of  gas  (--) . 

Covolume  (inVibm). 

Molecular  weight  (Ibm/lbmol) , 

Chemical  energy  of  propellant  (ibf-in/lbm) . 

Density  of  solid  propellant  at  zero  pressure  (lbm/ in'). 

Mass  of  propellant  (lbm). 

Time  de)  ay  following  first  exposure  of  increment  base  before 
combustion  begins  (msec). 


XDEL(XK) 


Time  interval  over  which  increment  combustion  rate  increases 
to  steady-state  value  (msec). 


File  [TC5]  "Burn  Rate  Switches”  (215)  One  Card  TC-Mandatory 


MNJRKXK)  0  -  Exponential  form  for  neatnred  burn  rate  description  of 
propollant  XX. 

-  Tabular  form  for  measured  burn  rate  description  of 
propellant  KX. 

NPORS(KK)  0  -  Measured  burn  rate  data  are  given  as  function  of  pressure 
on  reacted  side  of  flame.  Ideal  burn  rate  value  SIGMAX 
(see  File  [TC3])  pertains  to  reacted  side. 

1  -  Measured  burn  rate  data  are  given  as  function  of  pressure 
on  unreacted  side  of  flame.  Ideal  burn  rate  value  SIGMAX 
pertains  to  unreacted  side. 


File  [TC6j  " Tabular  Burn  Rate  Dpta”  (2I5.F10.0)  One  Card  TC-Contingent 

N.  B.  File  is  required  if  and  only  if  IDEAL  «*  0  and  MN)Rl(KX)  =  1, 

(See  Files  [TCI]  and  [TCS] )  . 

NBRF(XK)  Number  of  pairs  of  data  in  KX-th  tabular  description  of  burn 
rate  (maximum  of  20.) 

INTD  Number  of  table  entries  to  use  in  divided  difference 

interpolation.  The  same  value  is  assumed  for  all 
propellants. 

Bl(KX)  Burn  rate  additive  constant  for  KX-th  increment  (in/ sec). 


File  [TC7]  "Tabular  Burn  Rate  Data  (cont'd)"  (2F10.0)  J®RP(SX)  Cards 

TC-Contingent 

N. B.  File  is  required  if  and  only  if  IDEAL  =  0  and  MNBRl(KK)  -  1. 

(See  Files  [TCI]  and  [TC5]  .) 


PE ( 1, XX)  First  value  of  pressure  in  tsbular  description  of  burn  rate 

of  propellant  KX  (psi). 

RP(1,KK)  First  value  of  burn  rate  of  propellant  XX  corresponding  to 
re(l, XX)  ( in/sec) . 


re( NBRP(KK) , EX) 
RP(  N3RP(KX) , XX) 
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?ile  [TC8]  "Exponential  Burn  Rate  Data"  (3F10.0)  One  Card  TC-Contingent 


'J. B.  File  required  if  and  only  if  IDEAL  =  0  and  HNBRl(EE)  -  0. 

(See  Files  [TCI]  and  [TCS]  .) 

Jl(KE)  Burn  rate  additive  constant  for  EE-th  propellant  (in/sec). 

)2(EE)  Burn  rate  pre-exponential  coefficient  for  EE-th  propellant 

(  in/  se  c-psi®N)  , 

)N(EK)  Burn  rate  exponent  for  EE-th  propellant  (  ■)  . 


Pile  [TC9]  "Propellant  Analytical  Rheology  Data"  (2F10.0)  One  Card 

TC-Contingent 


N.B.  File  required  if  and  only  if  NHtC  -  _  (See  File  [TCI].) 

EADP(EK)  Compressive  wave  speed  at  zero  pressure  in  analytical 
description  of  propellant  rheology  (in/sec). 

XADVN(EE)  Expansion  wave  speed  in  analytical  description  of  propellant 
rheology  (in/ sec).  If  XAIRrN(EE)  is  entered  so  that  it  is 
less  than  the  nominal  compressive  wave  speed,  the  loading 
value  is  used.  By  entering  XADWN(KE)  “  0,  a  reversible  law 
is  defined. 


File  [TC10]  "Propellant  Tabular  Rheology  Data"  (15)  One  Card 

TC-Continge nt 


N.B.  File  required  if  and  only  if  NFRC  *  2.  (See  File  [TCI].) 

MNSS(EE)  Number  of  pairs  of  entries  in  tabular  description  of 
propellant  rheology  (maximum  of  20). 
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File  [TC11]  'Propellant  Tabular  Rheology  Data  (cont'd)'  (3F10.0) 

MNSS(KK)  Cards  TC-Contingent 


N.B.  File  required  if  and  only  if  NFRC  =  2.  (See  File  [TCI].) 

PSTA(1,K)  First  value  of  percent  strain  in  tabular  description  of 
propellant  rheology  (-). 

STRL(1,KK)  Corresponding  value  of  pressure  on  nominal  loading 
(compression  curve  (psi)) . 

STRU(l.KK)  Corresponding  value  of  pressure  on  nominal  unloading 
.  (expansion  curve  (psi)). 


P$TA(HNSS(EK) , EE) 
STRL(MNSJ(KK),KK) 
STRD(KNSS(KK) , KK) 


File  [TC12]  'Ablative  Film  Data”  (2F10.0)  One  Card  TC-Contingent 

N.B.  File  required  if  and  only  if  MfFR  <  0.  (See  File  [TCI].) 

VISLYR(KK)  Viscosity  of  gas  film  used  to  lubricate  propellant 
( lbm/  in-se  c)  . 

DELYR(EK)  Thickness  of  film  (in). 


File  [TC13]  'Propellant  Friction  Data'  (8F10.0)  One  to  Three  Cards 

TC-Coatinge  nt 

N.B.  File  required  if  and  only  if  M(FR  >  0.  (See  File  [TCI].) 

All  propellants  are  assumed  to  have  the  same  number  of  data  NWFR. 

AMDVd.KX)  First  value  of  velocity  of  propellant  (in/sec). 

AMU(l.KK)  Corresponding  coefficient  of  friction  on  tube  (-) . 


AMUV(  NffFR,  KK)  Last  value  of  velocity  (in/sec). 
AHU(MfFR,KK)  Corres ponding  coefficient  of  friction  (-). 
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File  [TCJ4]  "Tabular  Bore  Resistance  Data"  (8F10.0)  One  to  Three  Cards 

TC-Continge  nt 

N.B.  File  required  if  and  only  if  MJRES1  fc  0.  (See  File  [TCI].) 

BRX(l)  First  value  of  projectile  travel  (in). 

BR(1)  Corresponding  value  of  resistive  pressure  due  to  obturator 

(psi) . 

e 

a 

BRI(MlRESl)  Last  value  of  projectile  travel  (in). 

BR(NBRESl)  Corresponding  value  of  resistive  pressure  (psi). 

File  [TC15]  "Barrel  Shock  Resistance  Data"  (4F10.0)  One  Card 

TC-Contingent 

N.B.  File  required  if  and  only  if  NBRES2  t  0.  (See  File  [TCI].) 

A1RGAM  Ratio  of  specific  heats  of  air  (-)  . 

AIRPO  Pressure  of  air  in  barrel  (psi). 

AIkTO  Temperature  of  air  in  barrel  (R) . 

AIRMW  Molecular  weight  of  air  in  barrel  ( lbo/ lta>  ol )  . 

File  [TC16]  "Obturator  Setback  Resistance  Data”  (3F10.0)  One  Card 

TC-Contingent 

N.B.  File  required  if  and  only  if  NBRES3  0.  (See  File  [TCI].) 

FRJSI  Mass  of  projectile  ahead  of  midpoint  of  obturating  band 

( lbm)  . 

ELd  Length  of  bearing  section  of  obturating  band  (in). 

ANU  Poisson's  ratio  of  obturating  band  (-) . 


148 


File  [TC17]  ’Obturator  Friction  Data*  (8F10.0)  One  to  Three  Cards 

TC-Cont  intent 

N.  B.  File  required  if  and  only  if  NBBES3  f  0.  (See  File  [TCI].) 

BMJV(l)  First  value  of  velocity  of  projectile  (in/sec). 

BM0(1)  Corresponding  value  of  coefficient  of  friction  between 

-obturator  and  tube  (-). 


BMTV(N3BES3 )  Last  value  of  velocity  of  projectile  (in/sec). 
BMJ(NBRES3)  Corresponding  coefficient  of  friction  (-) . 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 
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